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TECHNICAL  REPORT  SUMMARY 

Geophysical  systems  consisting  of  electrical  and  acoustical  tools 
were  developed  and  special  mapping  techniques  applied  to  obtain 
geotechnical  data  in  Monterey  Bay,  California:  The  objective  of 
the  field  exercise  was  to  try  and  relate  geophysical  measurements 
to  the  mass  physical  properties  and  engineering  properties  of  sea¬ 
floor  sediments.  Measurements  consisted  of;  (1)  mapping  the  Ray¬ 
leigh  Reflection  Coefficient  continuously  by  a  ship  underway,  (2) 
measuring  in  situ,  the  propagation  velocities  of  compress ional  and 
transverse  (shear)  waves,  and  (3)  taking  both  underway  and  in  situ 
measurements  of  sediment  conductivity. 

P  and  S  wave  velocities  were  measured  using  a  prototype  seismic 
shear  wave  generator.  Measurements  of  sediment  conductivity  were 
taken  with  a  direct  current  resistivity  system  deployed  in  a  sur¬ 
face-towed  configuration  and  as  a  diver-operated  expandable  array 
implanted  in  the  seafloor.  These  data,  along  with  the  computed  re¬ 
flection  coefficients  (bottom  loss)  data,  were  correlated  to  analy¬ 
sis  of  samples  taken  using  a  unique  diver-operated  core  sampler  and 
vibracorer  sampler  that  were  specially  designed  for  the  field  exper¬ 
iment. 

The  surveys  and  sampling  programs  were  carried  out  in  Santa  Cruz 
Harbor,  Monterey  Bay,  California.  Surveys  were  conducted  on  a 
fine-grained  scale.  Grid  patterns  for  profiling  and  sampling  cov¬ 
ered  an  area  of  30  square  miles.  A  high-accuracy  precision  radio¬ 
location  navigation  and  positioning  system  was  employed  throughout 
the  field  exercise. 

The  experiments,  equipment  employed,  and  corroboration  of  data  to 
previous  studies  of  a  similar  nature  are  described.  It  was  con¬ 
cluded  that  each  of  the  systems  employed  have  potential  application 
for  measuring  remotely  and  directly  {in  situ)  parameters  of  mass 
physical  properties  of  seafloor  sediments,  and  that  only  the  shear 
wave  system  showed  potential  for  measurement  of  the  engineering 
properties  of  seafloor  sediments. 
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GEOPHYSICS  APPLIED  TO  GEOTECHNICAL  PROBLEMS  IN  A 
MARINE  ENVIRONMENT 

A  CASE  STUDY:  MONTEREY  BAY,  CALIFORNIA 

B.  B.  Barnes,  R.  F.  Corwin,  T.  G.  Hildenbrand,  L.  Jackson 
R.  Kessler,  W.  Takeyama,  M.  Hornick  and  R.  Jenkins 

1.  INTRODUCTION 

The  purpose  of  this  work  is  to  demonstrate  the  value  of  applying 
geophysical  remote  sensing  techniques  to  the  classification  of  engineer¬ 
ing  and  mass  physical  properties  of  seafloor  sediments  commonly  referred 
to  as  "geotechnical  properties"  and  to  suggest  specific  studies  in  which 
the  methods  could  be  further  utilized  to  obtain  geotechnical  information 
necessary  for  seafloor  excavation,  such  as  in  mining  marine  mineral  depo¬ 
sits  of  shallow  and  deep  sea  origin. 

Tools  and  techniques  recently  developed  in  this  laboratory  for  ob¬ 
taining  acoustical  and  electrical  resistivity  measurements  on  target  areas 
in  San  Francisco  Bay  (Barnes  et  al . ,  1972),  Baja  California,  Mexico  (Cor¬ 
win  and  Conti,  1973)  and  more  recently  in  Monterey  Bay,  California,  have 
greatly  enhanced  the  usefulness  of  geophysical  methods  in  determining 
subtle  variations  observed  in  a  wide  variety  of  sediment  types  from  both 
estuarine  and  continental  shelf  environments.  Acoustical  and  electrical 
resistivity  variations  within  individual  sediment  types,  detectable  by 
applying  fine-grained  surveying  practices  and  employing  advanced  instru¬ 
mentation,  can  reflect  information  useful  in  depicting  geotechnical  var¬ 
iations  occurring  in  these  two  environments.  Thus,  the  methods  described 
herein  are  powerful  tools  for  use  in  site  selection  for  undersea  founda¬ 
tion  emplacement  and  in  predicting  a  geologic  situation  before  seafloor 
excavation  and  mining.  To  demonstrate  this,  the  present  report  will  give 
as  an  example  geotechnical  variations  obtained  using  geophysical  methods 
which  are  verified  by  sediment  analyses  of  core  samples  in  selected  test 
sites.  The  analyses  depict  several  different  mass  physical  properties 
and  engineering  properties,  each  systematically  obtained  from  diverse 
kinds  of  sediments  located  within  a  geographical  province  in  the  primary 
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test  area  located  in  Monterey  Bay,  California.  An  additional  advantage 
of  applying  advanced  sampling  equipment  to  enable  more  accurate  and  re¬ 
liable  correlations  of  acoustical  and  electrical  resistivity  measurements 
to  sediment  analyses  is  also  described  in  the  report. 

The  distribution  of  changes  In  mass  physical  properties,  consist¬ 
ing  of  porosity,  density,  mean  grain  size,  and  sorting  coefficients,  and 
of  engineering  properties,  such  as  measured  by  penetrometer  and  vane 
shear  apparatus,  is  made  complex  by  a  host  of  oceanographic  operators  — 
pressure,  currents,  mineralogy,  biology,  geochemistry,  and  structural 
deformation,  to  mention  a  few.  The  complexity  is  further  enhanced  by 
the  interrelation  and  interdependence  of  parameters  related  to  the 
sediment  properties.  In  this  study  we  are  mainly  concerned  with  the 
static  situation  —  the  seafloor  as  it  is  in  one  place  in  space  and  time. 

During  the  last  decade,  interest  in  geotechnical  measurements  has 
been  steadily  growing  because  of  the  utilization  of  our  coastal  zone  and 
adjacent  shelf  for  engineering  purposes,  such  as,  for  example,  emplace¬ 
ment  of  oil  pipelines,  drilling  platforms,  bridge  and  dam  foundations, 
and  submerged  nuclear  power  plants.  More  recently  the  United  States  min¬ 
ing  Industry  has  Indicated  an  interest  in  developing  more  rapid  and  accu¬ 
rate  methods  to  obtain  geotechnical  data  in  deep  sea  areas  where  ferroman¬ 
ganese  deposits  are  located.  Although  the  feasibility  of  applying  remote 
sensing  instrumentation  to  geotechnical  problems  is  still  widely  debated, 
considerable  progress  has  been  made  in  the  development  of  methods  and  in¬ 
strumentation  for  this  purpose.  We  still  have  insufficient  knowledge  of 
the  full  potential  to  reliabily  predict  sediment  property  variations  and 
to  understand  the  mechanisms  that  cause  such  changes  to  occur.  Knowledge 
can  only  be  derived  from  making  many  empirical  measurements  in  the  actual 
environment,  utilizing  computed  analyses  and  the  best  in  sediment  sampl¬ 
ing  equipment  in  conjunction  with  advanced  geophysical  tools  and  methods. 

Elements  in  this  study  consist  of:  (1)  two  separate  acoustic  ex¬ 
periments  based  on  continuously  mapping  reflection  coefficients  by  a  ship 
underway  and  using  in  situ  apparatus  implanted  in  the  seafloor  capable  of 
generating  and  detecting  P  and  S  waves  in  a  saturated  medium,  (2)  an  elec 
trical  resistivity  experiment  utilizing  especially  configured  arrays  for 


obtaining  in  situ  and  surface  towed  direct  current  resistivity  measure¬ 
ments  that  indicate  sediment  conductivity  taken  in  place  and  by  a  shi? 
underway  and,  (3)  the  use  of  unique  sampling  equipment  for  obtaining 
cores  of  the  near  surface  seafloor  sediments.  The  report  defines  in  de¬ 
tail  the  interrelationship  and  correlation  of  data  obtained  from  each  of 
these  three  sources  and  relating  these  data  to  the  overall  objective  and 
purpose  of  the  project. 

1.1  Discussion  of  Problems 

Application  of  geophysics  to  the  solution  of  marine  civil  engineer¬ 
ing  problems,  or  "geotechnics",  constitute  the  subject  of  this  report. 

The  point  of  view  is  that  of  the  engineer,  and  the  earth  sciences,  par¬ 
ticularly  geophysics,  have  been  brought  into  the  engineering  pattern  only 
when  they  have  direct  bearing  upon  a  problem.  Our  goal  is  to  present 
only  those  basic  geotechnical  problems  connected  with  depicting  the  natu¬ 
ral  environment  of  an  engineering  structure  and/or  trafficabil ity  of  sea¬ 
floor  systems,  such  as  mining  implements,  that  may  be  used  for  excavation. 
The  case  history  of  Monterey  Bay  has  been  used  to  elucidate  the  value  of 
applying  geophysics  to  standard  practices  for  determining  seafloor  prop¬ 
erties. 

Several  cases  of  differential  settlement  of  oil  drilling  platforms 
in  the  Gulf  of  Mexico  and  huge  pipelines  shifting  from  their  anchorages 
In  the  North  Sea  and  Mediterranean  Sea  have  resulted  ir.  the  loss  of  pro¬ 
duction  and  a  considerable  amount  of  revenue.  Before  these  events,  not 
all  marine  engineers  could  see  clearly  that  the  design  of  a  structure 
should  be  preceded  by  a  careful  study  of  its  environment,  particularly 
the  foundation  material  on  which  the  structure  was  to  be  placed.  The 
problems  are  more  complicated  than  those  encountered  on  land  because  of 
the  overlying  water  column  and  unusual  behavior  of  marine  sediments.  The 
past  15  years  have  seen  great  emphasis  placed  towards  diagnosing  in  ad¬ 
vance  the  state  of  sediment  properties.  With  the  advent  of  offshore  min¬ 
ing  on  our  continental  shelves  and  in  the  deep  sea,  undoubtedly  more  re¬ 
quirements  for  rapid,  but  accurate  determinations  will  be  necessary;  thus. 


a  new  type  of  technology  is  needed  to  assist  standard  methods  of  geotech¬ 
nical  measurements. 

1.2  Operational  Constraints 

Certain  constraints  to  applying  geophysical  methods  must  be  identi¬ 
fied  before  proceeding  further.  Acoustical  measurements  taken  by  a  ship 
underway  are  particularly  affected  by  topogranhy,  acoustic  beam  width, 
power,  system  stability  in  a  hostile  environment,  and  water  depth,  in 
situ  measurements  are  likewise  affected,  but  to  a  lesser  degree,  their 
greatest  limitation  being  time  and  cost  factors. 

A  level  set  *  r  terrain  is  most  desirable  for  acoustic  measure¬ 
ments.  Micro  topoo.  ,phic  changes  have  greater  influence  on  acoustic 
sources  that  possess  characteristics  of  narrow  beam  width  and  high  fre¬ 
quency.  Ideally,  a  flat  bottom  with  little  variation  of  sediment  type 
produces  the  most  reliable  data.  Water  depth  becomes  an  important  factor 
if  source  and  receiver  transducers  are  on  the  surface.  Attenuation  of 
energy  by  spreading  and  absorption  are  negllble  in  shallower  water  but 
can  become  an  appreciable  factor  with  depth.  The  stability  of  the  sys¬ 
tems,  in  particular  the  source  and  receiver  orientation  to  the  reflect¬ 
ing  surface,  are  obviously  a  consideration  in  obtaining  accurate  signal 
ratios  between  the  transmitted  and  received  (reflected)  energy.  In  the 
case  of  the  towed  electrical  array,  similar  constraints  apply.  In  this 
instance,  the  sensitivity  of  the  measurement  is  affected  causing  the  in¬ 
formation  to  be  less  accurate  than  in  situ  measurement.  These  tolerances 
i,\nd  the  errors  induced  will  be  described  in  this  report. 

2.  PREVIOUS  INVESTIGATIONS 
2.1  Application  of  Acoustics 

The  earliest  work  in  applying  reflection  acoustical  methods  began 
during  World  War  II  when  the  Division  of  War  Research  and  several  Navy 
research  agencies  locked  into  submarine  warfare  applications.  Slant 
angle  measurements  at  various  frequencies  and  distances  provided  the 
first  data  on  the  seafloor  geological  environment.  Attention  was  focused 


on  two  aspects  relating  wave  propagation  phenomena  and  sediments.  In 
the  period  following  the  war  several  investigators  (Liebermann,  1948; 

Ui ick,  1954;  Urick  and  Saling,  1962;  Mackenzie,  1960;  McKinney  and  Ander¬ 
son,  1964;  and  Jones  et  al . ,  1964)  directed  their  efforts  toward  identi¬ 
fying  the  characteristics  of  sound  transmission  and  reflection.  Atten¬ 
tion  was  not  given  to  sediment  properties  until  the  mid- 1950 ’ s 3  when 
investigations  by  Hamilton  (1956),  Shumway  (1960);  Sarmiento  and  Kirby 
(1962),  Richards  (1962),  and  Nafe  and  Drake  (1963)  produced  much  of  the 
deductive  theory  relating  sediment  properties  to  acoustics.  Direct  mea¬ 
surements  made  by  a  ship  underway  did  not  begin  until  the  mid-1960's 
(Breslau,  1964,  1967;  Taylor-Smith  and  Li,  1966;  Taylor-Smith,  1968). 

Recent  work  by  Pawlowicz  (19711  indicates  that  it  is  possible, 
with  caution,  to  describe  the  porosity  of  seafloor  sediments  in  terms  of 
reflectivity,  Furthermore,  the  strong  correlation  which  exists  between 
reflectivity  and  porosity  can  be  used  to  infer  qualitatively,  the  poten¬ 
tial  load-bearing  properties  of  the  substrate.  If  a  more  quantitative 
description  is  desired,  computer  data  processing,  either  in  the  laboratory 
or  aboard  ship,  could  be  used.  Pawlowicz  suggests  that  it  should  be  feas¬ 
ible  to  relate  porosity  to  sediment  type  in  general  terms,  i.e.,  high  por¬ 
osity  relates  to  silts  and  clays,  and  low  porosity  to  sands  and  gravels; 
however,  for  porosities  greater  than  55  percent,  there  is  a  great  overlap 
in  sediment  type  and  precise  designations  may  be  inconclusive  (Taylor- 
Smith  and  Li,  1966;  Faas ,  1969). 

The  Mine  Defense  Laboratory  has  developed  a  Sea  Bottom  Classifier 
(SBC)  which  determines  the  softness  of  a  sea  bottom  by  measuring  the  elong¬ 
ation  of  reflected  acoustic  pulses  using  a  standard  depth  sounder  (Stanley 
and  Harris,  1968).  If  a  short  narrow-beam  acoustic  pulse  is  transmitted 
into  a  hard  bottom,  the  received  echo  is  about  the  same  length  as  the 
transmitted  pulse.  On  the  other  hand,  in  a  soft  bottom,  the  acoustic 
pulse  penetrates  the  bottoms,  and  volume  reverberation  within  the  bottom 
causes  elongation  of  the  echo  pulse  proportional  to  pulse  penetration. 

Li  and  Taylor-Smith  (1969)  believed  that  by  profiling  with  a  re- 
peatmg  broad-band  acoustic  pulse  source  over  areas  where  a  wedge  of 
sediment  was  underlain  by  a  strong  reflector  it  should  be  possible,  by 


comparing  a  suo-bottom  reflection  (and  in  some  cases  its  associated  mul¬ 
tiples)  to  the  pulse  reflected  from  the  sea-floor  and/or  to  its  adjacent 
reflections,  that  one  could  give  a  reasonable  analysis  of  the  bulk  prop¬ 
erties  of  the  sea-floor  medium  in  terms  of  variations  In  the  frequency 
spectrum.  The  values  so  obtained  from  Icboratory  and  in  eitu  probe  mea¬ 
surements,  would  thus  provide  data  about  the  mechanical  properties  of  the  » 

sediment. 

Sub-bottom  profiling  devices  as  conventionally  used,  r.o  matter  how 
complex,  can  only  provide  an  assessment  of  the  geology  of  the  sediments,  -  v 

not  the  sediment  types  or  their  properties.  While  it  should  be  possible 
to  obtain  a  feeling  for  the  hardness  of  the  sea-floor  from  observed  multi¬ 
ples  on  a  sub-bottom  record  and,  to  a  certain  extent,  the  material  sizes 
making  up  the  surface  shown  by  scattering  or  lack  of  penetration  observed 
on  a  record.  In  analysing  multiples,  few  attempts  have  related  the  energy 
decay  rate  of  multiples  to  varying  sediment  properties. 

Other  studies  by  Li  and  Taylor-Smith  (1966)  and  Breslau  (1965),  and 
Shumway  (1960),  have  also  pointed  out  that  frequency  analysis  of  seismic 
energy  reflected  from  the  bottom  may  provide  the  information  needed  to 
absolutely  fix  acoustic  reflection  values  to  sediment  types.  Experimental 
studies  by  Taylor-Smith  and  LI  (1966  and  1969),  and  Breslau  (1965)  have 
used  spectrum  analysis  measurements  utilizing  the  attenuation  factors  of 
various  sediment  types  versus  grain  sizes  with  the  comparison  of  the  vari¬ 
ous  frequency  components  of  an  acoustical  source. 

Results  of  work  done  by  McCann  (1966)  and  Buchan  et  al .  (1967)  show 
that  attenuation  plotted  against  grain  size  varies  for  various  frequencies 
within  the  low-frequency  sound  source  range  (fig.  1)  According  to  McCann 
(1966),  two  mechanisms  are  responsible  for  the  acoustic  energy  loss  in 
saturated  sediments:  (1)  A  solid  friction  loss  for  sediments  with  grain-  % 

to-grain  contact,  this  being  characterized  by  a  linear  variation  of  the 
attenuation  coefficient  with  frequency,  and  (2)  a  viscous  loss  caused  by 
Interaction  between  solids  and  liquids,  this  being  characterized  by  an 
attenuation  coefficient  which  is  : > oportional  to  the  square  root  of  the 
frequency. 
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Sediment  mean  diameter  (mm) 


Figure  1.  Sound  attenuation  plotted  against  mean  grain  diameter. 

(McCann,  1966). 

Reflection  coefficients  measured  continuously  by  a  ship  underway 
have  been  made  Breslau  (1965),  Li  and  Taylor-Smith  (1965),  and  Barnes  et 
al .  (1972).  These  consisted  of  measuring  the  amplitudes  and/or  energy  in 
the  transmitted  and  reflected  pulses  and  correlating  attenuation  to  vari¬ 
ous  bottom  sediment  properties. 

Another  way  of  ascertaining  sediment  properties  is  by  in  situ  mea¬ 
surement,  from  which  one  can  obtain  direct  information  on  the  compres- 
sional  wave  velocity,  shear  wave  velocity,  attenuation,  and  Poisson's 
V-  ratio.  Direct  measurement  of  the  shear  wave  velocity  in  ocean  sediments 

has  been  restricted  to  the  seismology  of  high  energy  sources  that  produce 
w  returns  from  the  deeply  buried  material  and  no  returns  from  the  soft  sup- 

m 

erficial  sediment.  This  is  because  shear  waves  do  not  propagate  at  all 
in  a  fluid  and,  consequently,  are  transmitted  poorly  in  a  saturated  por¬ 
ous  medium.  The  shear  wave  velocity  for  the  saturated  sediments  of  the 
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sea  floor  has  been  determined  indirectly  through  measuring  in  situ  the 
Stoneley  wave  velocity  (Hamilton  et  al . ,  1970;  Davies,  1965;  Jones, 

1958;  Cohen,  1968;  and  Bieda,  1970).  It  has  been  demonstrated  in  labor¬ 
atory  experiments  that  the  generation  of  torsional  waves  by  a  viscoelas- 
tometer  provides  Information  on  the  shear  wave  velocity  in  saturated  sed¬ 
iments. 

The  compressional  wave,  unlike  the  transverse  (shear)  wave,  is 
transmitted  readily  in  water  and  ir,  saturated  sediments.  Its  velocity 
in  marine  sediments  has  been  documented  in  several  reports  (Nafe  and 
Drake,  1963;  Hamilton,  1970a;  Lewis,  1971;  Horn  et  al.,  1968).  If,  in 
addition  to  the  compressional  wave  velocity  and  the  shear  wave  velocity, 
density  is  known,  one  can  calculate  the  dynamic  values  of  rigidity,  bulk 
modulus,  Young's  modulus,  and  other  elastic  constants  (Hamilton,  1970a; 

Nafe  and  Drake,  1963).  Due  to  the  difficulty  in  measuring  the  shear 
wave  velocity,  the  elastic  constants  have  been  derived  by  assuming  a 
value  of  the  sediment  bulk  modulus  from  an  aggregate  theory  (Hamilton, 

1971;  Laughton,  1957;  Lewis,  1971). 

Other  marine  sediment  properties  of  Interest  are  porosity,  grain 
size,  shear  strength,  and  percent  of  clay,  sand,  or  silt.  Acoustical 
properties  are  determined  by  these  physical  properties  and,  thus,  have 
been  empirically  related  to  them.  The  reader  is  referred  to  excellent 
discussions  on  this  subject  in  Hamilton  (1970b,  1971),  Horn  et  al .  (1968), 
Schreiber  (1968),  and  Lewis  (1971). 

2.2  Application  of  Electrical  Resistivity 

Apart  from  acoustics  there  has  been  an  appreciable  amount  of  effort 
devoted  to  resistivity  techniques  applied  to  prediction  of  porosity  in 
marine  sediments. 

The  earliest  marine  methods  evolved  from  the  application  of  theory 
and  technology  used  in  land  borehole  logging  to  marine  investigations, 
(Bouma  et  al.,  1971;  Erchul  and  Nacci,  1971;  Kermabor  et  al.,  1969;  Wyllie, 
1963).  A  horizontal  array  for  measuring  seafloor  resistivity  has  resulted 
from  this  evolution  (Corwin  and  Conti,  1973).  For  marine  use,  a  horizontal 
array  may  be  towed  over  the  sea  surface  or  deployed  on  the  seafloor.  The 


8 


. . 


- . .  riirniirttMiB-i 


jfc  ^  i  ^ 


first  application  provides  a  fast  and  convenient  method  of  mapping  hori¬ 
zontal  changes  of  average  bottom  resistivity  over  large  areas.  An  array 
deployed  in  the  second  manner,  however,  provides  a  more  detailed  picture 
of  the  bottom  resistivity  distribution  at  a  given  location. 

There  has  been  little  previous  use  of  horizontal  resistivity  arrays 
on  the  seafloor  An  array  of  fixed  spacing  dragged  over  the  seafloor 
off  the  coast  of  Cornwall,  England,  was  used  to  map  lateral  changes  in 
seafloor  resistivity  (Marke,  1965).  The  fixed  electrode  spacing  used  by 
Marke,  however,  precluded  the  acquisition  of  information  about  the  verti¬ 
cal  distribution  of  bottom  resistivity.  Schlumberger  and  Leonardon 
(1934)  used  a  bottom-deployed  array  to  measure  depth  to  bedrock  below 
seafloor  sediments.  Dunlap  and  Johnson  (1958)  mapped  offshore  faults  and 
salt  domes,  also  using  a  fixed  array  dragged  over  the  seafloor. 

The  deployment  of  a  surface  towed  resistivity  array  in  work  done  by 
this  laboratory  is  believed  to  be  the  first  such  application  using  this 
mode.  Several  other  electrical  studies,  including  self  potential,  direct 
current  resistivity  and  electromagnetic  sounding  have  been  carried  out  in 
this  laboratory  and  are  reported  in  Barnes  et  al.  (1972),  Beyer  (1972), 
and  Crowin  and  Conti  (1973). 


3.  THEORY  AND  TECHNIQUE 
3.1  Reflection  Coefficient  Mapping 


A  simplified  approach  to  a  very  complex  physical  process  involving 
wave  propagation,  attenuation,  and  absorption  has  been  used  in  develop¬ 
ing  the  model  of  the  sea  floor-water  interface.  The  sea  floor  is  con¬ 
sidered  a  plane  interface  between  two  fluids.  Energy  is  focused  on  a 
zone  of  about  one  square  meter  surface  area.  The  backscattering  strength 
associated  with  this  surface  area  depends  greatly  on  the  roughness  of  the 
sea  bottom.  A  smooth  bottom  acts  as  a  mirror  such  that  the  reflected 
energy  is  concentrated  in  a  specular  direction.  A  rough  bottom  will  re¬ 
flect  the  energy  more  evenly  in  all  directions.  It  is  important  to  note 
that  the  terms  "smooth"  and  "rough"  must  be  considered  in  relation  to 
the  acoustical  wave  length.  Where  the  wave  length  is  only  a  few  centi¬ 
meters,  it  is  seen  that  a  100  percent  smooth  bottom  is  rere.  Even  small 
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stones,  holes  or  ripples  will  make  the  bottom  rough  in  the  sense  consid¬ 
ered  here.  When  a  sound  wave  strikes  the  sea  bottom,  part  of  the  inci¬ 
dent  energy  will  be  reflected.  The  reflected  energy  is  spread  in  all 
directions  with  a  typical  distribution,  as  shown  in  figure  2.  The  dis¬ 
tribution  is  totally  dependent  on  the  roughness  of  the  sea  bottom.  Fig¬ 
ure  3  illustrates  how  the  roughness  of  the  topography  can  influence  the 
bottom  backscattering  strength  as  a  function  of  the  grazing  angle.  The 
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Figure  2.  Seafloor  scattering  —  specular  energy. 
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Figure  3.  Typical  curves  for  bottom  backscattering  strength. 

(After  Bodholt,  1969). 


angle  of  incidence  is  seen  to  also  play  an  important  role  in  determining 
the  strength  of  the  echo  signal.  The  experiments  described  in  this  re¬ 
port  were  all  conducted  with  the  source  and  receiver  set  up  to  record 
normal  incidence.  It  can  be  shown  that  theoretically  the  most  accurate 
measurements  of  bottom  backscattering  strength  should  be  obtained  with  a 
narrow  beam  transducer.  If  the  beam  width  is  5-10°,  the  bottom  backscat¬ 
tering  strength  will  vary  only  a  few  dB  or  less  within  the  beam,  of  course 


depending  on  the  type  of  bottom.  The  magnitude  of  the  bottom  sampling 
area  will  depend  on  the  pulse  length  of  the  transmitted  sound.  The 
shorter  the  pulse  the  smaller  the  sampling  area.  Figure  4  illustrates 
the  beam  width,  pulse  length,  and  area  concept. 

In  the  exercise  conducted  for  this  investigation,  two  sources  are 
employed.  One  is  an  uncontrolled  beam  width  emitting  a  .5  millisecond  v 

pulse  length  and  the  other  is  a  wide  beam  (35°)  echo  sounder  having  a 
pulse  length  of  .?i>  millisecond. 
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BW  =  Beam  Width 
PL  =  Pulie  Length 
d A  =  Incremental  Area 
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Figure  4.  Transmitted  sound  guise 
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Details  of  the  mathematical  description  for  acoustic  wave  propaga¬ 
tion,  reflection  and  absorption  (bottom  loss)  is  contained  in  a  previous 
/eport  by  Barnes  at  al .  (1972).  Several  references  were  consulted  and 
the  best  mathematical  expressions  relating  to  the  experiments  conducted 
in  this  investigation  were  taken  from  derivations  by  Breslau  (1567). 

From  the  derivations  the  most  useful  equations  applicable  to  data 
reduction  and  processing  are  represented  by  defining  the  relationship  be¬ 
tween  the  fractional  loss  of  intensity  at  the  sea  floor,  K,  and  the  Ray¬ 
leigh  Reflection  Coefficient,  R,  i.e., 

K  =  R2  (1) 

where  K  is  defined  on  a  peak  pressure  basis, 


K  = 


1 

-a2D 

e 


(2) 


where 

PR  =  pressure  of  the  echo  pulse, 

P$  =  pressure  of  the  source  pulse, 

D  =  transmission  loss  due  to  spherical  spreading,  and 

g-a2D  _  transmission  loss  due  to  dissipative  attenuation  of 
sound  in  seawater. 

An  expression  of  the  fundamental  equation  for  bottom  loss,  BL,  is  defined 
by  Breslau  (1967)  as 


BL  =  SpL  -  EpL  -  TLS  -  TLa  ,  (3) 

where 

BL  =  bottom  loss  (in  dB)  and  is  equal  to  -10  log  K, 

SpL  =  pressure  level  of  the  source  pulse  (in  dB//l  dyne/cm2) 
and  is  equal  to  20  log  P<., 

Ep.  =  pressure  level  of  the  echo  pulse  (in  dB// 1  dyne/cm2)  and 
is  equal  to  20  log  PR, 

TLS  =  spreading  loss  (in  dB)  and  is  equal  to  20  log  2D, 

TL^  =  attenuation  loss  (in  dB)  and  is  equal  to  10  2D  log  e. 
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A  final  expression  relates  bottom  loss,  BL,  to  reflectivity,  R,  expressed 
in  decibels; 


BL  =  20  log  R.  (4) 

Other  manifestations  of  these  basic  expressions  that  include  measurements 
dealing  with  the  energy  versus  time  waveform  of  the  echo  are  described  in 
Breslau  (1967).  Since  this  investigation  was  concerned  only  with  peak 
pressure  measurements,  the  derivations  relating  to  energy  in  the  echo 
strength  are  not  included  in  this  discussion. 

The  computational  basis  for  correcting  data  related  to  the  field 
configurations  (geometry)  of  acoustic  sources  and  their  respective  re¬ 
ceivers,  as  used  in  the  Monterey  Bay  experiment,  are  contained  in  Appen¬ 
dix  A-l.  The  mathematical  expressions  and  constants  derived  were  used 
in  computer  analyses  and  processing  of  the  data  (see  section  7.1.1). 

3.1.1  Relation  of  Bottom  Loss  to  Mass  Physical  Properties 

The  mass  characteristics  of  the  bottom  sediment  are  related  to  bottom 
loss  through  the  Rayleigh  Reflection  coefficient  and  the  acoustic  imped¬ 
ance  contrast  at  the  interface,  and  it  has  been  shown  that  the  Rayleigh 
Reflection  Coefficient  is  related  to  bottom  loss  through  equation  (4) 
(Breslau,  1967).  It  is  possible  to  establish  a  relationship  between  por¬ 
osity  and  bottom  loss  using  the  impedance  which  is  dependent  upon  density 
and  compressional  velocity,  both  factors  that  relate  to  the  porosity  of 
natural  sediments  (Barnes  et  a!.,  1972;  Breslau,  1964,  1967). 

In  summary,  several  investigations  have  ascertained  that  a  linear 
relationship  exists  between  porosity  and  density,  in  particular  that  bot¬ 
tom  loss,  as  deduced  from  many  measurements  made  on  sediments,  bears  a 
definite  relationship  to  the  porosity  of  the  sediment.  And  furthermore, 
since  the  porosities  of  natural  sediments  are  somewhat  related  to  sedi¬ 
ment  grain  size,  it  is  also  possible  to  establish  a  general  relationship 
between  bottom  loss  and  the  geological  properties  of  sediments. 

For  further  discussion  of  the  sediment  bulk  properties  and  varia¬ 
tions  within  the  physical  makeup  of  the  grains,  interstitial  properties, 
etc.,  and  more  specifics  on  how  these  factors  influence  porosity,  the 
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reader  is  directed  to  the  papers  by  Nafe  and  Drake  (1963),  Birch  et  al. 
(1942),  Hamilton  et  al.  (1956),  Shumway  (1960),  Sutton  et  al .  (1967), 
Breslau  (1964),  and  Barnes  et  al .  (1972).  The  information  available  on 
this  subject  is  too  voluminous  to  be  discussed  here. 

3.2  In  Situ  Shear  Wave  Generation 

In  an  earlier  report  (Barnes  et  al.,  1972),  it  was  concluded  that 
the  best  method  to  generate  and  directly  detect  shear  waves  in  saturated 
marine  sediments  is  by  implanting  both  receivers  and  a  SH~source  into  the 
sediments  and  by  timing  the  direct  c^-wave.  Assuring  a  large  response  to 
the  direct  S^-wave,  a  linear  array  of  horizontal  geophones  was  implanted 
with  thn  sensitive  axis  of  each  geophone  parallel  to  the  ground  motion 
of  direct  S^-wave.  The  S^-source  was  placed  at  one  end  of  the  geophone 
array.  Theoretically,  in  a  homogeneous  half  space  the  first  arrival  de¬ 
tected  by  each  horizontal  geophone  will  be  the  direct  SH-wave.  The  shear 
wave  velocity  is,  thus,  obtained  from  a  time  distance  graph  constructed 
from  the  first  arrival  at  each  horizontal  geophone. 

The  shear  wave  source  was  also  employed  to  generate  compressional 
waves.  Little  concern  was  given  to  the  detection  of  the  compressional 
wave  since  it  is  transmitted  and  detected  with  little  difficulty  in  sat¬ 
urated  sediments.  However,  it  was  thought  best  to  implant  source  and  re¬ 
ceivers  into  the  sediment.  This  would  avoid  the  confusion  that  might  ex¬ 
ist  in  differentiations  between  the  direct  compressional  wave  propagating 
through  the  sediment  and  the  direct  compressional  wave  propagating  through 
the  water. 

A  linear  array  of  vertical  geophones  was  used  to  record  the  arrival 
time  of  the  direct  compressional  wave.  Theoretically,  a  vertical  geophone 
will  not  respond  to  the  direct  compressional  wave  if  source  and  receiver 
are  at  the  same  depth;  however,  because  of  the  predictable  increase  in 
density  of  natural  sediments  with  depth,  the  direct  P-wave  will  propagate 
in  an  arcuate  path  and,  thus,  be  detected  by  a  vertically  sensitive  geo¬ 
phone.  The  advantage  of  using  a  geophone  sensitive  to  vertical  dis¬ 
placement  instead  of  a  standard  omnidirectional  geophone  or  a  horizontal 
geophone  is  that  information  on  the  direct  Sy-wave  might  be  obtained 
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from  the  vertical  geophone.  Barnes  et  al.  (1972)  describes  in  detail 
successful  applications  of  this  technique  in  moist  sand  and  saturated 
beach  sands. 

3.2.1  Viscoelastic  Model 

There  remains  the  question  of  what  model  to  employ  for  the  sea  floor. 
Recent  reports  by  Hamilton  (1970a,  1972)  have  clarified  this  problem.  From 
his  data  and  from  an  extensive  literature  search,  Hamilton  (1972)  arrived 
at  the  following  conclusions  for  a  water  saturated  porous  media:  (1)  for 
the  small  stresses  produced  by  sound  waves,  marine  sediments  act  essenti¬ 
ally  as  a  closed  system  (interstitial  water  and  solids  move  together); 

(2)  viscosity  of  the  interstitial  water  is  neglible;  (3)  dispersion  of 
velocity  with  frequency  is  insignificant  in  the  frequency  range  of  inter¬ 
est  for  geophysical  studies;  (4)  Hookean  model  and  the  equations  of  elas¬ 
ticity  that  follow  may  be  used  in  most  cases  to  compute  the  compressional 
and  shear  wave  velocity;  and  (5)  a  linear  viscoelastic  model  with  attenu¬ 
ation  proportional  to  frequency  and  energy  damping  independent  of  fre- 
cuency  can  be  used  if  wave  energy  losses  are  to  be  considered.  These 
conclusions  will  not  be  discussed  in  great  detail  except  the  usefulness 
of  the  Hookean  model  to  determine  wave  velocities. 

In  geophysical  and  soil  mechanic  studies  of  rocks  and  minerals,  the 
equations  of  Hookean  elasticity  have  been  shown  to  adequately  define  the 
compressional  and  shear  velocity  (Anderson  and  Liebermann,  1968)  (Hardin 
and  Richards,  1963).  The  widespread  use  of  the  elastic  model  is  primar¬ 
ily  due  to  the  fact  that  this  model  does  not  provide  for  wave  energy  losses 
and  that  these  losses  are  small  in  most  materials.  To  investigate  whether 
saturated  marine  sediments  can  be  considered  as  a  low-loss  material,  an  ap¬ 
propriate  model  providing  for  energy  absorption  was  studied  and  then  com¬ 
pared  to  the  elastic  model. 

An  appropriate  anelastic  model  is  to  consider  the  seafloor  as  a 
homogeneous,  isotropic  linear  viscoelastic  medium.  Linear  viscoelasticity 
has  been  discussed  in  various  forms  by  Borcherdt  (1972),  Ferry  (1961), 

White  (1965),  and  others.  The  basic  equations  of  linear  viscoelasticity 
were  derived  (Appendix  A-2)  following  the  approach  by  Borcherdt  (1972). 
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A  linear  viscoelastic  model  accounts  for  wave  energy  losses  by  re¬ 
placing  the  Lame  constants  of  the  Hookean  equations  with  complex  Lame  con¬ 
stants.  In  the  expressions  for  rigidity,  yR  +  iy,  and  bulk  modulus, 

KR  +  iKj,  yR  and  KR  denote  the  elastic  response;  and  y^  and  Kj  denote  the 
damping  of  wave  energy,  yj,  yR>  Kj  and  KR  are  all  real  and  positive. 

Anelasticity  can  be  measured  in  terms  of  the  dimensionless  quantity 
Q"1  (specific  attenuation  factor).  It  is  defined  as  the  (27r)-1  times  the 
ratio  of  the  energy  dissipated  per  cycle  of  forced  oscillation  to  the  peak 
energy  stored  during  the  cycle  (Borcherdt,  1972,  p.  12).  The  specific 
attenuation  factor  in  terms  of  velocity  (V),  attenuation  (a),  and  circu¬ 
lar  frequency  (w)  is  found  to  be 


q-i  =  — 2Va/to -  _  (5) 

1  -  (aV/tt)2 

The  specific  attenuation  factor  for  dilatation,  Q”1,  and  shear  Q-1,  are 

p  s 

found  by  substituting  the  appropriate  subscript  p  or  s  in  equation  (5). 
Note  that  because  velocity  is  independent  of  frequency  and  attenuation 
is  proportional  to  frequency,  Q"1  is  independent  of  frequency. 

For  a  low-loss  material  (0  <  Q"1  <,  1),  a V/w  «  1;  and  therefore, 
(5)  may  be  approximated  by 


Q-1 


2Va 

(0 


(6) 


The  expressions  for  the  phase  velocity  of  a  homogeneous  plane  com- 
pressional  wave  (V  )  and  shear  wave  (V$)  ir  a  linear  viscoelastic  medium 
is  given  by: 


(7) 


♦ 


(8) 
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where  $  =  sediment  density 


2(1  +  Q"2) 

q  =  - 1 — 

p  1  +  Q‘2  +1 


(9) 


2(1  +  Q-2) 
1  +  Qs2  +1 


(10) 


As  energy  losses  become  small,  Qp1  (or  Q^1)  -*■  0  and  qp  (or  q$)  approach 
unity.  Thus,  (7)  and  (8)  can  be  simplified  to  the  familiar  expressions 
for  velocities  in  an  elastic  medium: 


Kn  + 


V 


4 

T 


(11) 


(12) 


The  question  remains  whether  ocean  bottom  sediment  can  be  classi¬ 
fied  as  a  low-loss  medium.  If  so,  the  low-loss  approximations,  equations 
(6),  (11),  and  (12)  can  be  employed  to  describe  wave  propagation.  Fig¬ 
ure  5  illustrates  the  low-loss  error  percentages  as  a  function  of  Q"1. 
Values  of  Q"1  ranging  from  approximately  .5  to  2.0  have  been  reported  by 
Hamilton  (1970a,  p.  4048)  for  ocean  bottom  sediments  off  San  Diego,  Cali¬ 
fornia  [assumption  of  no  absorption  in  bulk  (Kj  =  0)  was  used].  Referring 
to  figure  5  and  to  the  viscoelastic  equations  discussed  above,  error  per¬ 
centages  ranging  from  8  to  43  percent  can  be  introduced  in  determining  the 
velocity  of  a  homogeneous  shear  wave,  and  the  error  percentage  for  the  cor¬ 
responding  absorption  coefficient  is  between  15  and  184  percent.  Because 
these  large  errors  are  not  permissible,  the  linear  viscoelastic  equations 
are  needed  in  determining  the  shear  wave  velocity  or  Lame  constants. 

Hamilton  (1972,  p.  635)  also  reports  Q"1  values  ranging  from  0.002- 
0.044  for  sediments  off  San  Diego.  The  low-loss  approximations,  equations 
(6)  and  (11)  can  be  used  without  any  significant  error. 
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Figure  5.  Low-loss  error  percentages  for  the  velocity  cmd  absorption  co¬ 
efficient  of  a  homogeneous  shear  wave  with  0.05  <  •  Q~l  <  1.0. 

Density,  Vp,  VgJ  ap,  and  ag  were  measured  to  calculate  Q"1 ,  yR,  y 
Kj,  KR,  Poisson's  Ratio,  and  Young's  modulus.  Appendix  A-3  lists  equa¬ 
tions  that  can  be  employed  to  determine  various  sediment  properties  from 
the  measured  properties. 


3.3  Direct  Current  Resistivity 

The  electrical  resistivity  of  a  geologic  material  is  related  to  the 
porosity  of  the  material  (Archie,  1942)  with  low  resistivity  implying  high 
porosity.  Resistivity  (p)  is  numerically  equal  to  the  electrical  resist¬ 
ance  between  opposite  faces  of  a  cube  of  the  material  of  unit  dimensions, 
and  usually  is  measured  in  ohm-centimeters  (ohm-cm)  or  ohm-meters  (ohm-m) 
(Grant  and  West,  1965).  As  porosity  is  related  to  sediment  composition 


and  engineering  properties,  in  situ,  measurement  of  the  vertical  variation 
of  the  electrical  resistivity  of  marine  sediments  can  provide  Information 
about  sea-floor  structure. 

Resistivity  techniques  are  well  established  on  land,  and  increasing 
use  Is  being  made  of  them  In  the  offshore  environment.  Measurements  of 
resistivity  and  porosity  of  sea-floor  sediments  Indicate  thit,  except  for 
some  highly  conductive  clays,  the  formation  factor  of  a  mrrine  sediment 
predicts  the  porosity  of  the  sediment  to  within  5  or  10  percent.  Forma¬ 
tion  factor  F  Is  defintl  by  the  equation: 


F  = 


» 


(13) 


where  p$  Is  the  resistivity  of  the  sediment  and  Is  the  resistivity  of 
the  overlying  water.  Figure  6  summarizes  some  of  these  measurements.  Mea¬ 
surement  of  the  vertical  resistivity  distribution  of  the  sea  floor,  there¬ 
fore,  provides  a  method  of  determining  sub-bottom  porosity,  and  thus  geo¬ 
logic  and  engineering  properties. 


Figure  6.  Porosity  vs  formation 
factor  for  marine  sediments. 
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Resistivity  is  measured  on  land  by  borehole  logging,  or  by  the  use 
of  an  expanding  horizontal  surface  array.  In  boreholes,  a  small  probe 
lowered  into  the  hole  provides  a  continuous  record  of  the  resistivity  of 
the  material  surrounding  the  probe  (Wyllie,  1963).  This  method  has  been 
adapted  successfully  to  marine  use  by  forcing  a  thin  probe  into  the  sea¬ 
floor  sediments  (Bouma  et  al.,  1971;  Erchul  and  Nacci ,  1971;  Kermabon  et 
al.,  1969),  as  shown  schematically  in  figure  7.  The  use  of  a  logging 
probe  in  marine  sediments  gives  a  relatively  precise  determination  of  re¬ 
sistivity  as  a  function  of  depth,  but  does  offer  some  drawbacks.  The  depth 
of  measurement  is  limited  to  the  depth  to  which  the  probe  can  be  inserted 
into  the  sea  floor,  which  may  be  small  where  sediment  porosity  is  low.  The 
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Current  I  I*  fed  Into  the  sediment*  through  current  electrode*  A,  and  cau*e* 
a  voltage  AV  to  appear  between  potentlol  electrode*  B  Sediment  resistivity 
o  le  then  determined  from  the  formulo 


where  K  It  o  foctar  depending  on  the  geometry  of  tne  el*. ‘rod#  arrongement. 


Figure  7.  Well  logging  probe  used  to  measure  resistivity  of 

marine  sediments . 


volume  of  bottom  material  sampled  is  relatively  small  for  a  single  inser¬ 
tion  of  the  probe,  so  interpretation  errors  may  result  if  the  sampled 
point  is  not  representative  of  the  surrounding  area. 

Surface-deployed  horizontal  arrays  also  are  used  to  measure  verti¬ 
cal  resistivity  profiles  on  land  (Grant  and  West,  1965;  Van  Nostrand  and 
Cook,  1966).  As  shown  in  figure  8,  this  method  employes  two  current  elec¬ 
trodes  to  feed  low-frequency,  commutated  direct  current  into  the  ground, 
and  two  (or  more)  potential  electrodes  to  measure  the  potential  difference 
generated  at  the  surface  of  the  earth  by  the  impressed  current  flow.  The 
potential  difference  is  a  function  of  the  electrode  configuration  and 
spacing,  and  of  the  subsurface  resistivity  distribution. 

For  the  equispaced  Wenner  array,  shown  in  figure  8,  the  "apparent 
resistivity",  pfl,  of  the  earth  is  given  by  the  formula  (Van  Nostrand  and 
Cook,  1966) 

pa  =  2?ra  T“  *  (I4) 

where 

pa  =  apparent  resistivity, 
a  =  electrode  spacing, 
aV  =  measured  potential  difference,  and 
I  =  impressed  current. 

By  expanding  the  array,  the  current  is  forced  to  flow  deeply  into  the 
earth,  and  the  apparent  resistivity  at  larger  electrode  spacings  is  indi¬ 
cative  of  the  material  at  greater  depth.  The  field  data  is  interpreted 
by  plotting  measured  apparent  resistivity  Pjj  against  electrode  separation 
a,  and  matching  these  curves  with  those  contained  in  standard  catalogs  of 
published  curves  (Mooney  and  Wetzel,  1956),  as  shown  in  figure  9. 

The  use  of  an  expanding  horizontal  array  for  sea-floor  resistivity 
measurements  (fig.  13)  offers  a  number  of  advantages  over  the  use  of  a 
probe.  The  depth  below  the  sea  floor  to  which  information  may  be  obtained 
is  limited  only  by  the  size  of  the  array  and  the  amount  of  current  avail¬ 
able,  and  the  measuring  process  does  net  disturb  the  sediments.  Also,  it 
is  a  simpler  operation  to  deploy  an  electrode  array  on  the  bottom  than  to 


t 


T 


V|  ond  V,  art  pofentlol  electrodes 
Xi  ond  I2  art  current  electrodes 


B)  Large  Separofian  a. 


Current  t  le  fed  Into  the  ground  through  the  current  electrodes  I|  ond  I2l  creating 
a  potential  field  on  the  turfoce  of  the  ground.  The  voltage  difference,  AV,  between 
the  potenflol  electrodes  V|  and  V2  Is  0  function  af  the  resistivity  of  the  ground 
below  the  electrodes.  The  apparent  resistivity  of  the  ground,  Pa,  Is  close  to  P, 
when  the  separation  a  Is  smoll ,  and  approaches  Pz  os  a  becomes  very  lorge 


Figure  8.  tenner  resistivity  array 
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Figure  9.  Standard  Warmer  resistivity  curves 
(after  Mooney  and  Wetzel,  1956). 

drive  a  p-obe  into  the  sea  floor,  and  a  larger  volume  of  material  is  sam¬ 
pled  at  each  data  point.  On  the  other  hand,  a  probe  gives  a  more  precise 
determination  of  resistivity  and  structure  at  a  given  point.  Also,  data 
interpretation  is  easier  and  more  accurate  with  a  probe  because  there  is 
no  need  to  match  curves  to  obtain  the  resistivity-depth  relationship  from 
the  field  data.  The  probe  and  horizontal  array  methods,  then,  are  compli¬ 
mentary,  with  the  probe  offering  greater  precision,  and  the  horizontal 
array  greater  penetration  and  convenience  and  less  sediment  disturbance. 

For  marine  use,  a  horizontal  array  may  be  towed  over  the  sea  sur¬ 
face  cr  deployed  on  the  sea  floor  (fig.  11).  A  surface-towed  array  pro¬ 
vides  a  fast  and  convenient  method  of  mapping  horizontal  changes  of  average 
bottom  resistivity  over  large  areas.  An  array  deployed  on  the  sea  floor, 
however,  provides  a  more  detailed  picture  of  the  bottom  resistivity  dis¬ 
tribution  at  a  given  location.  Both  surface-towed  and  bottom-deployed 
arrays  were  tested  for  this  study. 
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The  surface  array  Is  easier  ta  taw,  but  the  sea  flaar 
array  is  closer  to  the  bottom,  so  it  can  "see"  bottom 
resistivities  at  smaller  values  at  the  separation  a 


Figure  10.  Surface  and  seafloor  resistivity  arrays. 


4.  ELECTRONIC  SYSTEMS  AND  EQUIPMENT 
4.1  Reflection  Coefficient  Mapping  Experiment 

4.1.1  Acoustic  Sources 

The  types  of  sources  to  be  used  for  predicting  the  geology  of  the 
seafloor  and  its  substructure  vary  depending  on  the  type  of  sediments  and 
and  the  problem  being  examined.  For  reflectivity  measurements  it  is  de¬ 
sirable  to  cover  a  broad  spectrum  containing  varying  energy  levels,  band 
width,  pulse  width,  and  beam  concentration. 
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In  standard  sub-bottom  profiling  techniques,  the  aim  is  to  get  the 
highest  energy  level  possible  in  a  pulse  of  fairly  short  duration.  While 
these  requirements  are  not  entirely  divorced  from  frequency  considerations, 
frequency  enters  into  the  argument  mainly  in  terms  of  both  discrimination 
and  penetration.  It  is  a  common  belief  that  for  discrimination  a  high 
frequency  is  needed,  while  for  penetration  the  converse  is  true;  in  fact 
some  sub-bottom  profilers  now  commercially  available  are  designed  to  op¬ 
erate  with  a  fixed  frequency  pattern  to  achieve  a  desired  combination  of 
penetration  and  discrimination.  In  general  this  is  a  misapplied  concep¬ 
tion,  particularly  if  one  desires  to  consider  application  to  solving  geo¬ 
technical  problems.  A  seismic  pulse  should  have  a  short  time  duration 
for  resolution  and  a  large  energy  content  to  combat  noise.  Multi-channel 
stacking  is  one  way  of  achieving  an  effective  increase  in  signal  energy 
relative  to  uncorrelated  noise  for  a  given  pulse  width.  Hutchins  (1969), 
accomplished  the  same  effect  by  using  a  large  time-bandwidth  (TW)  product 
signal.  The  idea  behind  this  appraoch  is  to  generate  a  long  duration  sig¬ 
nal  having  a  bandwidth  appropriate  for  the  required  resolution,  to  obtain 
more  signs!  energy,  and  then  to  effectively  time  compress  the  reflected 
signals  by  matched  filter  detection.  Such  a  signal  can  be  band  limited 
random  noise,  frequency  modulated  chirp  or  any  other  sufficiently  compli¬ 
cated  sequence.  This  concept  has  been  applied  in  transmitting  large 
TW  signals  used  in  radar  and  sonar  to  get  maximum  information  from  a 
single  transmitted  pulse.  This  is  clearly  necessary  when  targets  are 
moving  and  probably  taking  evasive  action. 

When  using  seismic  methods  to  map  subsurface  geologic  structures 
buried  beneath  the  ocean  floor,  it  is  often  desirable  to  be  able  to  con¬ 
trol  the  energy  versus  frequency  spectrum  of  the  sonic  pulse  (Caulfield, 
1962).  Sediments  have  attenuation  factors  that  are  frequency  dependent 
(Shumway,  1960).  As  greater  penetration  is  required,  lower  frequencies 
must  be  used.  Since  the  ability  of  a  sonic  system  to  resolve  range  dif¬ 
ferences  depends  on  the  wave  lengths  being  employed,  it  follows  there 
will  be  some  pulse  spectrum  that  provides  the  best  compromise  between 
penetration  and  resolution. 


In  this  investigation,  a  discrete  frequency  wide-beam  trans¬ 
ducer  operating  at  12  kHz  was  selected  to  obtain  a  correlation  with 
the  values  obtained  by  Breslau  (1964),  and  a  broadband  multiple-elec¬ 
trode  sparker  source  was  tested  as  an  alternate  approach  using  a  low 
frequency,  omnidirectional  beam. 

Pawlowicz  (1971)  noted  the  principal  differences  between  the  con¬ 
tinuous  reflection  profiling  systems  and  those. systems  used  for  echo¬ 
sounding  of  water  depth  are  the  requirements  of  high  acoustic  power  out¬ 
put  and,  generally,  lower  frequency.  Although  penetrations  of  90  ft  into 
the  seafloor  have  been  obtained  with  frequencies  as  high  as  14  kHz  (Murray, 
1967),  lower  frequencies  are  more  commonly  used.  A  lower  limit  to  the  fre¬ 
quency  used  is  set  by  the  resolution  desired  --  layers  thinner  than  one 
half  the  wavelength  generally  cannot  be  resolved  clearly.  Pulse  duration 
is  another  important  factor  and  must  be  short  for  good  resolution.  In 
general , 

2d  =  vAt,  (15) 

where 

d  =  distance  between  reflecting  layers  (m), 

v  =  velocity  of  sound  in  water  (1500  m/sec),  and 

At  =  difference  between  arrival  times  of  the  reflected  wave  fronts. 

Thus,  when  the  pulse  length  becomes  comparable  with  the  time  of  2d/v  sec, 
the  two  layers  may  not  be  resolved. 

Each  type  of  reflection  equipment,  which  differs  principally  in  the 
type  of  energy  source  utilized,  has  its  own  uses.  For  a  comprehensive  dis¬ 
cussion  of  subbottom  systems,  see  Schlank  (1968).  Briefly,  a  transducer 
system  operating  at  frequencies  of  1.5  to  12  kHz  can  show,  under  favor¬ 
able  circumstances  (shallow  water,  silt  and  clay  sediments),  a  subbottom 
structure  to  depths  of  50  m,  althougn  penetration  of  20  m  is  more  common. 
This  penetration  is  more  than  adequate  to  satisfy  the  objectives  for  de¬ 
termining  surface  sediment  properties  as  done  in  this  Investigation. 

Problems  with  fixed  frequencies  and  various  sediment  types  are  as 
those  noted  by  Moore  and  Kulm  (1970),  who  describe  a  particular  case: 

"No  sub-bottom  information  was  obtained  with  the  3.5 
kHz  profiler  on  the  Oregon  continental  shelf.  Cores 
taken  off  Oregon  indicate  that  thick  sand  layers  cover 
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much  op  the  continental  shelf  and  are  overlain  by 
only  a  thin  layer  of  mud  (a  few  cm  thick)  on  the 
middle  and  outer  shelf.  Apparently  the  signal 
strength  of  the  3.5  kHz  system  is  not  sufficient 
to  penetrate  these  sandy  deposits." 

Our  experience  with  3.5  kHz,  4.5  kHz,  and  12  kHz  sources  also  show  little 

*  or  no  penetration  in  similar  sediment  even  at  the  shallow  depths  of  area 
surveyed. 

Discounting  the  need  for  subsurface  structural  data  and  consider¬ 
ing  only  the  surface  one  can  say  that  for  reflectivity  measurements  in¬ 
volving  the  sea-floor  surface  material  only,  it  seems  that  any  sound 
source  may  be  used  as  long  as  the  measurement  of  the  pulses  is  proper  for 
the  sound  source  being  used;  however,  for  sub-bottom  observations  it  ap¬ 
pears  that  basically  a  low-frequency  with  broad-band  characteristics  is 
needed  to  allow  both  penetration  and  resolution  of  the  sub-bottom  layers. 
Regarding  the  high  frequency  sources  (>  4  kHz),  it  has  been  acknowledged 
that  standard  echo-sounders  can  offer  simple  tools  for  the  geologist  to 
use,  and  with  careful  interpretation  can  provide  a  wealth  of  data  impos¬ 
sible  to  obtain  as  simply  by  any  other  means.  Any  conclusion  as  to  the 
presence  of  a  certain  sediment  tvpe  should  only  be  made  after  confirming 
evidence,  such  as  by  cores,  grab  samples,  or  photographs  examined.  The 
experiments  conducted  in  this  study,  and  in  previous  work  (Barnes  et 
al . ,  1972)  with  the  12  kHz  and  41  kHz  sources  tend  to  corroborate  this. 

Two  sources  were  used  in  the  reflection  coefficient  mapping  experi¬ 
ment  described  in  this  report.  They  consist  of  a  standard  12  kHz  trans¬ 
ducer  and  a  Directional  Multi-Electrode  Sparker  Sound  Source  (DMSS)  de¬ 
signed  and  built  in  this  laboratory,  and  described  in  Barnes  et  al.  (1972) 
and  by  Poston  (1972).  The  two  sources  were  used  simultaneously  in 
^  tests.  The  12  kHz  pinger  transducer  produced  a  pulse  length  of  1,625 

m/sec,  and  the  multi-electrode  sparker  sound  ;ource  transmitted  a  nominal 
pulse  length  ranging  from  .3  to  .5  m/sec.  The  higher  frequency  source 

*  was  discrete  in  bandwidth,  but  the  low  frequency  source  was  composed  of 
a  nominal  2.5  kHz  center  frequency  with  a  broad-band  width  with  frequen¬ 
cies  ranging  from  0  to  5  kHz. 
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Ten  of  the  twenty-five  available  electrodes  on  the  DMSS  sparker 
were  used  during  the  experiment  (fig.  12).  This  number  was  selected  af¬ 
ter  extensive  testing  of  impulse  output  versus  bubble  pulse  output  for 
various  numbers  of  grouped  electrodes.  For  this  experiment,  300  joules 
or  800  joules  of  output  could  be  selected  at  the  capacitor  bank.  800 
jouies  was  used  during  the  experiment.  Several  modifications  to  the  pre¬ 
vious  system,  described  in  Barnes  et  al .  (1972),  have  been  incorporated 
in  this  work.  The  sparker  is  mounted  on  a  surface  towed  vehicle  fabri¬ 
cated  from  a  metal  frame  and  2  surfboards  (fig.  13).  Adjustment  of  the 
mounting  bracket  and  DMSS  unit  makes  it  possible  to  orient  the  electrodes 
vertically  downward  into  the  water  or  upward  toward  the  air/sea  inter¬ 
face.  An  upward  orientation  proved  better,  since  there  was  less  bubble 
pulse  and  therefore  less  ambient  noise  interference.  Vertical  adjustment 
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Figure  12.  Directional  Multiple  Electrode  Sound  Source  (DMSS). 
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Figure  1*.  DMSS  Unit  and  Impulse  Hydrophone 
'deceiver  —  Tow  Configuration. 


of  the  source  also  enabled  the  position  of  the  sparker  to  be  varied  rela¬ 
tive  to  depth  below  the  water  surface.  Tests  conducted  on  various  other 
seismic  sources  (i.e.,  the  Uni  Boom,  the  3-electrode  sparker  manufactured 
by  Egerton  Germershausen  and  Grier,  Inc.,  and  a  12  kHz  transducer  manu¬ 
factured  by  Ocean  Sonics,  Inc.)  and  other  non  spark-type  discharge  sources 
in  this  geometrical  configuration  indicated  that  the  quality  of  the  trans- 
mi  oted  pulse  could  be  improved  if  the  source  was  moved  closer  to  the 
water  surface.  This  modification  in  geometry  differs  from  the  position 
used  in  the  previous  experiment  conducted  by  Barnes  et  al .  (1972).  A 
marked  improvement  in  the  quality  of  seismic  subbottom  profiles  was  noted 
when  the  position  of  the  electrode  tips  was  set  at  about  10  inches  below 
the  water  surface.  In  most  instance,  it  was  discovered  that  much  more 
information  could  be  obtained  from  the  broadband  frequency  source.  The 
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essential  requirements  of  any  acoustic  source  employed  for  reflection  co¬ 
efficient  mapping  is  that  it  should  provide  a  broad  frequency  band  acous¬ 
tic  pulse  and  be  capable  of  analysing  any  changes  produced  in  this  pulse 
by  the  effect  of  the  various  bottom  and  sub-bottom  reflectors. 

Tests  of  the  broad  band  multi -electrode  sparker  revealed  much  the 
same  results  as  obtained  by  Li  and  Taylor-Smith  (1969).  Their  comments 
about  experimenting  with  a  multi -electrode  sparker,  were,  "as  the  number 
of  sparks  in  parallel  increased  so  did  the  efficiency  of  producing  acous¬ 
tic  energy."  In  the  tests  conducted  by  this  laboratory,  the  energy  peaked 
and  then  decreased  after  the  number  of  electrodes  was  increased  at  a  given 
power  level.  Also,  the  duration  of  the  bubble  forming  period  and  pulse 
length  decreased,  leading  to  a  greater  concentration  of  energy  in  the 
higher  frequencies.  The  pulse  length  and  frequency  content  was  also 
found  to  be  dependent  on  the  spark  gap  distance.  The  DMSS  gap  was  norm¬ 
ally  held  constant,  except  when  the  electrodes  burned  too  short  and  arced 
to  the  base  of  the  sparker  unit,  producing  a  higher  frequency  output. 

The  two  sources  were  used  simultaneously  and  were  triggered  at  a 
1  sec  rate  by  the  recorder.  The  time  base  stability  of  the  recorder  is 
specified  to  be  10  parts  per  million.  It  was  intended  that  the  12  kHz 
plnger  source  be  triggered  at  a  .25  sec  rate  and  the  sparker  at  a  1  sec 
rate  with  the  use  of  an  external  countdown  circuit  and  relay  to  trigger 
the  capacitor  bank,  but  the  triggering  operation  was  too  erratic.  The 
addition  of  a  separate  power  supply  and  a  one-shot  multi  vibration  circuit 
did  not  produce  good  results.  After  many  trials  the  problem  was  found  in 
the  trigger  circuit  of  the  recorder.  Random  pulses  caused  the  two  sources 
to  be  triggered  out  of  synchronization.  A  preliminary  study  of  the  repe¬ 
tition  rate  was  done  on  the  first  tapes  made,  and  it  was  discovered  that 
a  steady  rate  of  1  sec  repetition  was  more  desirable  for  computer  process¬ 
ing  efficiency. 

4.1.2  Hydrophone  Receivers 

There  were  three  hydrophones  used  in  conjunction  with  the  sparker 
signal.  The  hydrophones  monitored  the  transmitted  and  reflected  signals. 
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The  sled -mounted  hydrophone  were  used  to  receive  the  transmitted  impulses 
from  the  sparker.  The  two  return  impulse  hydrophones  used  during  the  ex¬ 
periment  were  the  streamer  hydrophone  and  a  second  sled-mounted  hydrophone, 
not  shown  in  figure  13.  The  receiver  hydrophones  were  identical  in  elec¬ 
tronics,  but  of  different  lengths.  The  streamer  hydrophone  was  approxi¬ 
mately  11  ft  long  with  a  fluid  filled  section  of  10  ft,  8J5  in.  The  sled- 
mounted  hydrophone  was  approximately  6  ft  long  with  a  fluid  filled  section 
of  62  in.  The  hydrophones  including  the  impulse  unit  measured  1-5/8  in 
in  diameter.  Six  hydrophone  elements  make  up  each  of  the  streamers.  These 
are  wired  in  parallel  to  a  preamplifier  and  analog  recording  devices.  The 
six  hydrophone  elements  are  summed  for  a  gain  of  15.6  dB,  which,  added  to 
the  basic  hydrophone  element  sensitivity  of  96  dB/u/ubar,  yields  a  basic 
hydrophone  sensitivity  of  60.4  dB/u/ubar.  The  hydrophone  arrays  were  cal¬ 
ibrated  using  a  Massa  Model  115C  calibrated  hydrophone  and  a  Navy  J-9 
sound  projector.  The  measured  test  values  were  -60.5  dB/u/ubar. 

The  sled  mounted  return  impulse  hydrophone  did  not  produce  useable 
data  during  the  test  due  to  excessive  ringing  and  noise,  so  all  the  pro¬ 
cessed  data  was  obtained  from  the  streamer  hydrophone.  Various  attempts 
were  made  to  remedy  the  problem,  including  the  use  cf  a  less  sensitive 
hydrophone,  but  this  resulted  in  insufficient  amplitude  to  tape  the  re¬ 
turn  signal. 

All  three  hydrophones  are  connected  to  the  same  amplifier  chassis; 
however,  each  hydrophone  has  its  separate  amplifier  circuit  and  amplitude 
control.  The  streamer  and  return  impulse  hydrophone  signals  were  at  unity 
gain,  the  only  amplifier  being  the  preamplifier  in  the  hydrophone  array. 

The  average  effective  dynamic  range  of  the  hydrophone  is  35.4  dB.  The 
low-end  sensitivity  is  limited  by  the  towing  noise  observed  on  the  records. 
The  nominal  operating  dynamic  range  is  66  dB  when  the  low  end  sensitivity 
is  limited  only  by  ambient  noise  generated  by  the  electronics.  The  streamer 
hydrophone  output  was  also  connected  to  the  Ocean  Sonics  recorder  through 
the  rack  mounted  amplifier  which  has  a  gain  of  67  or  37  dB. 

The  location  and  position  of  the  monitoring  hydrophones  was  also 
improved  by  modification  of  the  tow  vehicle.  A  more  sensitive  receiver 
hydrophone  was  installed  in  addition  to  the  standard  unit  used  in  previous 
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experiments  (Barnes  et  al.,  1972).  The  receiving  and  impulse  hydrophones 
were  placed  27  in  and  32  in  below  the  water  surface,  respectively. 

In  conjunction  with  the  sparker,  an  EDO  Type  324  (AT/200  A-UQN) 
transducer  array  using  A.D.P.  type  crystals  was  used.  The  12  kHz  trans¬ 
ducer,  referred  to  as  the  pinger  source,  was  driven  by  an  Ocean  Sonics 
Model  GDR-T  precision  depth  recorder.  The  transmitted  pulse  is  rated  at 
600  w  peak  power  into  the  rated  130  ohm  load  of  the  transducer,  which 
yields  an  acoustical  power  output  of  110  dB/v  at  one  yard.  The  trans¬ 
ducer  has  a  beam  angle  of  35°  on  the  frontal  plane.  The  receiving  re¬ 
sponse  is  -72  dB/v  at  one  yard  (EDO  Western  Corp.). 

The  pinger  is  mounted  on  the  end  of  a  two  in  diameter  pipe  which 
slides  up  and  down  in  a  socket  at  the  end  of  an  A-frame  on  the  port  side 
of  the  boat.  The  pipe  is  locked  into  position  (fig.  14).  The  depth  of 
the  transducer  was  set  at  1.8  ft  below  the  water  surface.  The  distance 
from  the  side  of  the  boat  to  the  center  of  the  pinger  was  29  in,  and  the 
distance  from  the  center  of  the  pinger  to  the  stern  was  16h  ft.  This  lo¬ 
cation  was  picked  because  of  the  wake  characteristics  of  the  boat  and  con¬ 
venience  of  mounting.  The  A-frame  was  mounted  to  the  deck  with  pillow 
blocks,  allowing  the  transducer  to  be  raised  up  and  Inboard,  out  of  the 
way  for  portside  docking  and  fast  running. 

4,1.3  Transmission  Measurements 

During  field  tests,  the  transmission  impulses  and  return  impulses 
were  constantly  monitored  on  an  oscilloscope.  The  pulses  were  visually 
compared  to  one  another  for  the  reflection  measurement  and  to  observe 
signal/noise  ratios.  This  monitoring  is  very  important  since  the  ampli¬ 
tude  of  the  return  pulse  is  dependent  upon:  (1)  The  transmitted  ampli¬ 
tude  of  the  acoustic  source;  (2)  the  path  length  of  the  sound  signal  to 
and  from  the  reflecting  sea-floor;  and  (3)  the  reflecting  properties  of 
the  sea-floor. 

Because  the  amount  of  sound  energy  received  is  dependent  on  the 
amount  of  energy  transmitted,  and  mor.t  echo  sounders  have  no  method  of 
displaying  the  relative  energy  levels,  one  cannot  assume  that  the  trans¬ 
mitted  energy  is  constant.  A  variation  in  the  power  supply  voltage,  for 


instance,  can  cause  a  considerable  variation  in  the  electrical  power  sup¬ 
plied  to  the  transducer.  Also,  relationships  existing  between  the  physi¬ 
cal  properties  of  the  transducer  and  those  of  the  surrounding  water,  such 
as  salinity  will  affect  the  output  acoustic  power  even  if  the  electrical 
power  is  constant.  This  is  especially  true  of  a  sparker  device.  This  has 
been  calculated  by  Caulfield  (1962,  p.  342),  i.e.. 


=  R  B2  (  k2  ) 

Total  sp  2F  ^2.^2) 


(16) 


where 

^Total  =  Total  ener9y  t0  spark, 

Rsp  =  sPark  resistance, 

B2  =  V2/L2K2, 

b  =  R/2L , 

k  =  (R2/4L2  -  1/LC 

R  =  Total  resistance, 

L  =  Total  inductance,  and 

C  =  Total  storage  capacitance. 

Equation  (16)  describes  what  occurs  in  the  electrical  circuit.  For 

example,  the  spark  resistance  (R  )  can  be  affected  by  the  salinity  of  the 

sp 

water,  which  can  vary  from  3.9  mho/cm  for  a  2  percent  salinity  solution  to 
43.9  mho/cm  for  a  35  percent  salinity  solution. 

If  one  desires  to  know  acoustical  pressure,  the  equation  approxi¬ 
mating  the  first  pressure  pulse  as  a  function  of  time  is, 


P  =  Poe"mt  (17) 

where 

p  =  pressure  (dyne/cm2), 

p0  =  peak  pressure  (dyne/ cm3  at  1  meter), 

m  =  decay  constant  (determined  by  experiment),  and 

t  =  time. 
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Then  the  acoustic  energy  is  expressed  as 


lf(+)l2  dt, 


(18) 


where 

f(t)  =  p0e"mt, 
p  =  density  (g/cm3), 

c  =  velocity  of  propagation  in  the  medium,  and 
m  =  decay  constant  (determined  by  experiment). 

Caulfield's  formula  can  be  correlated  to  Breslau's  (1967,  p.  2) 
formula  in  terms  of  peak  pressure,  i.e.. 


i 

P*  dt  =  P*  (RMS)  t  (19) 

where 

Pg  =  pressure  of  the  source  pulse, 

t  =  duration  of  the  outgoing  pulse,  and 

T  =  time. 

The  formula  used  to  obtain  the  values  of  the  Rayleigh  Reflection 
Coefficient  from  the  digitized  values  of  V<-^  +  Vsx  may  be  expressed  as; 


(see  Appendix  A-l), 


i  v  cn 

R  =  .008  (4d2  +  832)  2  r~ 

VSX 


(20) 


where 

R 

=  Rayleigh  Reflection  Coefficient, 

d 

=  water  depth, 

VSR 

=  sparker  return  impulse  voltage,  and 

VSX 

=  sparker  transmitted  impulse  voltage. 
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The  12  kHz  pinger  transducer  transmitted  output  wcs  also  constantly 
monitored  during  field  operations.  The  voltage  applied  co  the  transducer 
was  monitored  on  an  oscilloscope  while  being  recorded  cr  an  analog  mag¬ 
netic  tape  recorder.  A  fixed  resistive  attenuator  network  to  reduce  the 
high  transmitted  voltage  to  a  lower  voltage  was  necessary  to  make  the  re¬ 
cording  level  suitable  for  input  to  the  tape  recorder.  The  return  Impulse 
voltage  was  amplified  by  the  internal  amplifier  of  the  recorder,  (section 
4.1.2)  and  then  connected  to  the  magnetic  tape  recorder. 

Many  problems  were  encountered  in  calibrating  the  pinger  trans¬ 
ducer  due  to  the  small  physical  size  of  the  tanks  available  for  calibra¬ 
tion  use  in  the  laboratory,  and  the  large  amplitude  of  the  transmitted 
signal.  Proper  interrogation  of  the  pulse  to  be  used  for  the  computation 
of  the  reflection  coefficient  was  threefor®  difficult.  Problems  were  en¬ 
countered  in  obtaining  repetitive  measurements  using  a  calibrated  hydro¬ 
phone  and  a  sound  source  In  the  laboratory  tanks.  The  calibrating  hydro¬ 
phone  and  source  are  linear  devices,  and  the  pinger  transducer  Is  a  tuned, 
or  non-linear  device;  the  ratio  of  the  transmitted  to  received  voltages 
of  the  transducer  Itself  could  have  been  used  to  obtain  an  absolute  cali¬ 
bration  value,  but  physical  limitations  of  calibrating  in  the  tank  re¬ 
sulted  In  relying  on  manufacturer's  specifications.  Data  for  obtaining 
absolute  values  of  the  reflection  coefficient  were  taken  from  the  manu¬ 
facturer's  calibration  specifications. 

Output  variations  based  on  reflectivity  measurements  have  been  ob¬ 
served  by  others  employing  various  kinds  of  sources.  Li  and  Taylor-Smith 
(1969,  pp.  244,  245)  noted  intensity  fluctuations  using  a  sparker  source. 
In  their  examination  of  the  intensities  of  sound  returned  to  the  surface, 
both  from  bottom  and  sub-bottom  reflections,  it  was  thought  that  the  loss 
in  intensity  might  be  attributed  to  a  large  number  of  factors,  of  which 
absorption  in  the  medium  is  only  one.  They  go  on  to  say  that,  the  effect 
of  the  unwanted  losses  may  be  reduced  from  a  statistical  point  of  view, 
by  a  procedure  used  In  spectral  analysis,  in  so  doing,  It  Is  essential 
that  each  set  of  pulses  analysed  is  taken  from  single  individual  transmis¬ 
sions.  However,  even  with  this  procedure,  a  fluctuation  In  the  intensity 
of  the  reflection  can  be  caused  both  by  the  irregularity  of  the  reflecting 


boundaries  and  by  the  irregular  structural  pattern  of  the  medium  between 
the  bottom  and  sub-bottom.  The  effect  is  mainly  one  of  random  scattering, 
although  interference  effects  probably  play  an  important  part  as  well. 

Such  fluctuations  in  intensity  may  be  of  diagnostic  value. 

A  statistical  study  of  various  reflected  signals  shows  that  the 
observed  fluctuations  appear  to  correlate  with  the  state  of  roughness  of 
the  reflector  (section  3.1,  fig.  3).  Li  and  Taylor-Smith  (1969)  found 
that,  the  standard  deviation  of  the  intensity  fluctuations,  expressed 
as  a  percentage  of  the  mean  peak  to  peak  amplitude,  for  the  signal  trans¬ 
mitted  directly  through  the  water  from  a  3-electrode  spark  is  about  4 
percent,  whereas  the  figure  obtained  for  a  reflection  from  a  sand-covered 
floor  is  about  10  percent  and  from  a  clay-covered  floor  5  percent.  While 
these  are  isolated  results  from  which  general  conclusions  cannot  be  drawn, 
it  is  certain  that  scattering  due  to  boundary  roughness  is  an  important 
factor  in  signal  fluctuations.  In  fact,  this  property  is  employed  in  the 
attempt  to  use  the  oblique  asdic  (or  sideways-looking  sonar)  to  define 
the  sedimentary  state  of  the  sea-floor  (Chesterman  et  al . ,  1967). 

Intensity  fluctuations  were  also  noted  by  Porter  and  Bell  (1972) 
in  their  reflection  studies  with  3.5  krlz  and  5  kHz  electromechanical 
transducers.  A  bottom  variation  of  ±  0.85  dB  was  noted  in  this  experi¬ 
ment,  which  was  conducted  with  the  source  mounted  in  a  rigid  frame  on 
the  bottom  of  an  inland  body  of  water.  They  found  that  "the  magnitude 
of  the  observed  variations  differed  significantly  and  exceeded  the  val¬ 
ues  calculated  for  the  component  attributable  to  source-receiver  motion," 
and  "concluded  that  other  error  sources  such  as  changes  in  incident  an¬ 
gle  and  undetected  spatial  variability  in  sediment  composition  exist." 

They  did  note  that  a  reduced  degree  of  variability  was  exhibited  in  the 
narrow-band  data,  which  they  attributed  to  an  inherently  higher  signal - 
to-noise  ratio. 

4.1.4  Data  Acquisition  and  Display 

The  shipboard  recording  equipment  consists  of  a  rack-mounted  Am- 
pex  SP-300,  7-channel  magnetic  tape  recorder.  The  inputs  to  the  seven 
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channels  are: 

Channel  1:  Pinger  transmitted  voltage  (fig.  15a) 

Channel  2:  Sparker  received  impulse  (Streamer  Hydro¬ 
phone)  (fig.  16a) 

Channel  3:  Pinger  received  impulse  (fig.  15b) 

Channel  4:  Voice 

Channel  5:  Sparker  transmitted  Impulse  (fig.  16b) 

Channel  6:  Time  Code,  IRIG  B 

Channel  7:  Sparker  received  Impulse  (sled  hydrophone). 

An  Ocean  Sonlcs  Series  GDR-T  precision  Sonar  recorder  and  transceiver  was 
used  as  a  visual  recorder  for  the  subbottom  profile  record  (fig.  17).  The 
recorder  was  modified  to  allow  the  presentation  of  either  the  12  kHz 
pinger  or  the  0.5  kHz  sparker  seismic  returns.  The  recorder  unit  also 
houses  the  energy  source  for  the  12  kHz  pinger  transudcer.  The  output 
voltage  of  the  transmitter  in  the  recorder  to  a  35  dB  attenuator  network 
to  reduce  the  output  voltage  to  a  suitable  taping  level. 

The  recorder  is  also  modified  so  that  the  amplified  output  of  the 
received  pinger  Impulse  Is  connected  to  the  tape  recorder.  The  amplifier 
gain  can  be  changed  from  a  setting  of  XI,  which  produced  a  measured  gain 
of  50,  to  a  setting  of  X50,  which  produced  a  measured  gain  of  920.  tlec- 
trlcal  cables  from  the  streamer,  and  sled-mounted  amplifier.  The  output 
of  the  streamer  hydrophone  was  fed  from  the  hydrophone  preamplifier  to 
the  tape  recorder.  The  streamer  hydrophone  output  was  also  connected 
to  a  Khron-HIte  Model  3200  Variable  High  or  Low  Pass  electronic  filter 
and  then  to  the  main  hydrophone  amplifier.  This  amplified  output  was 
connected  to  the  Ocean  Sonics  recorder  to  observe  the  2.5  kHz  sparker 
subbottom  Information.  The  output  from  the  receiver  hydrophone  was  fed 
through  an  attenuator  network  to  the  tape  recorder.  During  the  survey 
the  attenuation  was  set  at  unity.  The  sparker  impulse  hydrophor.a  outp 't 
was  fed  through  an  amplifier  to  the  tape  recorder.  Throughout  the  survey 
the  amplifier  gain  was  set  at  unity.  The  voice  channel  was  used  to  note 
the  same  events  noted  in  the  tape  roll  logs.  The  time  reference  for  the 
magnetic  tape  record  was  obtained  from  an  IRIG  B  Time  Code  generator. 

The  generator  and  the  digital  clock  used  in  the  navigation  system  were 


a.  0-5  kHz  transmitted  signal 
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b.  0-5  kHz  received  signal. 


Figure  16.  Sparker  signal  wave  forms. 
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Figure  17.  Typical  3ubbottom  profile  record  —  Monterey  Bay 

test  area. 


set  together  in  reference  to  WWV.  The  timing  lines  on  the  Ocean  Sonics 
recorder  were  also  controlled  by  the  same  digital  clock  marker  interface 
so  that  the  locations  of  the  seismic  records  could  be  accurately  deter¬ 
mined. 


4.2  Shear  Wave  Experiment 


4.2.1  Shear  Wave  Generator 


The  seismic  source  needed  for  the  shear  wave  experiment  had  to 
fulfill  three  important  requirements.  It  had  to  (1)  be  rich  in  shear 
energy,  (2)  generate  waves  of  low  frequency  (as  attenuation  is  propor¬ 
tional  to  frequency),  and  (3)  discharge  its  energy  below  the  water-sedi 
ment  interface.  Several  sources  were  considered  as  to  their  capability 
of  fulfilling  the  above  requirements.  A  3  cubic  inch  air  gun  was  the 
final  choice  as  a  shear-wave  generator. 


Attached  to  the  air  gun  is  a  probe,  a  hollow  1  in  diameter  pipe 
with  a  coned  point.  The  energy  from  the  air  gun  is  transported  through 
the  probe  and  then  out  through  an  exit  port  at  the  base  of  the  probe 
(fig.  18).  By  implanting  the  0.5  m  long  probe  into  sediment,  the  source 
energy  is  discharged  below  the  water-sediment  interface. 

The  exit  port  is  located  on  the  side  of  the  probe  to  generate  shear 
waves  (in  particular,  S^-waves).  With  the  probe  implanted  into  the  sedi¬ 
ment,  water  and  a  little  sediment  fill  the  hollow  probe.  When  the  gun 
is  fired,  the  water-sediment  mixture  is  forced  out,  producing  a  horizon¬ 
tal  jet  stream.  Thus,  the  basic  procedure  to  detect  the  direct  S^-wave 
involves  a  linear  array  of  transverse  geophones  with  the  S^-source  placed 
at  one  end  of  the  detector  array.  The  sensitive  axis  of  the  geophones 
and  the  water  scream  from  the  probe  are  parallel  and  both  are  transverse 
to  the  direct  S^-ray  path.  This  alignment  of  horizontal  geophones  and 
SH-source  produces  maximum  response  to  the  direct  SH-wave  (fig.  19). 

The  3  cubic  inch  air  gun  (Rix  Industries,  Emeryville,  California) 
is  cylindrically  shaped  with  a  height  of  5  in  and  diameter  of  3  in.  The 
high  pressure  air  escapes  from  a  large  diameter  single  opening  which 
provides  for  a  concentration  acoustical  source.  High  pressure  is  sup¬ 
plied  to  the  air  gun  from  a  single  Scuba  diver's  tank.  Employing  a  Scuba 
tank  provides  a  high  pressure  range  from  300  to  2200  psi;  however,  a  high 
pressure  setting  of  600-1000  psi  is  used  for  tests.  The  air  gun  is  re¬ 
motely  triggered  by  a  switch  which  activates  a  solenoid,  which  in  turn 
discharges  the  high  pressure  air  in  the  piston  chamber  of  the  air  gun. 

There  are  two  inherent  disadvantages  concerning  this  S^-source. 

One  obvious  problem  is  sediment  disturbance.  Although  the  water  stream 
from  the  probe  is  directed  away  from  the  sediment  between  the  source  and 
receivers,  the  air  discharged  into  the  sediment  will  eventually  produce 
a  cavity  after  several  firings.  This  reduces  the  efficiency  of  produc¬ 
ing  SH~waves.  The  second  disadvantage  arises  when  the  piston  of  the  air 
gun  returns  after  discharging  the  air.  This  draws  excessive  amounts  of 
sediment  into  the  probe,  causing  a  blockage.  The  latter  problem  was 
rectified  by  attaching  a  fine  wire  mesh  across  the  exit  port,  allowing 
only  fines  to  enter  the  probe. 
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Figure  18.  Marine  S ^-source. 
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Figure  19.  Probe  for  measurement  of  shear -wave  and 
compress ional-wave  velocity. 


Altogether,  the  source  proved  to  be  successful  in  producing  low 
frequency  shear  waves  in  the  frequency  ranges  from  80-100  Hz. 

4.2.2  Receiver  Array 

Four  geophone  units  are  employed  to  record  the  arrival  of  the 
direct  shear  wave  and  direct  congressional  wave  (denoted  hereafter  as 
direct  SH-wave  and  direct  p-wave,  respectively).  Each  unit  contains  a 
vertical  geophone  to  respond  to  the  direct  p-wave  and  a  horizontal 
geophone  to  respond  to  the  direct  S^-wave  (receivers  manufactured  by 
Geospace,  Houston,  Texas;  HS-5  Model  K,  standard  frequency  of  30  Hz). 
Each  geophone  unit  was  equipped  with  a  1  in  brass  stake  used  to  im¬ 
plant  units. 
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4.2.3  Signal  Processing  and  Display 


The  four  horizontal  geophone  outputs  are  fed  into  a  four  channel 
voltage  gain  amplifier.  The  amplified  signals  are  then  displayed  on  a 
dual  beam  oscilloscope  which  was  converted  to  four  independent  traces 
with  a  Tektronix  1A-4  plug-in.  The  sweep  on  the  oscilloscope  is  trig¬ 
gered  from  the  same  circuit  that  triggers  the  solenoid  of  the  air  gun. 
Unfortunately,  the  air  gun  has  a  52  m/sec  delay  time  between  the  activa¬ 
tion  of  source  trigger  switch  and  the  actual  firing  of  the  air  gun.  Thus, 
it  became  necessary  to  delay  the  sweep  on  the  oscilloscope  in  order  that 
the  signals  cor i d  be  viewed  on  the  oscilloscope  screen.  The  signals  are 
recorded  by  means  of  a  single  exposure  oscilloscope  camera  equipped  with 
Polaroid  type  47  film,  3000  ASA  speed.  The  process  is  repeated  by  feed¬ 
ing  the  four  vertical  geophone  outputs  into  the  amplifier  and  then  re¬ 
cording  the  signals  displayed  on  the  oscilloscope  with  the  camera.  A 
functional  block  diagram  of  this  simple  system  is  shown  in  figure  20. 
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4.3  Electrical  Resistivity  System 
4.3.1  Expandable  Array 

The  system  used  to  measure  bottom  resistivity  was  designed  for 
operation  by  two  Scuba-equipped  divers,  and  consists  of  separate  current 
and  voltage  units,  a  length  of  line  on  which  electrode  separations  are 
marked,  and  an  underwater  slate  for  recording  data  (fig.  21).  The  out¬ 
put  of  the  current  source  is  a  10  Hz  square  wave,  with  amplitude  adjus¬ 
table  between  0  and  500  milliamps  in  10  steps.  Power  is  supplied  by 
16  D-size  rechargeable  nickel -cadmium  batteries,  in  a  series-parallel 
arrangement  which  provides  a  nominal  10  V  output.  Current  electrodes 
are  1/4  in  (0.64  cm)  diameter  copper  rods,  1/2  in  (1.27  cm)  long. 

The  voltage  unit  is  a  high-impedance  (106  megohms)  rectifying 
voltmeter,  with  ranges  of  1,  10,  100,  1000,  and  10,000  MV  full  scale,  and 
is  powered  by  mercury  batteries.  The  potential  electrodes  have  a  non¬ 
polarizing  silver-silver  chloride  element,  housed  in  a  1  in  (2.54  cm) 
outside  diameter  cylinder  of  acrylic  plastic,  and  are  filled  with  2.7 
molar  potassium  chloride  saturated  with  silver  chloride.  Electrical 


Figure  21.  Schematic  diagram  -  diver  operated  direct 
current  resistivity  apparatus. 


contact  to  the  sea  water  is  made  through  a  rod  of  porous  ceramic  material 
(Corwin,  1973).  Each  range  of  the  voltmeter  is  calibrated  by  measuring 
a  known  voltage  across  a  fixed  resistor. 

Both  current  and  voltage  units  are  housed  in  clear  acrylic  plastic 
cylinders  of  6  in  (15.2  cm)  inside  diameter  and  1/4  in  (0.64)  wall  thick¬ 
ness,  which  allows  direct  reading  of  the  current  and  voltage  meters.  All 
openings  are  sealed  with  0-rings,  and  the  housings  are  pressurized  to  a 
few  psi  before  immersion  to  check  for  leaks.  The  length  of  the  current 
unit  is  12  in  (30.5  cm);  the  voltage  unit,  9H  in  (24.1  cm).  Sufficient 
cable  is  provided  for  an  electrode  separation  a  of  300  cm  in  the  Wenner 
configuration . 

In  use,  a  line  marked  for  Wenner  electrode  spacings  of  10,  30, 

100,  and  300  cm  is  laid  on  the  bottom  with  the  ends  anchored  by  diving 
weights.  The  current  and  voltage  electrodes  then  are  inserted  into  the 
sediment  at  the  locations  for  a  =  10  cm  and  the  value  of  the  voltage, 

AV,  recorded  for  several  values  of  the  current,  I.  The  electrodes  then 
are  moved  to  the  next  larger  electrode  spacing  and  the  procedure  repeated. 
As  skin  contact  resistance  is  low  in  salt  water,  care  is  taken  by  the  di¬ 
vers  not  to  touch  the  current  electrodes  during  operation,  to  avoid  the 
possibility  of  a  dangerous  shock. 

Electrode  spacing,  current,  voltage,  and  water  depth  and  tempera¬ 
ture  are  recorded  on  the  slate,  and  a  water  sample  is  taken  for  later 
salinity  measurement.  Although  slow  and  cumbersome  to  use,  especially 
in  water  of  poor  visibility,  this  prototype  system  proved  reliable  in 
field  operation,  and  yielded  data  which  appears  to  confirm  the  feasi¬ 
bility  of  the  sea-floor  horizontal  array  concept. 

4.3.2  Surface  Towed  Array 

The  equipment  used  for  the  towed  resistivity  array  field  test  is 
shown  schematically  in  figure  22.  A  12  V  lead  acid  automotive  storage 
battery  was  used  as  the  current  source,  with  current  I  controlled  by  the 
rehostat  and  monitored  on  a  0-10  amp  ammeter.  The  switch  allows  the  cur¬ 
rent  to  be  turned  off  to  monitor  the  zero  level  of  the  potential  elec¬ 
trodes,  or  for  polarity  reversal.  The  current  electrodes  were  strips  of 


zinc-coated  pipe  hanger  strap,  1  in  (2.54  cm)  wide  by  about  2  ft  (60  cm) 
long.  The  potential  electrodes  were  of  the  silver-silver  chloride  type 
described  above,  and,  potential  aV  was  read  on  an  Esterline-Angus  model 
T171B  battery  operated  strip  chart  recorder. 

In  operation,  the  current  electrodes  are  energized  continuously 
over  areas  of  interest,  and  are  turned  off  for  about  30  sec  every  5  min 
to  check  the  zero  level  of  the  potential  electrodes.  Typical  chart  re¬ 
corder  output  is  shown  in  figure  23.  Current  polarity  is  reversed  every 
so  often  to  equalize  electrolytic  dissolution  of  the  current  electrodes. 
Although  continual  reversal  of  current  electrode  polarity  would  have,  in 
effect,  doubled  the  signal-to-noise  ratio,  it  was  found  that  this  proce¬ 
dure  resulted  in  a  degraded  current  waveform,  with  several  minutes  re¬ 
quired  for  the  current  to  reach  a  stable  value.  As  the  current  was  not 
recorded  continuously  but  was  read  from  the  meter  and  noted  on  the  po¬ 
tential  record,  it  was  decided  to  hold  current  polarity  constant  as  des¬ 
cribed  above. 

4.4  Sediment  Sampling  Equipment 

To  verify  reflectivity  measurements,  it  was  necessary  to  have  a 
core  sampler  which  could  rapidly  sample  the  first  meter  of  surface  sedi¬ 
ments  and  obtain  relatively  undisturbed  samples.  This  was  required  for 
reliable  engineering  and  mass  physical  property  determinations.  A  core 


Figure  22.  Schematic  of  towed  array  —  direct 
current  resistivity  system. 
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Figure  23.  Towed  resistivity  —  chart  recorder  output. 


length  ranging  from  1  m  to  3  m  was  decided  upon  because  this  depth  of 
penetration  was  within  the  wave  length  of  either  of  the  two  acoustic 
sources  employed  for  the  reflection  coefficient  mapping  experiment  (sec¬ 
tion  4.1.1).  This  also  made  it  possible  to  determine  to  what  degree 
changes  in  the  near  surface  lithology  or  small  scale  structural  varia¬ 
tions  would  influence  the  reflection  coefficient.  Several  samplers  of 
the  gravity-type  or  piston-type  gravity  corer  variety  were  tested  but 
found  unsatisfactory.  It  was  therefore  necessary  to  design,  construct, 
and  test  coring  devices  that  would  provide  the  desired  capability.  It 
was  also  decided  that  the  energy  source  for  penetration  should  not  be 
percussion  or  rotary  since  these  types  tend  to  induce  sediment  distor¬ 
tion.  Therefore,  a  steady  thrust  energy  produced  by  pneumatic  power 
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was  incorporated  in  one  system  for  shallow  (1  m)  penetration,  and  a  vibra¬ 
tory  energy  system  was  incorporated  in  another  design  to  obtain  deeper 
penetration  (3  m). 


4.4.1  Diver  Operated  Core  Sampler 

Simplicity  of  operation  and  a  minimal  number  of  mechanical-electri¬ 
cal  features  were  other  criteria  used  for  design  of  a  shallow  penetration 
corer.  A  study  of  figures  24  through  27  gives  details  of  the  design  and 
hardware  incorporated  in  the  tool.  Details  covering  the  laboratory  and 
field  test  of  the  diver  operated  corer  are  given  in  a  report  by  Jenkins 
and  Takeyama  (1971).  The  field  testing  and  evaluation  of  the  prototype 
corer  were  done  in  conjunction  with  the  sampling  program  for  the  reflec¬ 
tivity  experiment  conducted  in  San  Francisco  Bay  (Barnes  et  al.,  1972). 


Figure  24.  Diver  operated 
sampler  (retracted  posi¬ 
tion)  . 
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Figure  26.  Sample  retainer  —  modified  Shelby  tube 
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Figure  27.  Diver  operated  sampler  piping  diagram. 


4.4.2  Bottom  Sitting  Vibracorer 

An  impact-vibrator  hammer  and  drill  pipe  mounted  in  a  breakaway 
tripod  (O'Brien  and  Duley,  1971)  was  used  to  sample  sediments  to  depths 
of  3  m  (fig.  28).  These  samples  were  tested  in  the  same  manner  as  the 
diver  operated  core  samples,  but  their  primary  purpose  was  to  provide 
lithologic  information  beyond  the  depth  capability  of  the  diver  core 
sampler.  The  energy  component  of  the  system  was  a  National  Vibrator, 
Model  BH-6.  The  weight  of  the  vibrator  is  350  pounds  with  an  impact 
force  of  240,000-in  pounds.  It  uses  60  cfm  of  air  at  70  psi  and  has  a 
blow  count  of  110  blows  per  min.  The  exhaust  ports  are  fitted  with  check 
valves  to  prevent  water  entering  the  cylinder. 


The  drill  pipe  is  made  from  round  mechanical  tubing  drawn  over  a 
mandrel  with  a  3.50  in  OD,  3.25  in  ID,  and  wall  thickness  of  0.125  in. 

The  drill  bit  is  made  from  regular  schedule  80  pipe,  flush  jointed  to  the 
pipe,  with  an  OD  of  3.50  in,  and  an  ID  of  2.875  in.  The  bit  was  shaped 
with  a  crowd-out  relief  angle  of  15°  and  is  fabricated  to  permit  the  in¬ 
stallation  of  polyethylene  preloaded  cartridge  (fig.  29).  The  cartridge 
is  used  as  a  lining  for  the  drill  pipe  to  reduce  wall  friction,  and  be¬ 
cause  it  is  collapsible  it  also  serves  as  a  core  catcher. 

The  vibrator  and  drill  pipe  are  coupled  together  by  a  Victaulic 
coupling.  Grooved  ends  of  the  coupling  are  welded  to  the  drill  pipe  and 
vibrator. 

To  provide  the  support  needed  for  near  perpendicular  penetration, 
and  to  ensure  stability,  a  tripod  was  designed  which  supports  the  drill 
pipe  in  the  upright  position  during  the  initial  1-2  m  of  penetration. 

Once  penetration  has  proceded  to  this  point,  a  coupling  on  the  drill  pipe 
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strikes  the  top  of  the  tripod,  causing  the  tripod  to  collapse  and  allows 
the  drill  pipe  to  continue  to  penetrate  to  full  length.  Sea  floor  slopes 
greater  than  15°  could  cause  the  drill  to  overbalance  and  fall  on  its  side. 
Complete  details  covering  the  design  and  fabrication  of  the  tripod  are 
covered  by  O'Brien  and  Duley  (1971). 

The  combination  of  components  used  in  this  core  sampler  make  a 
tool  that  is  light  in  weight  and  simple  to  construct  and  operate. 

4.5  Precision  Navigation  and  Positioning  System 

A  precision  radio  location  system  was  used  in  all  the  field  exper¬ 
iments  and  sampling  programs  described  in  this  report.  The  main  system, 
along  with  a  specially  designed  subsystem,  is  used  in  referencing  ship's 
position  to  geodetic  control  and  for  automatic  real-time  tracking  during 
survey  operations.  The  electronic  subsystem  enables  real-time  on-line 
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digital  and  graphical  conversion  of  lane  identification  coordinates  to 
x-y  grid  coordinates.  The  subsystem  was  developed  at  this  laboratory 
(Barnes  and  Newman,  1972).  The  ship's  position  and  track  are  accurately 
located  to  within  ±  3  m  based  on  daily  repeatability  checks  of  the  cali¬ 
brated  value  which  is  fixed  by  a  land  survey  and  represents  a  geodetic 
control  point  at  the  ship's  berth. 

Three  peripheral  input/output  units  have  been  added  to  the  system 
since  the  work  reported  in  Barnes  et  al.  (1972).  These  include  a  digital 
fathometer,  a  line  follower  meter,  and  a  teletype.  The  programmable  cal¬ 
culator/computer  allows  interface  capabilities  of  these  subsystems  with 
the  basic  navigation-positioning  data.  A  block  diagram  of  the  system  is 
shown  in  figure  30. 

4.5.1  Automatic  Line  Follower 

The  line-follower  meter  indicates  the  position  of  the  vessel; 
whether  it  is  left  or  right  of  the  prescribed  track  line  and  is  cali¬ 
brated  to  read  how  far  off  track  the  vessel  is,  from  0  to  100  ft.  The 
incorporation  of  the  line-follower  into  the  navigation  subsystem  allows 
the  boat  operator  to  follow  a  more  precise  course  along  the  track  line, 
and  alleviates  constantly  monitoring  the  x-y  plotter  during  survey  oper¬ 
ations. 


4.5.2  Recording  and  Display 

A  U-Tech  Model  109  digital  depth  converter  was  incorporated  into 
the  subsystem,  along  with  a  teletypewriter.  The  "red”  and  "green"  lane 
counts,  x-y  coordinates,  time  of  print-out,  and  digital  depth  is  printed 
out  by  a  digital  printer  at  2  min  intervals  on  command  from  the  program 
in  the  calculator.  At  the  same  time,  the  teletypewriter  prints  out  the 
time,  digital  depth,  and  x-y  coordinates.  The  format  of  the  teletype¬ 
writer  printout  is  conducive  to  quick  observation  of  time,  depth  and 
position  data.  Greenwich  Mean  Time  was  used  for  all  recorded  events  in 
the  navigation  and  positioning  operation  and  for  all  magnetic  tape  re¬ 
cordings  of  acoustical  data. 
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Figure  SO.  Block  diagram  of  Basic  navigation  system  and  subsystems. 
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The  track  lines  of  the  survey  area  were  pre-plotted  on  transparent 
overlays  referenced  to  conventional  maps  using  the  California  State  Grid 
Coordinate  System  (Lambert  Coordinates).  The  scale  of  the  overlays  was 
1  in  =  1,000  ft.  Each  overlay  sheet  covered  an  area  15,000  ft  by  10,000 
ft,  with  the  size  of  the  sheet  being  limited  by  the  size  of  the  x-y  plot¬ 
ter  (fig.  31).  The  boat  operator  is  guided  by  observing  the  line  meter 
or  th?  x-y  plotter  with  the  appropriate  overlay  installed.  The  vessel 


4.5.3  Indexing 


A  marker  interface  which  simultaneously  marks  the  records  of  the 
fathometer  seismic  recorder  and  Raydisl  strip  chart  recorder,  and  which 
also  controls  the  Raydist  digital  printer,  is  activated  by  a  digital 
clock  contained  in  the  subsystem  interface.  The  digital  clock  gives 
selectable  time  outputs  ranging  from  1  to  30  min.  Records  are  period¬ 
ically  annotated  with  time  and  station  numbers,  referenced  to  the  ves¬ 
sel's  location  in  lane  counts  which  are  printed  out  on  the  Raydist  digi¬ 
tal  printer.  In  addition  to  outomatic  operation  of  the  marker  interface, 
the  marker  can  also  be  manually  operated  to  index  events  such  as  those 
that  occur  during  maneuvering,  marking  of  core  sites  with  buoys,  and 
calibration  points. 

4.5.4  Shore  Stations 

The  two  shore  station  sites  required  for  operating  the  system 
were  located  approximately  15.5  miles  apart,  southeast  of  the  operations 
area.  The  north,  or  "red",  shore  station  was  a  35  ft  whip  antenna  powered 
by  two  thermal -electric  generators  and  was  located  at  a  site  south  of 
Watsonville,  California,  directly  over  a  USC&GS  benchmark  designated  Holm. 
This  site  was  at  an  elevation  of  378  ft,  and  about  6  miles  inland  from  the 
coast,  in  line  of  sight  with  the  calibration  point  at  the  ship's  berth  in 
Santa  Cruz  Yacht  Harbor.  The  south  or  "green"  station  was  also  a  35  ft 
whip  antenna  but  was  powered  by  110  V  ac.  This  station  was  located  on 
the  Fort  Ord  military  base  over  a  USC&GS  benchmark  designated  Reserve, 
at  an  elevation  of  485  ft,  and  approximately  2.6  miles  inland  from  the 
coast,  also  in  line  of  sight  with  the  calibration  point. 

The  red  and  "green"  lane  count  —  each  lane  equals  148.2282  ft  — 
from  each  of  the  respective  shore  stations  to  the  vessel  is  automatically 
connected  to  x-y  grid  coordinates  on  a  Hewlett-Packard  Model  7100A  pro¬ 
grammable  calculator.  These  coordinates  are  then  fed  to  a  Hewlett-Pack¬ 
ard  9125A  digital  plotter  to  provide  a  real  time,  automatic  plotting, 
navigation  and  positioning  system. 


Mis-ties  of  an  average  of  0.058  lanes  (8.59  ft)  from  the  "red" 
station,  and  0.087  lanes  (12.89  ft)  from  the  "green"  station  were  noted 
at  the  calibrate  point.  Since  the  Raydist  is  a  non  line-of-sight  system 
using  ground  wave  transmissions,  it  was  felt  these  mis-ties  must  be  at¬ 
tributed  to  one  or  more  of  3  causes,  namely,  (1)  changes  in  the  velocity 
of  wave  propagation,  (2)  insufficient  warm-up  time  before  commencing  op¬ 
erations,  and  (3)  characteristics  of  the  propagation  path  from  both  shore 
stations,  taking  into  consideration  the  overland  distance,  soil  conduct¬ 
ivity,  surrounding  vegetatijn,  and  topography. 

5.  FIELD  OPERATIONS 
5.1  Test  Sites 

5.1.1  Monterey  Bay 

Following  the  1971  San  Francisco  Bay  Study  (Barnes  et  al . ,  1972), 
three  areas  in  Monterey  Bay,  California,  were  selected  for  continuing  in¬ 
vestigations  (fig.  32).  A  program  was  laid  out  to  conduct  several  sur¬ 
veys  including  sampling,  subbottom  profiling,  bathymetry,  photo  reconnais¬ 
sance,  and  tasks  associated  with  the  acoustical  and  electrical  resistivity 
research  studies.  Only  one  area  (Area  A)  was  surveyed  due  to  time  restric¬ 
tions  and  limitation  of  funds.  The  underwater  photo  reconnaissance  was 
not  attempted  because  of  poor  visibility.  Actual  production  on  the  area 
selected  coincides  well  with  the  program  intended  (fig.  33).  Sampling 
with  the  diver  operated  corer  was  limited  to  water  depths  less  than  140 
ft,  thus  reducing  the  number  of  core  sites  occupied.  Possible  hazard 
of  decompression  sickness  was  the  main  reason  for  not  extending  core  op¬ 
erations  into  deeper  water. 

The  bottom  sediments  in  the  survey  area  are  of  several  types, 
consisting  of  gravels,  medium  to  fine  sands,  silts  and  muds.  Hard  lime¬ 
stone  reefs  extend  seaward  from  Point  Santa  Cruz  and  Point  Soquel ,  and 
a  generally  northeast-southwest  trending  reef  outlier  parallels  the  east 
edge  of  the  Point  Santa  Cruz  extensions.  The  bottom  topography  is  gen¬ 
erally  smooth  with  a  gentle  southerly  gradient  except  in  the  reef  areas 
(Tig-  34).  The  wide  variability  of  sediment/rock  type,  shallow  water  and 
even  topography  were  factors  influencing  the  selection  of  the  area. 
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ACOUSTIC  COVERAGE-CORE  SITE  INDEX  MAP 
Santa  Cruz  Harbor  —  Monterey  Boy 
Santa  Cruz,  California 


Figure  33.  Acoustic  and  electrical  resistivity  coverage  —  core  site 

index  map. 
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BATHYMETRIC  MAP 
Sonia  Cruz  Morbor  -  Monterey  Ba» 
Sonia  Cruz,  California 

i 


Figure  34.  Bathymetric  map3  Santa  Crus  Harbor  —  Monterey  Bay 


The  program  in  Area  A  comprises  30  square  miles  and  was  designed 
to  apply  fine  grained  survey  practices.  A  seismic  grid  pattern  of  600  ft 
line  spacing  was  laid  out  for  the  continuous  reflection  coefficient  map¬ 
ping  survey  and  towed  electrical  resistivity  survey.  Core  locations  were 
designated  In  two  densities,  the  shallow  water  area  grid  was  laid  out 
with  hole  spacings  of  1,200  ft  by  1,200  ft,  and  the  deep  water  area  re¬ 
flects  hole  spacings  of  3,000  ft  by  3,000  ft  (fig.  34).  Two  samplers 
were  used,  but  the  data  presented  in  this  study  was  mostly  confined  to 
correlations  based  on  analyses  from  the  1  m  core  samples  taken  using  the 
diver  operated  corer. 


5.1.2  Other  Test  Sites 

Field  tests  for  the  shear  wave  and  electrical  resistivity  experi¬ 
ments  were  conducted  at  two  other  sites  besides  Monterey  Bay.  The  first 
test  area  was  located  in  lake  sediments  near  Fort  Cronkite,  San  Francisco, 
California.  This  body  of  water,  which  was  once  open  to  the  ocean,  is  now 
separated  by  a  beach  of  gravelly  coarse  sand.  Shear  wave  experiments  were 
carried  out  on  the  seaward  side  of  the  lake  in  homogeneous  beach  material. 
The  test  site  offered  the  advantages  of  an  unlayered  sediment  and  no  hind- 
erance  from  water  wave  action  and  tides.  Electrical  resistivity  tests 
were  conducted  at  another  test  site  located  in  a  small  bay  off  the  Gulf 
of  California,  Mexico  (section  5.4.1). 

Experiments  employing  both  the  shear  wave  generator  and  electrical 
resistivity  apparatus  were  conducted  in  Monterey  Bay,  California.  Shear- 
wave  measurements  and  core  samples  were  obtained  at  ten  stations  (refer 
to  f’g.  33  for  location  of  stations  94,  103,  118,  121,  122,  123,  X-l,  X-2, 
X-3,  and  mooring).  These  stations  were  selected  on  the  basis  of  providing 
a  wide  variability  of  sediment  type.  Although  the  method  of  measuring 
shear-wave  was  essentially  the  came  at  the  Fort  Cronkite  and  Monterey  Bay 
test  sites,  the  field  procedures  were  different  and  will  be  treated  sepa¬ 
rately.  Figure  33  also  shows  locations  and  track  lines  where  the  direct 
current  resistivity  experiments  were  conducted. 
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5.2  Reflection  Coefficient  Mapping  Experiment 

The  reflectivity  experiment  was  performed  aboard  the  R/V  Doodlebug 
(fig.  35)  with  a  three  man  crew,  consisting  of  a  navigator,  an  electron¬ 
ics  technician,  and  a  boat  operator.  The  navigator  monitored  the  preci¬ 
sion  navigation  system  and  attended  to  operation  of  subsystems  for  accu¬ 
rate  tracking  and  position.  He  also  annotated  the  various  location  and 
depth  recording  printouts,  as  well  as  preplotting  and  selecting  the  lines 
:o  be  traveled  as  indicated  on  the  x-y  plotter  and  by  the  line  follower 
(fig.  36). 

The  electronics  technician  operates,  calibrates,  and  monitors  the 
subbottom  and  reflectivity  measuring  and  recording  apparatus  and  annotates 
the  subbottom  profile  records  and  magnetic  recording  tape  logs  (fig.  37). 

The  first  operation  of  the  day  was  to  calibrate  and  set  all  the 
navigation  devices  while  tied  up  at  the  berth,  which  was  also  the  navi¬ 
gation  calibration  point.  The  tape  recorder,  signal  generator  and  ac 
voltmeter  also  were  turned  on  and  allowed  tc  warm  up  before  commencing 
operations.  The  boat  then  proceeded  to  the  selected  track  lines  and/or 
core  sites  for  the  day.  Buoys  marking  the  core  sites  were  numbered  in 
reference  to  the  site  number,  and  track  lines  were  accordingly  annotated 
(fig.  33).  Buoys  left  overnight  at  a  prescribed  core  site  would  be  re¬ 
checked  to  insure  that  they  were  still  in  position  on  the  site. 


Figure  35.  R/V  Doodlebug. 
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The  seismic  and  reflectivity  operation  consisted  of  loading  the 
tape  recorder  with  a  pre-numbered  magnetic  tape  roll.  The  channels  were 
calibrated  and  the  signal  generator  set  to  either  2.5  kHz  or  12  kHz  and 
adjusted  to  a  1  V  rms  output.  The  two  pinger  monitoring  channels  on  the 
recorder  were  calibrated.  The  12  kHz  signal  was  observed  for  60  to  90 
sec  or  more  while  the  tape  output  level  was  monitored  by  the  ac  voltmeter 
and  oscilloscope.  Each  tape  channel  gain  control  was  adjusted  for  a  1  V 
rms  indication  on  the  voltmeter.  The  same  procedure  was  used  to  cali¬ 
brate  the  sparker  monitoring  channels,  except  that  a  2.5  kHz  signal  was 
used.  The  calibration  times  are  noted  in  a  tape  log,  which  was  kept  con¬ 
stant  iy  for  purposes  of  noting  regular  gain  changes  and  malfunctions. 

Shortly  before  the  line  to  be  profiled  was  reached,  the  boat  was 
stopped  and  the  pinger  transducer  lowered  into  the  water  and  pinned  into 
position.  The  sparker  sled  was  lowered  from  its  secured  position  on  the 
stern  with  the  deck  crane  and  allowed  to  drift  to  a  position  about  50  ft 
astern  of  the  boat,  where  it  is  secured  for  towing. 

The  streamer  hydrophone  was  removed  from  its  protective  tube  and 
streamed  behind  the  boat  until  it  was  about  20  ft  astern.  An  outrigger 
holding  the  cable  was  then  extended  so  the  hydrophone  was  approximately 
14  ft  from  the  starboard  side  of  the  boat.  This  aft  and  side  tow  posi¬ 
tioned  the  hydrophone  outside  of  the  turbulence  caused  by  the  boat's  wake, 
thus  reducing  the  noise  level. 

While  the  boat  headed  to  the  start  of  the  survey  line,  the  preci¬ 
sion  recorder,  power  supplies,  and  capacitor  bank  are  energized.  The 
boat  operator  then  set  the  throttle  to  a  prescribed  rpm  for  attaining 
constant  boat  speed  (3-4  knots). 

The  subbottom  profiling  record  was  annotated  in  selected  time  mea¬ 
surements  of  2-5  min  and  the  record  scanned  periodically  to  determine  that, 
all  was  in  proper  working  order.  The  sound  source  impulse  and  receiver 
waveforms  were  monitored  on  the  oscilloscope  so  that  all  the  amplifier 
gains  could  be  adjusted  and  logged  before  the  tape  recording  commences. 
Upon  reaching  the  starting  point  of  the  desired  track  line,  the  tape  re¬ 
corder  and  paper  tape  punch  is  started  and  the  time  Vs  noted  in  the  tape 
roll  log,  as  well  as  on  the  voice  channels  of  the  tape  recorder. 
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During  the  profiling  operation  the  sparker  and  pinger  transducer 
receiver-signals  were  constantly  displayed  on  the  dual  channel  Tektronix 
Type  422  oscilloscope,  except  for  occasional  monitoring  of  the  transmitted 
impulse  waveforms.  The  monitoring  of  all  these  signals  was  done  to  in¬ 
sure  that  the  amplitudes  were  of  the  proper  level  for  taping  and  that  no 
distortion  was  present  because  of  noise  or  a  malfunction  of  any  system. 

The  presence  of  all  signals,  including  the  IRIG  B-Time  Code  and  voice, 
could  be  observed  on  the  recording  level  meters  on  each  recording  ampli¬ 
fier  of  the  tape  recorder. 

The  tape  roll  logs  and  verbal  comments  on  the  voice  channel  in¬ 
cluded  the  tape  roll  number,  data,  line  traveled,  direction  of  travel, 
calibration  times,  any  change  in  recording  amplitudes,  signal  source 
changes,  ar.d  the  start  or  stopping  of  the  tape  recorder  or  sources.  Also 
noted  were  sea  conditions,  and  occasionally  noting  sparker  impulse  wave¬ 
forms  or  any  other  events  the  technician  thought  may  influence  the  wave¬ 
forms. 

5.3  Shear  Wave  Experiment 
5.3.1  Field  Procedure  -  Lake  Test  Site 

The  electronics,  powered  by  a  portable  generator,  were  stationed 
on  the  beach  and  the  probe  placed  in  sediments  that  were  under  about  2  ft 
of  water.  Both  the  source  probe  and  receivers  were  implanted  by  hand  to 
a  depth  of  approximately  0.43  m  (some  experiments  were  conducted  with  the 
receivers  at  a  shallower  depth).  The  four  geophone  units  forming  the 
linear  receiver  array  were  inserted  at  predetermined  distances  from  the 
source  through  the  use  of  a  1  ft  scaled  rod.  Various  geophone  unit 
spacings  were  employed,  but  the  end  geophone  unit  of  the  array  was  never 
further  than  5  ft  from  the  source.  The  experiments  conducted  consisted 
of  recording  the  signals  from  the  four  horizontal  geophone  outputs,  the 
four  vertical  geophone  outputs,  or  two  horizontal  and  two  vertical  out¬ 
puts.  Each  experiment  required  2  to  4  source  shots  to  obtain  the  correct 
amplitude,  delay  time,  and  time  scale  settings.  After  acquiring  data  for 
a  particular  experiment,  the  source  and  receivers  were  relocated  to  a  near¬ 
ly  station  of  the  same  sediment  type,  and  a  new  test  was  then  conducted. 

70 


mill  I.  ,  I  i  ,»-:  LJ'UJ^PW .  J  ^ ?  ~ 


water,  an  additional  system  was  required  to  implant  the  source  probe.  The 
system  employed  was  attached  to  the  diver  operated  sampler  (fig.  38).  The 
system  basically  includes  3  air  operated  pistons  that  displace  an  inner 
frame  in  an  up-down  vertical  motion.  The  source,  attached  to  the  frame 
moves  up  or  down  a  track  mounted  on  an  outer  frame.  When  on  site,  the 
sampler  with  SR  generator  is  lowered  to  the  bottom  with  the  inner  frame 
in  its  highest  position  (i.e.,  the  pistons  are  retracted).  On  the  bottom, 
the  divers  actuate  valves  located  on  the  outer  frame.  This  high  pressure 
air  released  from  a  scuba  tank  forces  the  pistons  to  move,  resulting  in 
the  inner  frame  and  the  attached  source  with  probe  to  be  forced  downward 
into  the  sediment.  To  retract  the  probes,  the  valves  are  reversed  caus¬ 
ing  the  pistons  to  raise  the  inner  frame  and  source  probe  from  the  sedi¬ 
ment.  Mo  difficulty  was  experienced  in  employing  this  source  insertion 
method,  even  in  gravelly  coarse  sand. 

Before  Ibwering  the  entire  system,  the  four  geophone  units  were 
correctly  aligned  and  then  clamped  into  four  circular  rings  on  a  split 
bar.  With  the  corer  resting  on  the  bottom  and  after  inserting  the  source 
probe,  a  diver  removed  the  split  bar  from  the  corer  assembly.  Two  holes 
on  one  end  of  the  split  bar  were  positioned  over  two  corresponding  verti¬ 
cal  rods  located  on  the  corer  baseplate  near  the  source  (fig.  38b).  The 
divers  pushed  the  split  bar  down  over  the  rods  with  the  simultaneous  pen¬ 
etration  of  the  stakes  on  the  geoprione  units  into  the  sediment.  The 
divers  then  detached  the  split  bar  leaving  the  geophone  units  implanted. 
This  procedure  provided  a  linear  array  of  evenly  spaced  and  correctly  or¬ 
iented  receivers  (1  ft  spacings  with  the  end  unit  4  ft  from  source).  In 
some  cases  the  divers  implanted  the  geophones  without  the  split  bar. 

After  inserting  the  geophone  units,  the  divers  surfaced,  and 
acoustical  measurements  were  taken.  The  first  two  stations  went  well, 
but  a  current  drain  into  the  geophone  outputs  developed  at  the  third 
station.  This  caused  the  traces  on  the  oscilloscope  to  wander  extens¬ 
ively  when  recording.  The  current  drain  was  traced  to  the  30  V  battery 
triggering  the  air  gun.  The  problem  was  partially  rectified,  although 
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.  Inner  and  Outer  frame  f  .  Scuba  Tank 

.  Shear  Wave  Source  G.  Split  Bar  with  Circular  Rings 

.  Source  Probe  H.  Two  Vertical  Rods 


there  remained  some  wandering  traces.  Several  source  shots  were  required 
to  position  and  focus  in  on  the  desired  wave  arrivals.  An  average  of  24 
shots  per  station  were  needed  to  obtain  the  necessary  information  on  the 
direct  SH-wave  >nd  direct  p-wave;  therefore ,  sediment  disturbance  could 
be  extensive,  even  though  the  water  stream  plus  air  from  the  source  probe 
,  we>e  directed  away  from  the  propagation  path  of  the  acoustical  waves. 

Once  the  data  acquisition  was  completed,  the  divers  submerged,  re¬ 
tracted  the  source,  and  removed  the  geophone  units.  Core  samples  were 
•  then  obtained  at  the  site.  The  complete  cycle  of  positioning  the  ship, 

inserting  the  probes,  acquiring  acoustical  data,  removing  the  probes,  and 
core  sampling  took  between  3  and  4  hours. 


5.4  Direct  Current  Resistivity  Experiment 
5.4.1  in  Situ  Array 

The  electrical  resistivity  system  was  tested  in  two  locations:  a 
small  bay  off  the  Gulf  of  California  at  Puerto  Escondido,  Baja  California, 
Mexico;  and  Monterey  Bay,  California.  The  tests  at  Puerto  Escondido  were 
conducted  in  warm,  clear  water  over  a  carbonate  sand  bottom;  those  at 
Senta  Cruz  were  in  cold,  murky  water  over  a  bottom  of  silty  sand.  The 
in  situ  array  was  tested  at  both  sites. 

In  actual  operation,  a  line  marked  for  Wenner  electrode  spacings 
of  10,  30,  100,  and  300  cm,  is  laid  on  the  bottom,  with  the  ends  anchored 
by  diving  weights.  The  current  and  voltage  electrodes  then  are  first  in¬ 
serted  into  the  sediment  at  the  locations  for  a  =  10  cm  and  the  value  of 
the  voltage;  A V,  recorded  for  several  va-ues  of  the  current,  I.  The  elec¬ 
trodes  then  are  moved  to  the  next  larger  electrode  spacing  and  the  proce¬ 
dure  repeated.  As  skin  contact  resistance  is  low  in  salt  water,  care  is 
taken  by  the  divers  not  to  touch  both  current  electrodes  at  the  same  time, 
to  avoid  the  possibility  of  a  dangerous  shock. 

Electrode  spacing,  current,  voltage,  and  water  depth  and  tempera¬ 
ture  are  recorded  on  the  slate,  and  a  water  sample  taken  for  later  salin¬ 
ity  measurement.  Although  slow  and  cumbersome  to  use,  especially  in  water 


of  poor  visibility,  this  prototype  system  proved  reliable  in  field  opera¬ 
tion,  and  yielded  data  which  appears  to  confirm  the  feasibility  of  the 
sea-floor  horizontal  array  concept  (fig.  39). 

5.4.2  Surface  Towed  Array 

The  simple  towed  resistivity  system  described  in  section  4.3.2  was 
tested  only  in  Monterey  Bay.  The  purpose  of  the  test  was  to  determine 
whether  data  consistent  with  bottom  properties  which  were  measured  di¬ 
rectly  or  were  determined  by  other  geophysical  methods  could  be  obtained 
by  a  towed  resistivity  system.  A  Wenner  array  was  used,  with  electrode 
spacing  (a)  of  20  m.  While  it  would  have  been  desirable  to  run  several 
array  spacings,  extending  -he  array  spacing  to  greater  distances  proved 
impractical . 

The  array  was  towed  from  the  R/v  Doodlebug  at  about  3.5  knots,  us¬ 
ing  the  Raydist  navigation  system  described  in  section  4.5.  The  track 
covered  is  shown  in  figure  33. 


Figure  39.  Formation  factor  vs  a-spacing:  diver  operated  Wenner  array. 
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5.5  Sampling  Program 


m&mems 


t 


The  objective  of  the  sampling  program  was  to  determine  by  core  an¬ 
alysis,  the  mass  physical  and  engineering  properties  of  the  sediments  at 
the  same  locations  where  acoustic  reflectivity  measurements  were  made. 

The  sampling  program  was  conducted  aboard  the  R/v  roar’s  Avk,  a 
converted  LCM6,  56  ft  LOA  and  14  ft  beam  (fig.  40).  An  8-ton  telescoping 
crane  mounted  on  the  bow  of  the  vessel  was  used  to  raise  and  lower  equip¬ 
ment  over  the  side.  Air  for  the  pneumatic-hydraulic  diver-operated  samp¬ 
ler  and  the  vibratory  sampler  was  supplied  by  a  compressor  delivering  85 
cfm  of  air  at  100  psi . 

Since  the  sampling  ship  did  not  have  an  accurate  positioning  sys¬ 
tem,  the  R/v  Doodlebug  would  drop  marker  buoys  at  each  drill  site  for  the 
projected  day's  work  program.  The  choice  of  the  mooring  used  at  each  lo¬ 
cation  (one  or  two  point  anchoring),  was  dependent  upon  the  water  currents 
wave  action  and  the  velocity  and  direction  of  the  wind.  After  anchoring. 


) 


the  sampler  to  be  used  at  a  particular  site  was  prepared  and  put  over  the 


side  of  the  ship. 

A  vibratory  bottom-sitting  drive  sampler  (open  type)  was  used  on 
several  occasions.  This  unit  was  capable  of  obtaining  a  sample  300  cm 
long  and  7.0  cm  in  diameter.  Sampling  with  this  equipment  was  limited 
to  nine  core  sites  due  to  operating  difficulties  (section  4.4.2,  fig.  28). 


/ 


Figure  40.  R/V  NOAA  ’s  Ark. 
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The  procedure  for  obtaining  a  core  with  the  vibracorer  was  to  lower  the 
vibracorer  and  tripod  slowly  over  the  side  of  the  vessel  until  the  tripod 
reached  the  bottom. 

This  was  indicated  by  a  float  attached  to  the  tripod  with  a  length 
of  line  the  same  depth  as  the  water.  When  the  float  sank  into  the  water 
the  tripod  was  on  the  bottom.  At  this  point,  the  air  was  turned  on  to 
the  impact-vibrator  and  the  winch  operator  played  out  the  wire  rope  at¬ 
tached  lo  the  upper  part  of  the  vibrator.  When  there  was  slack  in  the 
wire  rope  it  was  determined  by  experience  that  the  drill  pipe  had  pene¬ 
trated  to  its  maximum  depth.  The  winch  operator  would  then  raise  the  vi¬ 
brator  and  tripod  slowly  onto  the  deck. 

When  the  diver  operated  sampler  was  used,  divers  were  sent  down  to 
activate  the  valves  for  coring,  remove  and  seal  the  sample  tube,  and  re¬ 
place  another  tube  into  the  sampler.  A  cycle  of  operation  including  an¬ 
choring,  sampling,  retrieving  the  anchor,  and  moving  to  the  next  sample 
site  would  take  approximately  1-2  hrs.  When  wind  and  water  conditions 
were  ravorable  and  malfunctions  or  equipment  breakdown  occurred,  five 
sample  sites  could  be  cored  in  one  working  day.  Depth  of  water  and  the 
number  of  divers  available  also  affected  the  number  of  samples  taken  per 
day.  All  dives  were  made  within  the  U.S.  Navy,  no  decompression  limits. 

At  125  ft  a  sample  could  be  taken  within  a  total,  surface-bottom-surface, 
time  of  6  min. 

Most  of  the  samples  were  obtained  with  the  diver-operated,  bottom¬ 
sitting  drive  sample’’  (piston  type)  described  in  section  4.4.  The  sam¬ 
pler  was  capable  of  obtaining  a  sample  70  cm  long  and  7.9  cm  in  diameter. 
The  length  of  samples  taken  ranged  from  a  maximum  of  65.6  cm  to  a  minimum 
of  4.5  cm.  When  a  sample  was  taken,  the  upper  piston  was  sealed  and  a 
plastic  cap  forced  over  the  bit  end  under  water  -  thus  insuring  in  situ 
water  contents.  The  tube,  with  sample,  was  then  handed  aboard  the  re¬ 
search  vessel  NOAA’s  Ark  and  placed  in  racks  in  a  vertical  position,  or¬ 
iented  with  top  of  sediments  up.  When  the  ship  returned  to  its  mooring 
at  Santa  Cruz,  the  sample  was  taken  ashore  and  shipped  by  truck  to  Moss 
Landing,  California,  where  a  soils  laboratory  had  been  set  up. 
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6.  GEOLOGICAL  DATA  ANALYSIS  AND  RESULTS 
6.1  Laboratory  Procedures  -  Sample  Analysis 

Generally,  there  was  a  4  to  6  hr  interim  between  actual  taking  of 
the  sample  and  the  processing  in  the  laboratory.  The  sample  tube  was 
cleaned,  and  weighed  to  0.1  g.  The  seals  were  then  removed  and  the  ex¬ 
cess  water  was  siphoned  off,  weighed  and  recorded.  Preweighing  of  all 
sample  tubes  facilitated  the  calculation  of  the  total  weight  of  the  sample. 
The  volume  of  the  sample  was  calculated  knowing  the  cross  sectional  area 
of  the  sample  tube  and  the  length  of  sample  retained  in  the  tube.  The 
sample  tube  was  then  clamped  into  a  horizontal  core-ejector  and  extrusion 
of  the  sample  commenced. 

Eighty-two  samples  were  collected  from  the  upper  Monterey  Bay  Area 
with  71  engineering  samples  obtained  by  the  bottom-sitting  drive  sampler, 
piston  type,  and  9  lithologic  samples  with  the  bottom-sitting  vibratory 
sampler,  open  type.  Both  types  of  samplers  were  designed  and  developed 
at  MMTC  (Jenkins  and  Takeyama,  1971;  O'Brien  and  Duley,  1971). 

Table  A-4.1  in  Appendix  A-4,  summarizes  the  equipment  used,  cor¬ 
rected  depth  and  other  information  related  to  the  samples.  Mathematical 
equation  used  to  calculate  sample  volumes  were: 


(21) 


where  Cg  is  the  area  ratio,  Dw  is  the  maximum  outside  diameter  of  the  sam¬ 
pler  cutter,  and  Dg  is  the  inside  diameter  of  the  sample  cutter; 


C 


o 


(22) 


where  CQ  is  the  outside  clearance  ratio  controlling  outside  friction  and 
Dt  is  the  outside  diameter  of  the  sample  tube; 


(23) 
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where  is  the  inside  wall  friction  and  Dg  is  the  inside  diameter  of  the 
sample  tube.  The  gross  recovery  ratio, 


(24) 


where  Lg  is  the  length  of  the  sample  obtained  and  H  is  the  depth  of  pene¬ 
tration  of  the  sampler.  This  factor  was  not  determined  in  the  sampling 
because  the  sampling  equipment  was  not  monitored  to  determine  penetration 
depth.  Table  1  summarizes  the  parameters  and  ratios  for  the  two  types  of 
samplers  used  in  this  project. 


6.1.1  Wet  Unit  Weight,  ywet 

Wet  unit  weight  was  determined  by  direct  measurements 


W  -  w 
_  total  wtaro 


total 


where  WtQta^  is  the  weight  of  the  sample  plus  the  sample  tube  and  all 
seals;  Wtare  is  the  weight  of  the  sample  tube  and  all  seals  alone,  and 
Vtotal  1S  tlle  volume  sediments  in  the  sample  tube,  y ■  t  repre¬ 

sents  the  in-place  unit  weight  or  bulk  density. 

Wet  unit  values  ranged  from  a  maximum  of  2.87  g/cm3  at  core  site 
No.  50,  to  a  minimum  of  1.85  g/cm3  at  core  site  No,  103.  Higher  values 
Ywet  were  noted  in  the  use  of  the  vibratory  drive  sampler  as  compared 
to  the  diver  operated  drive  sampler  at  the  same  location.  This  is  prob¬ 
ably  caused  by  plugging  of  the  sample  tube  and  subsequent  compacting  of 
the  sediment  within  the  tube. 

6.1.2  Specific  Gravity  of  Solids,  Gg 

A  representative  sample  of  approximately  100  g  was  retained  and 
overdried  at  110°C  for  10-12  hrs  or  until  a  constant  weight  was  observed 
This  sample  was  then  allowed  to  cool  in  a  desiccator,  then  it  was  split 
and  weighed.  The  volume  of  the  solids  (V  )  was  determined  directly  by  a 
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Table  1.  Sampling  Equipment  Specifications 


Piston-type 

Drive  Sampler  3.34  3.37  3.44  3.44  5.19  0.0  0.72 

(Diver  operated) 


Open-type 

Drive  Sampler  2.75  3.69  3.75  3.75  86.5  1.7  23.5 

(Vibracorer) 

Beckman  Model  930  Air  Comparison  pycnometer.  The  weight  of  solids  (Ws) 
was  determined  with  a  Sartorius  balance  to  0.01  g. 

Specific  gravity  of  the  solids,  Gs,  was  determined  from 

W_ 

=  y—  .  (26) 

Corrections  for  salt  content  in  the  specific  gravities  were  not  made. 
Specific  gravities  of  the  solids  for  over  100  samples  ranged  from  2.57  to 
2.75.  The  majority  of  samples  were  between  2.62  -  2.70  and  an  average  of 
2.66. 

6.1.3  Water  Content,  to;  Void  Ratio,  e;  Porosity,  <j> 

Water  content  is  a  ratio,  in  percent,  of  the  weight  of  water  in  a 
sediment  sample,  Ww,  to  the  weight  of  the  over-dried  solids,  W$,  and  is 
defined  as 

Ww 

to  =  ^  x  100  .  (27) 

s 

In  this  analysis,  determination  of  water  content  was  found  by  the  calcium 
carbide  gas  pressure  method  and  follow  AASH0  Designation  (see  T217-67  I, 
A2-2,  ASTM  1970,  p.  456-459).  Using  local  soils  from  Monterey  Bay  and 
distilled  water,  a  conversion  curve  was  calculated  by  regression  analysis 
with  a  correlation  coefficient  of  .99  and  standard  error  of  estimate 
equal  to  ±  1.93.  Equation  28  was  used  to  calculate  water  content,  or, 

a)  =  1.33  (calcium  carbide  method)  -  3.34.  (28) 
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Another  measure  for  water-saturated  sediments  is  the  void  ratio,  e, 
which  was  the  ratio  of  the  volume  of  the  void  spaces,  Vy,  to  the  volume  of 
the  solid  sediments,  Vg,  or 


e 


(29) 


Void  ratio  was  calculated  in  the  laboratory  from. 


e 


G  y  V 
s  Yw 


-  1  , 


(30) 


where  v  is  the  unit  weight  of  sea  water. 

*4 

Porosity  is  the  relation  expressed  in  percent  of  the  volume  of  the 
voids,  Vv,  to  the  total  volume,  V,  which  is  defined  by, 

V 

<|>  =  ~  x  100.  (31) 

In  this  report  porosity  was  computed  from  other  measured  parameters 
and  expressed  as 

♦  =  ile  x  100  .  (32) 

6.1.4  Grain  Size  Analysis 

A  representative  sample  from  any  one  lithologic  zone  observed  in 
the  core  sample  was  weighed,  put  into  a  mixing  bowl  of  water,  stirred  and 
allowed  to  stand  for  one  hour.  The  sample  was  wet  sieved  in  a  U.S.  Stan¬ 
dard  No.  200  sieve  and  all  material  retained  was  dried  overnight  at  110°C, 
or  until  a  constant  weight  was  observed.  The  wash  water  was  evaporated 
in  the  drying  oven  at  110°C  and  both  fractions  were  cooled  in  desiccators. 
The  coarse  fraction  was  sieved  using  a  nest  of  sieves  (openings  12.5,  4.69, 
1.981,  1.397,  .991,  .701,  .495  mm)  for  7  min  and  the  fraction  passing  No. 

35  (.495  mm)  was  sieved  in  a  nest  of  sieves  (openings  .351,  ,246,  .175, 
.124,  .088,  .074  mm)  for  10  min  in  a  gyrotating  sifting  machine.  The  size 
fraction  finer  than  No.  200  (.074  mm)  was  split,  and  a  representative 
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sample  of  approximately  35-40  g  was  used  for  the  determination  of  the 
grain  size  distribution  of  the  fine-grained  soil  by  the  hydrometer  analy¬ 
sis  method.  The  sample  was  placed  in  a  beaker  with  250  cc  of  distilled 
water  and  5  cc  of  a  10  percent  by  weight  solution  of  sodium  hexmetaphos- 
phate  (calgon)  and  hand-mixed.  The  suspension  was  transferred  to  a  dis¬ 
persion  cup  (ASTM  D422-63),  aid  mechanically  mixed  for  10  min.  The  dis¬ 
pensed  soil  suspension  was  then  transferred  to  a  1000  cc  graduated  cylinder, 
distilled  water  added  to  increase  the  volume  to  1000  cc,  and  placed  in  a 
constant  temperature  bath.  Approximately  1^-2  hrs  in  the  bath  were  needed 
to  bring  the  solid  suspension  up  to  test  temperature  at  25°C.  All  hydro¬ 
meter  tests  were  run  4  hrs  to  find  the  percentage  break-off  point  of  the 
silt  to  clay  fraction  at  5  microns  (U.S.  Bureau  of  Soils  classification). 

Results  of  sieve  and  hydrometer  analysis  were  graphed  on  a  cumula¬ 
tive  frequency  curve  by  weignt  from  which  the  mean  diameter  was  selected 
from  average  diameter  of  the  16,  50,  and  84  percentiles  (Folk  and  Ward, 
1957). 

6.1.5  Shear  Strength  (Penetrometer) 

A  standard  Proctor  penetrometer  with  interchangeable  needles  was 
first  jsed  by  the  divers  to  obtain  in  situ  measurements  but  was  discarded 
as  beirg  unmanageabte  and  not  enough  bottom  time.  Penetration  tests  were 
conducted  in  the  laboratory  using  a  pocket-size  penetrometer  which  was 
calibrated  in  T/ft2  and  Kg/cm2.  All  measurements  were  taken  as  the  sam¬ 
ple  was  slowly  ejected,  but  still  within  the  core  barrel. 

6.1.6  Shear  Strength  (Vane  Shear-Torvane) 

Attempts  were  made  to  use  a  hand  vane  shear  tester  operated  by 
divers  to  obtain  in  situ  measurements  of  shear  strengths,  but  this  was 
discarded  because  of  low  bottom  visibility  and  consumption  of  bottom  time. 
Vane  shear  measurements  were  made  in  the  laboratory  with  a  hand  Torvane 
device  calibrated  at  1  rev  =  1.0  Kg/cm2.  Here  again  measurements  were 
taken  within  the  sample  tube  as  the  sample  was  ejected.  Most  of  the  mea- 
. i rements  were  taken  when  a  visual  change  in  lithology  was  evident. 
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6.2  Interpretation  of  Sediment  analyses 

TaL'e  A-4.3  in  Appendix  A-4  represents  a  tabulation  In  sediment 
analyses  made  on  approximately  85  core  sites  contained  in  Santa  Cruz  har¬ 
bor,  Monterey  Bay  (fig.  33).  The  datums  for  parameters  of  mass  physical 
properties,  such  as,  porosity,  density,  mean  grain  diameter  and  sorting 
coefficient,  as  well  as  engineering  properties  of  shear  strength,  were 
posted  at  respective  core  sites,  and  contoured  interpretations  prepared. 
Figures  41  through  46  show  the  results  of  the  interpretations.  Summary 
of  the  parameters  datums  at  each  core  site  and  the  acoustic  values  are 
contained  in  table  A-4. 4  in  Appendix  A-4. 

6.2.1  Interrelationship  of  Sediment  Properties 

Computer  analyses  of  the  parameters  representing  the  mass  physical 
and  engineering  properties  were  done  using  a  regression  analysis  program. 
Cross-correlations  of  the  parameters  resulted  in  the  correlation  coeffi¬ 
cients  given  in  table  2.  The  degree  of  correlation  is  reflected  in  per¬ 
cent  with  a  value  of  1.000  representing  a  perfect  correlation.  The  rela¬ 
tive  interdependence  or  interrelationship  of  one  parameter  to  another  is 
reflected  in  the  correlation  coefficient.  As  predicted  in  other  studies 
cited  previously  (section  2.1),  the  correlation  of  porosity  and  density 
(.952)  is  very  good.  Also,  the  correlation  of  sorting  coefficient  to 
mean  grain  diameter  is  within  reasonable  expectations  (.830).  Evidence 
of  correlation  between  penetrometer  and  vane  shear  measurements  (.317) 
is  poor,  and  may  be  attributed  to  the  predominance  of  sand,  which  has 
always  been  an  inherent  problem  in  obtaining  reliable  shear  strengths. 

Other  correlations  appear  to  be  negligible,  thus  their  significance  in 
diagnosing  the  interdependence  or  interrelationship  of  these  parameters 
appears  to  be  minor. 

6.2.2  Sediment  Dynamics  and  Deposition  Characteristics 

The  correlation  cf  acoustical  properties  to  sediment  properties  is 
given  in  section  7.1.3.  The  value  of  relating  the  two  measurements  in  a 
positive  correlation  is  obvious  considering  the  economics  and  time  involved 
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Figure  42.  Density  distribution  —  Santa  Cruz  harbor 


Grain  size  (mean  diameter)  distribution  -  Santa  Cruz  harboi 
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Santa  Cruz  Harbor  -  Monterey  Bay 
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Figure  44.  Sorting  coefficient  distribution  —  Santa  Cruz  harbor 
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SHEAR  STRENGTH-  PENETROMETER  (kg /cm2) 
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Figure  46.  Shear  strength  ( penetrometer )  distribution  —  Santa  Cruz  harbor. 
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Table  2.  Correlation  Matrix  —  Sediment  Data 
(Linear  Degression  Analysis) 


Mass  Physical  and  engineering  Properties 

Correlation 

Coefficient 

Standard 
Error  of 
Estimate 

Porosi ty  vs : 

Dens i ty 

0.952 

±0.075 

Shear  Strength  (Penetrometer) 

0.075 

Shear  Strength  (Torvane) 

0.110 

Grain  Size  (mean  dia) 

0.022 

Sorting  Coefficient 

0.026 

Bulk  Density  (wet)  vs : 

Shear  Strength  (Penetrometer) 

0.117 

Shear  Strength  (Torvane) 

0.091 

Grain  Size  (mean  dia) 

0.069 

Sorting  Coefficient 

0.030 

Shear  Strength  (Penetrometer)  vs: 

Shear  Strength  (Torvane) 

0.317 

±0. 197 

Grain  Size  (mean  dia) 

0.062 

Sorting  Coefficient 

0.005 

Shear  Strength  (Torvant)  vs: 

Grain  Size  (mean  dia) 

0.086 

Sorting  Coefficient 

0.098 

Grain  Size  (mean  diameter)  vs: 

Sorting  Coefficient 

0.830 

±0.301 

in  sampling  by  core  drill  from  a  surface  platform  at  sea.  The  Monterey 
Bay  study  would  be  incomplete  if  the  sediment  dynamics  and  depositional 
characteristics  were  not  mentioned.  Figure  47  is  &  map  showing  the  dis¬ 
tribution  of  sediment  types  in  the  test  area.  As  can  be  seen,  consider¬ 
able  variation  occurs.  Figure  48  represents  3-component  diagrams  of  the 
distribution  of  sediment  types.  The  distribution  was  quite  different 
from  that  found  in  the  San  Francisco  Bay  study  (Barnes  et  al.,  1972,  p. 
125),  in  which  a  marked  propensity  toward  finer  grained  size  fractions 
was  noted.  Obviously,  the  forces  affecting  sediment  transport  and  depo¬ 
sition  are  different  in  the  two  areas.  Wolfe  (1970),  described  the  ef¬ 
fects  of  the  coastal  current  and  swell  induced  forces  on  the  distribution 
of  sediments  in  Monterey  Bay.  The  results  of  his  study  are  shown  by  vec¬ 
tor  diagrams  that  give  the  direction  of  major  sediment  transport  in  the 
northern  radius.  The  dominant  wave  refraction  patterns  have  also  been 
added  to  provide  a  complete  concept  of  the  dynamics,  as  shown  in  fiqure 
49.  The  deflection  of  flow  by  the  promontory  at  Santa  Cruz  headland 
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MAJOR  SEDIMENT  TYPES 
(AFTER  VCHUPI  a  GAAL-1963  ) 


Sonfo  Cruz  Horbor  —  Monterey  Boy 
Santa  Cruz,  California 


Figure  47.  Bottom  sediment  composition  —  Santa  Cruz  harbor 
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Figure  48.  Three -component  diagrams  sediment  types  -  Santa  Cruz  harbor. 
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<b) 


Figure  48.  Continued. 
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Hgure  49.  Sediment  transport  pattern  and  wave  refraction  pattern 
Monterey  Bay  (after  Wolfe,  1970). 
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coupled  with  swell -induced  bottom  currents  over  submarine  topographic 
obstructions  (rock  outcrops)  have  resulted  in  an  erratic  distribution  of 
many  sediment  types  in  the  near  shore  regions  of  Santa  Cruz  harbor.  Very 
poor  sorting  of  sediment  grain  size  combine  with  a  high  variability  of 
sediment  type.  This  has  resulted  In  many  of  the  problems  of  sand  deposi¬ 
tion  in  the  vicinity  of  the  Yacht  Harbor  basin  entrance  (fig.  33),  as  well 
as  the  absence  of  deposition  further  south  in  the  beach  areas  of  the  city 
of  Capitola. 

The  area  within  the  30  ft  contour  on  the  bathymetric  map  (fig.  34) 
and  inward  to  the  shore,  could  be  defined  as  a  secondary  circulation  cell. 
It  is  a  type  of  occurrence  more  commonly  found  in  estuaries  and  in  closed 
bays  than  the  offshore  areas  of  the  coast.  It  is  only  because  at  the 
unique  set  of  topographic  circumstances,  both  the  marine  environment  and 
the  shore  line  that  cause  such  a  situation  in  Santa  Cruz  harbor.  Except 
for  the  storm  season,  very  little  material  is  deposited  from  the  rivers 
that  enter  the  harbor.  Nearly  all  the  sediment  load  Is  borne  In  by  the 
longshore  current  and  swell -induced  bottom  currents.  The  forces  found  in 
the  so-called  secondary  cell  existing  in  the  shallow  water  region,  are 
created  out  of  components  of  the  dominant  forces.  It  is  the  activity  of 
this  subordinate  current  within  the  secondary  cell  that  plays  the  impor¬ 
tant  role  in  deposition  of  the  sediment  load  or  ultimate  transportation 
to  other  localized  areas  along  the  harbor  coastline. 

7.  GEOPHYSICAL  DATA  ANALYSIS  AND  RESULTS 
7.1  Continuous  Reflection  Coefficient  Mapping  Experiment 

7.1.1  Data  Reduction  Program 

Program  tasks  under  a  subcontract  to  the  Raytheon  Company,  Ports¬ 
mouth,  Rhode  Island,  specified  that  automatic  data  processing  and  reduc¬ 
tion  be  performed  on  a  series  of  analog  acoustic  signals  acquired  in  the 
reflection  coefficient  mapping  experiment.  Original  specifications  of 
data  volume  consisted  of  222  track  miles  of  acoustic  traverse  recorded  on 
7-channel  magnetic  tape  (fig.  32  -  Area  A).  Processing  of  this  data  con¬ 
sisted  of  analyzing  sequential  sets  of  15  acoustic  pulses,  spaced  along  a 


track  line  at  600  ft  intervals  for  each  of  4  data  channels  (2  channels  of 
0-5  kHz  signals  and  2  channels  of  12  kHz  signals).  Each  4-channel  set  in¬ 
cluded  therefore,  60  discrete  pulses  resulting  in  2200  sets  entailing  ex¬ 
amination  of  132,000  individual  pulses  (see  figs.  50-54). 

Statistical  analysis  also  was  included  on  each  separate  set  of  data 
values  for  the  0-5  kHz  and  the  12  kHz  signals  with  supporting  data  involv¬ 
ing  primarily  calculation  of  water  depth,  recorder  gain  variations  (acqui¬ 
sition  and  digital  processing  gains)  and  time  code  indices. 

The  acoustic  data  was  further  processed  to  yield  printout  listings 
and  digital  tapes  delineating  the  following  parameters  for  each  data  set: 

1.  Beginning  and  ending  time  code  indices. 

2.  Water  depth  -  15  individual  values  with  the  average 
range  and  standard  deviation. 

3.  Computer  spherical  spreading  loss  used  to  compute 
bottom  loss  values  with  the  average,  range  and  stand¬ 
ard  deviation. 

4.  Sound  speed  in  water  column  used  for  depth  and  spheri¬ 
cal  spreading  loss  computations. 
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Figure  50.  Over-all  program  flow  -  computer  processing  routine. 
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Figure  51.  Detail  program  next  time  code 

flow:  initialization.  _ „ 

SAMPLE  DATA 
(FIG.  3) 

PROCESS  DATA 
1  (FIG.  A. 5) 

REPEAT  FOR  SET 
OF  TIMES 


COMPLETE 


5.  Acquisition  system  gain. 

6.  Digitizing  system  gain, 

7.  Amplitude  of  individual  pulses  with  the  average, 
range  and  standard  deviation. 

8.  Bottom  loss  or  reflection  coefficient  with  the 
average,  range  and  standard  deviation. 

Initially,  processing  entailed  locating  in  the  seismic  wave  train, 
the  peak  pressure  pulse  transmitted  by  the  acoustic  source  and  the  first 
seismic  reflection  event  received  from  the  seafloor.  The  events  sampled 
were  located  by  determining  the  time  interval  from  transmission  to  recep¬ 
tion  of  the  water  path  (ray  path  travel  time  based  on  analysis  of  the  12 
kHz  wave  train)  from  time  domain  analysis. 

Regarding  the  quantities,  volume,  values  to  be  found,  and  the  pres¬ 
entation  format,  it  was  mutually  agreed  with  the  contractee  that  due  to 
the  experimental  nature  of  the  project  some  flexibility  was  required.  As 
originally  conceived,  the  analog  data  acquisition  phase  was  designed  to 
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Figure  52.  Detail  program  flow: 
sample  data  and  save  on  tape. 


REPEAT  FOR  BLOCK1 
OF  PINGS  <  ' 


NOT  I TY  OPERATOR 
OF  NEXT  TIME 
TO  BE  PR0CESSF.0 


PROCESS  DATA 
(FIGS.  4,5) 


acquire  simultaneous  wave  trains  (0-5  kHz  signals  at  one  per  second  and 
12  kHz  signals  at  4  per  second)  with  ships  speed  adjusted  in  such  a  manner 
that  a  nominal  600  ft  track  line  interval  would  equal  a  1  min  processing 
cycle.  Accordingly,  machine  processing  was  set  to  first  treat  0-5  kHz 
signals  and  then  the  12  kHz  signal  ,  generating  15  pulse-sets  for  each, 
every  minute.  The  actual  data  acquisition  program  resulted  in  both  the 
0-5  kHz  and  12  kHz  signals  transmitted  simultaneously  on  a  one  per  second 
basis  (section  4.1.3).  A  slower  ships  speed  also  resulted  in  the  600  ft 
track  line  interval  being  equivalent  to  a  2  min  processing  cyc?e,  thus 
nearly  trebling  the  number  of  data  points  acquired. 

Further,  because  of  the  high  variability  encountered  in  the  acous¬ 
tic  wave  trains  and  unforeseen  problems  which  resulted  from  variaticn  in 
the  data  acquisition  system  operation,  a  preliminary  editing  cycle  on  all 
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Figure  S3.  Detail  program  flow:  a 

process  data. 


data  was  required.  This  preliminary  editing  cycle  was  necessary  to  in¬ 
sure  that  all  time  intervals  specified  could  in  fact  be  computer  pro¬ 
cessed.  Original  program  scheduling  was  predicted  upon  this  preliminary 
editing  cycle,  having  been  performed  prior  to  initiation  of  the  process¬ 
ing,  but  this  turned  out  to  be  impractical  during  field  operations.  The 
impact  of  these  data  acquisition  changes  resulted  in  higher  cost  per  pro¬ 
cessed  unit.  It  was  neceisary,  because  of  budget  limitations,  to  select¬ 
ively  reduce  the  volume  of  data  to  be  processed.  Only  40  percent  of  the 
total  data  volume  was  processed  as  a  result. 

An  example  of  the  processed  data  format  is  presented  in  table  3. 

As  the  volume  of  processed  data  in  this  format  is  extremely  extensive,  a 
data  summary  was  prepared  and  is  included  in  Appendix  A-5. 
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FOR  USE  IN 
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Figure  54.  Detail  program  flow;  maxima  location. 


Digital  Processing  Objective.  Digital  processing  objectives  in 
volved  determination  for  each  set  of  pulse  pairs  (i.e.,  i  =  1,2  ---15, 
see  table  3)  of: 

(a)  V  .  =  peak  value  of  Pinger  Transmit, 

Ua  I 

(b)  V  .  =  peak  value  of  Pinger  Receive, 

(c)  t.  =  time  difference  between  Pinger, 

transmit  and  Pinger  Receive 

(d)  V  •  e  peak  value  of  Sparker  Transmit, 

b  A  I 

(e)  V  .  =  peak  value  of  Sparker  Receive, 
and  calculate 

(f)  D.  =  water  depth. 


lyVl.SSHVpr^Vpxc),  and 


(2)  Sparker  reflection  coefficient,  where  Gp  and  Gs 
are  defined  as  respective  system  acquisition  and 
digitizing  gains: 


Gs(4D^832)’*(Vsn./Vsx1). 


In  addition  calculate: 

(3)  Pinger  Bottom  Loss: 


-20  Log10  Rpi- ,  and 


(4)  Sparker  Bottom  Loss: 


-20  Log10  Rs.. 


Having  completed  these  determinations,  statistically  treat  each  of  the  re¬ 
sulting  nine  arrays  as  follows: 

Defining  each  array  as  A,  (j  array), 

J 

Aj  j  =  1,2,... 9 

iL 

and  each  element  within  these  arrays  as  a-,  (itn  element),  calculate 

*  J 

,  15 

Rj  =  X5-  £  aij  (average),  (37) 

Maximum  and  Minimum  (range),  (38) 

a^i  =(lT  ^  (standard  deviation).  (39) 

j  i  —  i  ^  J  J  a 
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For  these  analyses,  a  processing  cycle  time  of  1  min  was  used 
where  the  first  pulse  pair  (i=l)  occurred  7.0  sec  before  the  listed  time 
code.  The  results  of  these  determinations  and  calculations  were  output 
in  the  format  presented  in  table  3. 

Digital  Processing  Techniques.  As  outlined  above,  there  were  5 
basic  wave  form  parameters  to  be  determined  for  each  of  the  15  pulse  pairs 
within  each  pulse  set  (V  ^ ,  V  ^ ,  Vsxj,  Vsr . ,  t^).  Having  determined 
these  parameters,  a  series  of  subsequent  data  reduction  steps  were  per¬ 
formed  to  obtain  the  remaining  arrays  of  required  output  information  (i.e., 

°i*  D'  00'  Rp1’  v  Rs.---etO 

The  general  program  structure  was  constructed  to  process  a  basic 
block  of  15  pulse  pairs  per  set;  each  pulse  pair  was  comprised  of  a  trans¬ 
mit  pulse  and  a  received  pulse  from  cwo  distinct  acoustic  systems.  For 
each  of  the  two  acoustic  systems,  (a  0-5  kHz  system  and  a  12  kHz  system) 
we  ascertained  the  maximum  value  of  the  transmitted  wave  forms  envelope, 
the  maximum  value  of  the  received  waveforms  envelope,  and,  In  the  case  of 
the  12  kHz  system,  the  time  difference  between  these  transmit  and  receive 
max 'mums.  Having  determined  these  values,  the  reflection  coefficients 
and  oottom  losses  as  defined  in  equations  33  through  36,  were  computed. 

Data  Format.  The  original  analog  signal  format  suggested  two 
acoustic  systems  transmitting  at  a  repetition  rate  of  1:4  with  the  pinger 
system  having  the  higher  repetition  rate.  Employing  this  format,  the  max¬ 
imum  available  analog-to-digital  conversion  rate  of  30  kHz  could  be  apr 
plied  to  each  of  the  four  waveforms  of  interest  on  a  non-interfering 
basis.  Specifically,  the  bandwidth  of  the  0-5  kHz;  is  nominally  5  kHz; 
the  12  kHz  approximately  2  kHz.  A  30  kHz  sampling  ra+e  was  thus  suffi¬ 
ciently  above  the  Nyquist  rate  to  insure  individual  maxima  to  be  deter¬ 
mined  with  little  or  no  error. 

The  suggested  analog  signal  acquisition  format  was  not  achieved 
as  system  design  constraints  necessitated  both  systems  transmitting  syn¬ 
chronously  at  one  per  sec  repetition  rate.  As  a  consequence,  the  30  kHz 
sampling  rate  had  to  be  alternated  between  the  Pinger  and  Sparker  channels, 
which  resulted  in  an  equivalent  15  kHz  sampling  rate  per  channel.  For  the 
12  kHz  system  that  had  a  2  kHz  bandwidth,  the  15  kHz  sampling  rate  still 
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allowed  the  maximum  to  be  determined  directly.  For  the  5  kHz  system,  how¬ 
ever,  it  was  necessary  to  use  an  FFT  interpolation  technique  to  obtain  the 
required  degree  of  time  resolution. 

As  a  further  complication,  the  signal  of  the  analog  data  was  dif¬ 
ficult  to  decipher  because  of  either  transient  noise  or  multiple  reflec¬ 
tions.  The  simple  routine  of  finding  the  maximum  within  a  given  observa¬ 
tion  window  could  not  be  routinely  applied.  To  insure  data  quality,  a 
modified  receive-search  procedure  was  developed  in  which  the  received 
data  was  scanned  in  time  until  a  preset  threshold  value  was  exceeded, 
indicating  the  leading  edge  of  the  returning  pulse.  The  procedure  then 
selected  the  largest  value  occurring  in  the  next  t  sec  (t  approximately 
equal  to  the  transmit  pulse  deviation),  storing  this  value  as  the  re¬ 
ceived  maximum. 

Further,  provisions  were  made  to  add  a  small  sample  of  the  pinger 
return  during  the  time  it  should  have  negligible  amplitude  (i.e.,  during 
the  "off'1  period).  The  maximum  of  these  samples,  essentially  system  back¬ 
ground  noise,  also  is  presented.  Large  values  suggest  periods  of  spuri¬ 
ous  data. 

Processing  Parameters.  To  maintain  processing  flexibility,  at¬ 
tempts  were  made  to  make  all  key  parameters  detailed  selectable  at  run¬ 
time. 

Sampling  Data  —  The  sample  rate,  sample  size,  and  onset 
threshold  level  are  selectable. 

Check  Data  —  The  check  sample  size  and  the  channel  ex¬ 
amined  (of  the  four)  can  be  specified. 

Block  Size  —  The  number  of  pings  examined  can  be  varied. 

Receive  Parameters  -  The  onset  threshold  (within  the  data 
acquisition  window)  and  the  subsequent  time  to  find 
a  maximum  are  selectable.  This  permitted  observation 
of  the  waveforms  on  the  scope  and  established  a  thres¬ 
hold  which  lies  above  the  noise  level  yet  assures  de¬ 
tection  of  the  return. 

Anticipated  Depth  —  The  receive  data  observation  window  must 
be  wide  enough  to  accomodate  any  depth  variations  from 
one  block  of  pulse  sets  to  the  next.  This  window  is 
automatically  positioned  for  the  next  block  of  pulse 
sets  by  using  the  previously  calculated  depth.  Its 
location  is  held  constant  for  a  given  block  of  pulses. 


For  the  initial  block,  th;  program  accepts  the  oper¬ 
ator's  estimate  of  this  initial  depth,  requiring  for 
input  the  observed  return  time  as  seen  on  the  scope. 

Miscellaneous  Constants  -  Also  accessible  at  run-time  are 
several  other  constants.  Included  are:  the  speed  of 
sound;  a  pinger  gain  factor  related  to  attenuator  set 
tings  during  record,  and  a  scale  constant  for  each  re- 
ce i ve  channel . 

The  following  parameter  values  (one  exception  discussed  below)  were 


30  kHz 
256  samples 
15  samples 
15 

1  minute 
10  samples 
4800  ft/sec. 


used  during  processing: 

sample  rate  = 

sample  size  = 

check  sample  size  = 

number  of  pings  = 

observe  block  every 
receive  peak  sample  size  = 
speed  of  sound  = 

Occasionally,  shallow  water  required  a  smaller  sample  size  (128  in¬ 
stead  of  256).  This  was  necessary  since  the  transmit  observe  windows  would 
be  of  length  1/30,000  (sample  rate)  by  2  (alternate  sampling)  by  256  (num¬ 
ber  of  samples)  =  17.7  m  sec  or  a  depth  of  41.5  ft.  Thus,  one  could  not 
begin  observing  the  receive  channels  until  after  the  primary  returns  had 
arrived.  By  decreasing  the  sample  size  to  128,  depths  as  small  as  20.8 
ft  cculd  be  processed.  Under  this  condition,  however,  the  observation 
window  is  only  ±  10.4  ft  (as  opposed  to  ±  20.8  ft)  so  that  the  depth  var¬ 
iations  had  to  be  closely  monitored. 

Detail  Program  Flow  Diagrams.  Figures  40  through  44  present  a  de¬ 
tailed  breakdown  of  the  overall  program  flow  presented  in  figure  40. 


7.1.2  Summary  of  Processing  —  Comments  and  Findings 

The  tables  4  and  5  outline  in  summary  form  first  the  requested 
times  to  be  processed  along  with  the  actual  times  processed,  and  second, 
a  copy  of  the  log  made  out  at  the  time  of  processing.  As  a  general  com¬ 
ment,  the  actual  time  required  to  process  any  given  block  of  required 
times  varied  from  1:1  to  2:1  with  a  general  average  of  1.5:1  necessitated 
by  the  required  preprocessing  editing. 
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Table  4.  Processing  Specifications 


Hell  No. 

'  ' 

Requested  l  imes 

Total 

Actual  Times 

Processed 

Iota  1 

Comments 

001 

175553-190353 

09 

175553-190353 

09 

002 

10  1753-173553 

49 

Possibly  10  mins 

processable 

00,'i 

173953-180953 

31 

173953-17 1453 

0 

183153-190353 

33 

175253-180953 

18 

183253-184153 

10 

184453-190353 

20 

00  1 

153353-102953 

55 

153353-102953 

r>5 

000 

104553-173953 

55 

105753-173953 

43 

175153-182353 

33 

175153-132353 

33 

017 

143753-153053 

00 

143753  -153053 

00 

018 

102553-104353 

19 

102953-104353 

15 

174153-175753 

17 

174153-175753 

17 

019 

180953-181953 

11 

180953-181953 

11 

020 

151753-100153 

45 

151853-155353 

30 

100753-104553 

39 

100753-104553 

39 

021 

152353-102553 

03 

152353-102553 

03 

022 

105653-181953 

84 

105053-191953 

84 

029 

190353-195753 

55 

190353-195753 

55 

024 

102553-170553 

41 

102553-170553 

41 

Lontf  Form  only 

170953-174753 

39 

170953-174753 

39 

025 

180053-181353 

14 

1H0053 -181353 

14 

100953-105153 

43 

100953-105153 

43 

020 

172353-184453 

92 

172353-184453 

82 

027 

184753-190753 

21 

185053-190753 

12 

155753-104153 

45 

155853-104153 

44 

028 

105753-174753 

51 

Inverted  Tape 

180753-182153 

15 

Bad  Tape 

029 

182553-193353 

09 

183353-193353 

01 

020 

103153-174153 

71 

103153-174153 

71 

031 

174553-181753 

33 

175153-181753 

27 

192153-201753 

57 

192153-201753 

r>7 
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Roll  No. 

Requested  Times 

032 

171753-182153 

033 

182453-184553 

185153-200053 

034 

1G5553-170953 

035 

171153-180353 

185153-193153 

03(5 

193553-195953 

174253-175953 

037 

192553-201353 

03  3 

201753-202753 

203553-204153 

204953-214353 

039 

171453-180353 

041 

1G0753-165953 

042 

170953-180353 

04G 

154453-162353 

047 

1G3553-170353 

170953-180653 

048 

181953-185153 

185953-191753 

192353-195753 

049 

154553-164153 

1G4G53-165553 

050 

170753-180753 

181553-183553 

051 

184753-193853 

194753-201553 

052 

202553-211153 

054 

055 

182753-183953 

185153-190553 

190753-194053 

194353-201753 

rotal 

Actual  Times 

Processed 

Total 

Comments 

G5 

171753-182153 

65 

22 

183153-184553 

15 

70 

185153-200053 

15 

105553-170953 

53 

172153-180353 

41 

185153-193153 

25 

194053-195953 

20 

18 

174253-175953 

18 

49 

192553-201353 

49 

11 

202353-202753 

5 

7 

203553-204153 

7 

55 

205053-214353 

54 

50 

171453-180353 

50 

53 

0 

Sparker  channel 

55 

0 

Not  processable 

40 

0 

II 

29 

0 

II 

58 

0 

II 

33 

19 

0 

II 

31 

57 

0 

H 

10 

61 

0 

II 

21 

52 

o 

II 

19 

47 

0 

II 

13 

0 

II 

_ 1  n _ 1 _ _ 

055 


15 

34 


Had  Sparker 


Table  4.  Processing  Specifications  (Con't) 


Boll  No. 

Requested  Times 

Total 

Actual  Times 

Processed 

Total 

Comments 

058 

182153-185353 

33 

182153-183553 

m 

185753-133553 

39 

185753-193553 

183753-185353 

17 

Long  Form  only 

059 

194753-202153 

.35 

194753-202153 

35 

two 

170953-173753 

29 

Bad  Sparker 

Data  Summary j  Short  Form.  A  summary  form  of  all  the  processed  data 
generated  per  table  3  is  contained  in  the  Appendix  A-5.  The  summary  sheets 
present  the  numeric  values  of  the  average  reflection  coefficient  and  one 
standard  deviation  for  both  Pinger  and  Sparker  signal  systems,  together 
with  the  average  bottom  loss  values  as  a  function  of  data  and  time  code 
index. 

7.1.3  Correlation  of  Reflectivity  to  Sediment  Properties 

Regression  analyses  were  done  by  computer  processing  to  determine 
the  correlation  coefficients  for  the  12  kHz  and  0-5  kHz  acoustical  data 
to  the  sediment  data.  The  Rayleigh  Reflection  Coefficients  and  bottom 
loss  values  taken  over  areas  circumscribed  about  each  core  site  were, 
respectively,  averaged  for  a  representative  value  that  was  applied  to 
each  of  the  core  site  locations.  The  core  sites  each  were  thus  referenced 
by  respective  acoustical  and  sediment  property  values  or  datums.  These 
data  were  then  cross-correlated  and  correlation  coefficients  obtained 
through  the  regression  analyses.  The  objective  was  to  relate  the  degree 


Table  5.  Shipboard  Data  Analysis  and  Commentary  Log 


Roll  No. 

Findings 

Comments 

001 

Cal's  ok.  Good  data.  Few  drop-outs.  Added 
-10  dll  to  pinger  channels 

Good  run  with  very  few  few  drop-outs. 

003 

Calibration  signals  on  during  time  required 
for  processing.  Added  -10  dB  to  Pinger 
channels.  Reasonable  data. 

Segmented  Runs:  a/.  173953  174433 
Cal's  1745  -  1751 
b/.  175253  -  180953 
c/.  183253  -  184153  Sparker  bad 
1842,43 

d/.  184453  190353 

004 

Cal's  ok.  Deep  water,  multiple  drop-outs. 
Pinger  transmit  low.  Added  -10  dB  to 
Pinger  channels 

Processed  complete  run  drop-outs 
could  not  be  picked  up  due  to  poor 

S/N. 

000 

Calibration  signals  on  during  time  requested 
for  processing.  Not  possible  under  present 
program  configuration.  Added  -10  dB  to 
Pinger  channels.  Large  interference  spi- 
ker  in  Pinger  Rec.  Sparker  Rec. drop-outs. 

Segmented  Runs:  a/.  165753  173953 
Cal's  1646  1656  Inter.  1710  1712 
b/.  175153  182353  Many  data  drop¬ 
outs. 

017 

Cal's  ok.  Reasonable  data.  Few  drop-outs. 

10  dB  in  Pinger  channels. 

Processed  complete  run 

018 

Cal's  on  during  requested  times.  Time-code 
drop-outs.  Multiple  Pinger  Rec.  drop-outs. 
Added  -4.0  dB  to  Pinger  channels 

Segmented  Runs:  a/.  162953  164353 
Cal's  on  to  1629 

b/.  174153  175753  Multiple  drop-outs 

019 

Cal's  ok.  Reasonable  data  added  -4.0  dB  to 
Pinger  channels. 

180953  181953 

020 

Cal's  ok.  Data  starts  150805  Time-Code  skip 
Sparker  Transmit  funny  155853.  End  on 
155353  Problems  after  with  Sparker.  Addec 
-4.0  dB  to  Pinger  channels 

151953  155353  Reasonable  data.  Will 
review  problems  at  later  date. 

020 

155433-160153  Sparker  very  erratic.  -10  dB 
in  Pinger  channel 

Unable  to  auto  process 

020 

160753-164653  Sparker  problems.  Time 

Code  drop-outs.  Very  difficult  to  process 
-10  dB  in  Pinger  channel 

Segmented  Runs:  a/.  160753  161253  ok 
b/.  Missed  161353-162953  Sparker 
Transmit  and  Time  Code  problems 
c/.  163053-164653  ok  Sparker  acting 
funny 

021 

Cal's  ok.  Time  Code  out  161253,  13.  Funnies 
on  tape  -10  dB  Pinger  channels. 

Processed  complete  run 

Table  5.  Shipboard  Data  Analysis  and  Commentary  Log  (Con't) 


Roll  No.  Findings 

Comments 

022 

Cal's  ok.  Finger  transmit  loss  Sparker 
low.  Added  -10  dB  Finger  ehannels. 

Processed  complete  run 

022 

Cal's  ok.  Reasonable  data.  Few  drop-outs. 
-10  dll  Finger  ehannels. 

Processed  complete  run 

02  1 

Cal's  ok.  Reasonable  data.  Few  drop-outs. 
-10  dB  in  Finger  ehannels. 

Processed  complete  run 

025 

Cal's  checked.  Finger  Xmlt  Low  (3.0  Vp) 
Spavker  Xmit  high  (5. 0  Vp)  800j  source. 

-10  dB  in  Pinger  channels.  Sparker  data 
erratic  and  multiple  drop-outs,  lime 

Code  drop-outs. 

Segmented  Runs:  a/.  180053  181353 
high  sparker,  low  Pinger 
b/.  100953  -  105153.  Levels  ok,  but 
returns  erratic.  Questionable  sparker 
data  nm  a/. 

020 

Cal's  checked.  Both  Pinger  and  Sparker 
Xmit's  low  at  beginning  of  run.  -10  dB  in 
Pinger  ehannels.  Pinger  Xmit  about  nor¬ 
mal  at  end  of  run. 

Processed  complete  run 

027 

Cal's  cheeked.  Cal's  on  during  time  required 
to  be  processed.  -10  dB  Pinger  ehannels. 
Highly  variable  sparker  data  Time  Code 
skips.  Good  levels. 

Segmented  Runs:  a/.  185653  -  190753 
Cal's  on  in  beginning 
b/.  155853  164153  Tape  on  155803 
variable  data,  missed  162053,  and 
162553 

028 

Tape  from  net  re-wound.  Rewound  tape. 
Findings;  channels  1.  Noise 

2.  Time  Code 

3.  Noise 

5.  Pulsed  data,  but  not 
of  the  normal  form. 

Tape  to  be  re-examined  later. 

029 

Cal's  eheeked  ok.  On  until  183300  Pinger 
Xmit  low.  -10  dB  in  Pinger  ehannels. 

Time  Code  Skips  190453. 

Segmented  Runs:  a/.  183353  190353 
b/.  190553  193353 

030 

Cal's  ok.  Pinger  low  -10  dB  Finger  ehannels 

Processed  complete  run 

031 

Cal's  eheeked  but  Pinger  still  low.  -10  dB 
in  Pinger  ehannels.  Cal's  on  during  data. 
175153  variable  returns 

Segmented  Runs:  a/.  175153  -  181753 
b/.  192153  -  201753 

032 

Cal's  cheeked,  Pinger  low.  -10  dB  in  Pinger 
ehannels.  Very  erratic  sparker. 

Segmented  Runs:  a/.  171753  -  175253 
b/.  175353  -  182153 

033 

Cal's  eheeked.  Pinger  low.  Cal's  on  during 
data  183153.  Very  poor  sparker  data? 
Highly  variable.  Pulse  Rep.  rate  erratic. 
Time  Code  skips.  Difficult  tape. 

Segmented  Runs:  a/.  183153  -  184553 
b/.  185153  -  200053  Difficult  Tape! 
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Table  5.  Shipboard  Data  Analysis  and  Commentary  Log  (Con't) 


Roll  No*.  Findings 

Comments 

034 

Cal's  checked,  Ptnger  low,  Pulse  Rep.  er¬ 
ratic.  -10  dB  in  Pinger  channels. 

Completed 

035 

Cal's  checked.  Levels  wavered  in  time  Ch  #1 
and  # 5.  Rep.  rate  erratic.  Cal's  on  during 
times  requested.  -10  dB  in  Pinger. 

Completed.  Segmented  Runs. 

03  (i 

Cal's  checked.  -10  dB  in  Pinger,  Cal's  on 
during  time  requested.  Pinger  drop-outs 

Completed,  Segmented  runs 
a/.  194053  -  195953 
b/.  174253  -  175953 

037 

Cal's  checked.  Erratic  Rep.  Rate  -10  dB 
Pingers 

Processed  complete  run. 

038 

Cal's  checked.  On  during  time  requested. 

-10  dB  Pinger  very  low  Pinger  Xmit. 

Difficult  tape.  Segmented  runs. 

039 

Cal's  checked.  Pinger  low.  Pinger  Xmit 

2. 0  Vp  off  tape.  Big  change  174853.  Time 
Code  skips. 

Difficult  tape.  Completed  Run. 

041 

041  052  Sparker  Xmit  Missing 

054 

No  Sparker  Transmit  Total  Tape 

Pinger  channels  ok.  could  process 
with  additional  software  development. 

055 

No  Sparker  Transmit  Total  Tape 

Same  as  054 

058 

Calibration  ok.  Reasonable  data  same  pinger 
Rec.  drop-outs  =184000.  Pinger  Rec. 
channel  too  high.  Added  -10  dB  to  both 
Transmit  and  Receiver  data  Pinger  Gain 
change  at  183700. 

Segmented  Runs:  a/.  185800-193600 
Wrong  gain  factor  used,  should  have 
been  XI.  Requires  sealing, 
b/.  182153—183553  Time  Zero  of 
processing  000053. 

c/.  183653—185353  Some  Pinger  Rec. 
drop-outs. 

059 

All  channels  present  with  Cal's 

Water  depth  became  less  than  40'. 
Required  Segmented  runs. 

0G0 

No  Sparker  Transmit  Total  Tape 

Same  as  054 

of  correlation  between  acoustical  properties  (bottom  loss  and  reflection 
coefficient)  calculated  from  the  peak  amplitudes  generated  in  the  input- 
response  function  of  the  two  sources  to  parameters  measured  in  the  sedi¬ 
ment  analyses  (i.e.,  porosity,  density,  mean  diameter,  sorting  coefficient, 
bearing  capacity,  and  shear  strength).  No  edit  of  the  data  was  accorded 
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as  to  number  of  samples  taken,  other  than  to  restrict  the  area  of  influ¬ 
ence  about  a  core  site.  This  area  was  defined  by  drawing  circles  about 
each  core  site  to  a  radius  that  would  produce  tangency  to  circles  scribed 
about  adjacent  core  sites.  Since  the  core  hole  spacing  was  uniform  in 
pattern  within  the  high  density  area  and  low  density  area,  the  areas  of 
influence  for  the  statistical  sample  were  uniformly  the  same  within  the 
two  densities.  Corrections  were  made  on  the  acoustic  data  for  obvious 
wild  point  values  because  of  errors  incurred  in  the  analog  to  digital  con¬ 
version.  In  such  cases,  the  reflection  coefficient  and  corresponding  bot¬ 
tom  loss  values  were  recomputed. 

The  correlations  obtained  from  the  first  analysis,  which  incorpo¬ 
rated  the  entire  area  surveyed,  are  reflected  in  the  correlation  matrix 
given  in  table  6.  Plots  of  the  data  represented  in  the  mathematical  cor¬ 
relation  are  shown  in  figures  55  through  66.  One  can  immediately  see  the 

Table  6.  Correlation  Matrix  —  Acoustic  and  Sediment  Properties 


VARIABLE 

1  2 

3 

4 

fj 

6 

7 

6 

9 

10 

NUMBER 

1 

1.000  -.963 

.111 

-.085 

.120 

.014 

-.006 

-.009 

-.133 

-.131 

2 

1.000 

-.171 

.039 

-.041 

.066 

.030 

.004 

.OBO 

.0B7 

3 

1.000 

.266 

-.114 

-.074 

-.294 

-.OBB 

.138 

.145 

4 

1.000 

-.083 

-.047 

-.051 

.142 

.046 

.067 

5 

1.000 

-.217 

-.195 

.093 

.097 

6 

1.000 

.141 

.129 

-.27B 

-.262 

7 

1.000 

•  B41 

-.657 

-.625 

8 

1.000 

-.642 

-.620 

9 

1.000 

.989 

3.0  1.000 


BMW12R  -  STEPWISE  REGRESSION  -  VERSION  OF  APRIL  12,  1965 
HEALTH  SCIENCES  COMPUTING  FACILITY,  UCLA 


PROBLEM  CODE  ARPA 

NUMBER  OF  CASES  42 

NUMBER  OF  ORIGINAL  VARIABLES  10 

NUMBER  OF  VARIABLES  ADDED  0 

TOTAL  NUMBER  OF  VARIABLES  10 

NUMBER  OF  SUB-PROBLEMS  24 


VARIABLE  MEAN  STANDARD  DEVIATION 


PROSTY 

1 

44.03357 

4.80310 

DENSTY 

2 

2.08262 

.19319 

PNITMR 

3 

.27602 

.20130 

T-VANE 

4 

.09076 

.10524 

M-DIA 

5 

.20355 

.32998 

SORT  C 

6 

1.45900 

.32749 

SKR  BL 

7 

5.42912 

.73977 

SKR  WO 

8 

5.67864 

.85959 

PGRRFL 

9 

1.13074 

.22659 

PGR  WO 

10 

1.11633 

.23224 

results  show  very  poor  correlations  exist.  Another  representation  of  the 
data  *s  given  in  the  contoured  maps  of  bottom  loss  versus  sediment  proper¬ 
ties  in  figures  67  and  68.  Comparison  of  the  surface  distribution  of 
acoustical  properties  derived  from  the  two  sources  to  any  one  of  the  mass 
physical  properties  and/or  engineering  properties  again  gives  little  in¬ 
dication  of  positive  correlations.  Greater  resolution  of  the  acoustical 
distribution  over  the  entire  survey  area  is  shown  in  figures  69  and  70 
for  selective  comparison  to  figures  41-46. 

The  results  of  the  correlation  represent  a  marked  contrast  to  the 
previous  results  obtained  in  the  San  Francisco  Bay  study  (Barnes  et  al . , 
1972)  and  other  studies  (Breslau,  1964,  1967).  Correlation  coefficients 
derived  from  other  studies  of  acoustics  using  a  pinger  acoustic  source 
(12  kHz)  and  the  broad  band  sparker  (0-5  kHz)  source  are  shown  in  table  7. 
The  difference  in  the  results  of  previous  studies  and  the  Monterey  Bay 
study  may  be  attributed  to  several  causes,  each  of  which  bears  more  data 
processing  of  existing  data  than  could  be  accomplished  in  the  time  span 
of  this  study.  Some  of  these  probable  causes  can  be  enumerated  as;  (1) 
too  large  an  area  was  taken  for  the  representative  statistical  sample; 

(2)  sediment  velocity  variability  in  the  near  shore  region  may  be  more 
severe  than  assumed;  (3)  the  changing  water  depth  coupled  with  topography 
changes  over  short  distances  may  have  altered  normal  incidence  path  be¬ 
tween  source  and  reflector;  (4)  stability  of  the  source-receiver  systems 
in  a  surface  towed  configuration. 

An  attempt  was  made  to  delve  into  each  of  these  probable  causes. 

The  best  explanation  of  the  analysis  can  be  made  by  drawing  an  analogy 
between  the  San  Francisco  Bay  study  and  this  study.  The  first  of  which 
can  be  drawn  on  the  basis  of  size  of  area  surveyed.  The  previous  study 
encompassed  a  relatively  small  area  containing  approximately  1.5  sq  miles 
(1  mile  by  1.5  miles),  while  this  study  encompassed  nearly  30  sq  miles 
(5  miles  by  6  miles).  Although  the  control  for  navigation-positioning 
and  line  density  and  core  site  patterns  was  constant  in  both  stations, 
variability  of  sediment  type  was  much  greater  in  Monterey  Bay  as  well  as 
the  variability  of  water  depth  and  the  range  of  water  depth.  The  previous 
study  contained  sediments  that  consisted  mainly  of  fine  grained  silt  and 
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Porosity  (%) 


Figure  56.  Plot  of  bottom  loss  (12  kHz)  vs  porosity. 


Figure  57.  Plot  of  bottom  loss  (0-5  kHz)  vs  density. 
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Penetrometer  (Kg/cm2) 


Shear  (Kg/cm2) 


Mean  Grain  Diameter  (mm) 
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Figure  62.  Plot  of  bottom  loss  (0-5  kHz)  vs  mean  grain  diameter. 
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66.  Plot  of  bottom  loss  (12  kHz)  vs  sorting  ooeffioient. 


Figure  67.  Acoustics  vs  sediment  physical  propertiet 
Santa  Cruz  harbor. 
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Acoustics  vs  sediment  engineering  properties 
Santa  Cruz  harbor . 


ACOUSTIC  DISTRIBUTION  OF 
BOTTOM  LOSS  { dB ) - SPARKER  SOURCE 
Sonia  Crux  Harbor  -  Monterey  Boy 
Sonia  Crux.  Californio 

\ 


Aooustzo  ctistrrbutzon  of  bottom  loss  (0-5  kHz) 
Santa  Cruz  harbor. 


Fzgvre  69 
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ACOUSTIC  DISTRIBUTION  OF 
BOTTOM  LOSS  !dB  )  -  PINGER  SOURCE 

Sonto  Cruz  Horbor  -  Monterey  Boy 
Sonto  Cruz.  Colifornlo 


Figure  70.  Acoustic  distribution  of  bottom  loss  (12  kHz) 

Santa  Cruz  harbor. 


Table  7.  Comparison  of  Correlative  Values  of 
Barnes  et  al.  (1972)  vs  Breslau  (1967)* 


Correlation 

Standard  Error 

★  ★ 

Coefficient 

dB/Unit 

Barnes  et  al . 

Breslau 

Barnes  et  al . 

Breslau 

Poros i ty 

.73 

.70 

±  .^55 

±.039 

Grain  Size 

•  92 

.646 

±.227 

±.292 

Sort i ng 

.92 

.343 

±.252 

±.660 

Fine  Material 
(Si  is  i  Clay) 

.89 

.786 

±•305 

±.015 

rc 

Bottom  loss  calculations  are  based  on  peak  pressure  amplitude  measurements  and  Rayleigh 
Reflection  Coefficient  determinations. 

The  standard  error  of  estimate  (S)  Is  the  standard  deviation  about  a  line  of  average  re¬ 
lationship.  Assuming  a  normal  distribution,  68  percent  of  all  cases  will  i  le  within 
a  range  of  ±S ,  95  percent  will  fall  within  ±2S,  and  99.7  percent  wiil  fall  within 
±3S.  ir  this  instance,  the  relationship  of  bottom  loss  to  the  other  variabie  is  of 
interest.  The  values  of  S  are  therefore  In  terms  of  decibeis  (dB)  per  unit  of  paired 
variable. 


clay  fractions  with  very  little  material,  the  Santa  Cruz  harbor  sediments 
are  highly  variable  ranging  in  grain  size  from  gravel  to  silt  with  admix¬ 
tures  that  are  poorly  sorted  and  show  great  diversity  in  physical  proper¬ 
ties  and  engineering  properties.  According  to  D.  L.  Bell  (Raytheon  Com¬ 
pany,  Portsmouth,  Rhode  Island,  personal  communication),  it  was  learned 
that  such  high  variability  of  sediment  type  can  have  decided  effects  on 
the  range  of  velocity  of  sound  propagation  through  the  sediments.  Bell's 
recent  studies  of  nearshore  sediments  of  the  continental  shelf  off  New 
England,  show  that  sediment  velocity  changes  greatly  over  short  spacial 
distances  and  that  this  has  an  appreciable  influence  on  corrections  ap¬ 
plied  to  acoustic  attenuation.  Using  advanced  state-of-the  art  tools  and 
techniques  developed  by  Raytheon  Company,  he  has  been  able  to  accurately 
measure  marine  sediment  velocities  to  within  2  percent  of  the  standard 
velocity  for  sea  water.  This  accuracy  has  resulted  in  obtaining  new  in¬ 
formation  on  the  variability  of  marine  sediment  velocities.  According  to 
Bell,  40  percent  of  marine  sediments  in  the  nearshore  environment  range 
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in  velocity  within  5  percent  below  to  2  percent  above  the  standard  velo¬ 
city  of  sea  water.  He  further  believes  that  sediments  can  also  be  cate¬ 
gorized  into  8  or  9  types  through  velocity  analyses. 

Water  depths  in  the  San  Francisco  Bay  study  area  ranged  from  25  to 
40  ft,  as  compared  to  depths  ranging  from  18  to  250  ft  In  Monterey  Bay. 
Tests  for  correlations  as  a  function  of  water  depth  were  made  on  the  re¬ 
sults  obtained  from  the  data  taken  for  this  study.  The  results  are  shown 
in  the  correlation  matrix  given  in  table  8.  Sediment  analyses  and  acous¬ 
tical  data  were  correlated  from  core  sites  along  line  numbers  16  (sites 
38-44);  26  (sites  57-66);  36  (sites  81-88);  43  (sites  100-105;;  and  45 
(sites  109-115)  (figs.  33,  34).  Figures  71-90  show  the  cross-sections 


Table  8.  Correlation  of  Mass  Physical  and  Engineering  Properties 


V8  AC0U8tiC8 

(Linear  Regression  Analysis). 

Mass  Physical  and  Engineering 
Properties 

Correlation 

Coefficient 

Standard 

Error  of 

Estimate  (i) 

Spark er  BL 

Pinqer  BL 

Sparker  BL 

Pinqer  BL 

Line  No.  16 

(sites  38-44) 

Porosity 

0.393 

0.576 

0.3472 

0.7072 

Oensi ty 

0.602 

0.627 

0.3771 

0.7233 

Shear  Strenqth  (Penetrometer) 

0.554 

0.343 

0.3694 

0.6479 

Shear  Strer,  ith  (Torvane) 

0.610 

0.861 

0.3785 

0.3117 

Mean  Grain  L'iareter 

0.582 

0.6217 

0.3740 

0.7217 

Sorting  Coefficient 

0.0618 

0.5789 

0.3225 

0. 7082 

Line  No.  26 

(sites  57-66) 

Porosity 

0.506 

0.275 

0.4814 

0.7601 

Oensi ty 

0.3961 

0.195 

0.5125 

0.8061 

Shear  Strenoth  (Penetrometer) 

0.4998 

0.394 

0.6235 

0.8506 

Shear  Strenqth  (Ttrvane) 

0.304 

0.421 

0.5834 

0.8665 

Mear  Grain  Oiamete 

0.406 

0.8665 

0.5102 

0.674 

Sortinq  Coefficient 

0.271 

0.7032 

0.5935 

0.07882 

Line  No.  36 

(sites  81-88) 

Porosity 

0.485 

0.508 

0.6219 

0.4170 

Oensi ty 

0.570 

0.404 

0.639 

0.4011 

Shear  Strength  (Penetrometer) 

0.577 

0.545 

0.6388 

0.4235 

Shear  Strenqth  (Torvane) 

0.559 

0.524 

0.6339 

0.4198 

Mean  Grain  Oiameter 

0.382 

0.217 

0.5113 

0.3805 

Sortinq  Coefficient 

0.571 

0.425 

0.6370 

0.4040 

Line  No,  43  (sites  100-105) 

Porosity 

0.744 

0.534 

0.3849 

0.5128 

Density 

0.611 

0.489 

0.4561 

0.5035 

Shear  Strenqth  (Penetrometer) 

0.328 

0.779 

0.6064 

0.5733 

Shear  Strength  (Torvane) 

0.925 

0.696 

0.743 

0.696 

Mean  Grain  Diameter 

0.288 

0.673 

0.5518 

0.5450 

Sortinq  Coefficient 

0.7/8 

0.551 

0.7299 

0.5163 

Line  No.  45  ( 

sites  109-115) 

Porosl ty 

0.837 

0.080 

0.6270 

0.9056 

Oensi ty 

0.914 

0.226 

0.4643 

0.9377 

Shear  Strenqth  (Penetrometer) 

0.289 

0.600 

0.6759 

1.1112 

Shear  Strenqth  (Torvane) 

0,797 

0.713 

0.8018 

1.2802 

Mean  Grain  Oiameter 

0.618 

0.618 

1.3458 

1.0751 

Sorting  Coefficient 

0.615 

0.485 

1.3440 

1.0163 
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73.  X-seotion  mean  diameter  and  sorting  Figure  74.  X-seation  penetrometer  and  torvane 
coefficient  vs  acoustics  —  Line  16.  vs  acoustics  —  Line  16. 
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obtained  of  the  sediment  properties  and  acoustical  data,  as  well  as  geo¬ 
logic  cross-sections  along  each  line.  Analyses  show  no  definite  trends, 
except  that  correlation  coefficients  seem  to  run  much  higher  as  the  sta¬ 
tistical  sampling  number  is  reduced  from  the  ove.all  area  surveyed,  to 
line-by-line,  point-to-point  sample  areas.  Geologic  structure,  such  as 
noted  in  the  exposed  bedrock  outcrops  and  thinly  covered  subcrop  are  indi¬ 
cated  by  the  seismic  subbottom  profiles.  These  profiles  and  the  cross- 
sections  show  the  influence  on  the  correlations  obtained  (figs.  17,  41-47, 
67  and  68).  Such  structural  variation  over  short  distances  may  have  di¬ 
rect  bearing  on  the  results  obtained  using  the  high  frequency  (12  kHz) 
source  as  compared  to  the  lower  broad  band  frequency  (0-5  kHz)  source. 

It  was  previously  noted  in  table  6  that  the  correlation  of  the  input-re¬ 
sponse  functions  of  the  two  sources  was  relatively  high  (.657)  despite 
their  inherently  different  power  and  frequency  characteristics.  As  men¬ 
tioned  previously  (section  4.1.1),  the  depth  of  penetration  and  resolu¬ 
tion  of  acoustical  signal  reflected  off  the  bottom  and  subbottom  horizons 
vary  according  to  the  source  employed,  consequently  near  surface  struct¬ 
ural  as  well  as  sediment  property  variations  (within  the  wave  length  of 
the  signal)  can  have  a  decided  influence  on  the  reflection  coefficient 
obtained. 

Stability  of  the  surface-towed  source  and  receiver  systems  was 
not  considered  a  problem.  Obviously,  sea  conditions  are  a  limitation, 
and  should  be  considered  an  important  factor  in  obtaining  reliable  data. 

The  experiments  in  both  the  Monterey  Bay  study  and  the  San  Francisco  Bay 
study  were  conducted  under  the  best  prevailing  we ather  conditions  with 
due  regard  to  the  adverse  circumstances  inherent  in  sea  state,  swell, 
wave  height  and  currents. 

If  the  trend  analyses,  which  shows  that  Higher  correlation  coeffi¬ 
cients  in  shallower  water  exist  as  compared  to  those  noted  in  deeper  water, 
has  any  validity,  it  may  be  concluded  that  a  submerged  tow  vehicle,  which 
would  house  the  source  and  receiver  systems,  could  operate  with  even 
greater  effectiveness.  This  concept  ought  to  be  tested. 


141 


7.2  Shear  Wave  Experiment 

From  the  records,  the  arrival  time  of  the  direct  wave  was  taken  to 
be  the  point  determined  from  the  intersection  of  the  base  line  and  the 
slope  of  the  waveform.  To  obtain  the  shear  wave  velocity  or  compressional 
wave  velocity,  the  time  of  the  first  arrival  from  each  of  the  four  traces 
(produced  by  either  the  four  vertical  geophone  outputs  or  four  horizontal 
geophone  outputs)  were  plotted  on  a  time  distance  graph. 

Before  discussing  the  results,  the  sources  of  error  inherent  in  the 
seismic  method  should  be  mentioned.  One  source  of  error  is  in  the  distance 
between  geophone  units.  Implanting  the  receivers  by  hand  could  result  in 
an  inch  of  misplacement.  For  a  compressional  wave  velocity  of  1500  m/sec, 
an  error  of  approximately  8  percent  in  velocity  can  result  assuming  1  in 
separation  between  two  receivers,  when  in  fact  one  receiver  is  displaced 
an  inch. 

Another  source  of  error  is  in  the  time  scale  on  the  oscilloscope. 

An  appropriate  time  scale  for  a  wave  velocity  of  1500  m/sec  is  .1  m  sec/cm; 
however,  0.5  m  sec/cm  was  the  lowest  time  scale  achieved.  The  possible  er¬ 
ror  that  can  be  introduced  using  this  time  scale  is  9  percent  for  a  1500 
m/sec  wave  velocity.  On  account  of  the  large  errors  introduced,  the  error 
for  a  particular  measurement  is  not  discussed.  For  the  same  reason,  re¬ 
gression  equations  were  not  determined  for  the  acoustical-mass  physical 
property  relationship. 

It  should  be  mentioned  that  the  primary  concern  thus  far  has  been  to 
develop  an  inexpensive  tool  to  Investigate  the  feasibility  of  detecting 
$H-wave  in  saturated  marine  sediments.  Of  secondary  importance  was  con¬ 
struction  of  a  high  precision  monitoring  system.  It  was  felt  that  if  the 
detection  of  shear  waves  was  found  feasible,  a  refined  system  could  then 
be  designed. 

7.2.1  Results  from  Lake  Sediments 

The  determined  acoustical  and  mass  physical  properties  are  listed 
in  table  9.  Figure  91  is  an  example  of  the  signals  received  by  a  geophone 
implanted  2  in  from  the  source  (trace  1,  horizontal  geophone  output;  trace 
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Table  9.  Acoustic  and  Mass  Physical  Properties  at  Beach  Test  Site 


Sediment 

Type 

■  -  -  ■  --  -  -  = 

Mass  Physical  Properti 

es 

Porosity 

m _ 

Grain  Size  (mm)  Gravel 

Median  Mean  (%) 

Sand 

ill 

Clay  and  Silt 

_ ID _ 

gravelly 

46.2 

0.800  0.833  33.7 

66 

.3 

coarse  sand 

Acoustic  Properties  (In  eitu  Measurements) 

- ro 

v 

e 

as 

U) 

(m/sec) 

(m/sec! 

1  (qm/cm3) 

(m*1) 

(Hz) _ 

1690 

125 

1.68 

.921 

604 

Computed  Properties'* 

Q  _1 
vs 

yR 

UI  kr  xr 

0 

E 

0.40 

.02 

.09  4.77  4.75 

.498 

.04 

vi,  K,  X,  and  E  have  units  of  dynes/cm2  X  1010. 


2  vertical  geophone  output)  and  by  another  geophone  unit  implanted  4  in 
from  source  (trace  3,  horizontal  geophone  output;  trace  4,  vertical  geo¬ 
phone  output).  The  receivers  were  inserted  1  in  into  the  sediment,  whereas 
the  exit  port  on  the  source  probe  was  at  a  depth  of  lh  in.  Although  the 
voltage  gain  was  high,  the  signal  quality  was  poor.  This  is  attributed  to 
the  high  attenuation  of  the  shear  wave  and  to  the  low  response  of  the  ver¬ 
tical  geophone  to  the  direct  p-wave. 

In  figure  91,  the  first  arrival  detected  by  each  horizontal  geophone 
is  assumed  to  be  the  direct  S^-wave.  The  polarity  reversal  between  the 
first  arrival  detected  by  the  two  horizontal  geophones  is  caused  by  revers¬ 
ing  the  sensitive  axis  of  the  horizontal  geophone,  whose  output  is  trace  1, 
180°  with  respect  to  the  other  horizontal  geophone.  The  polarity  reversal 

9 

is  a  necessary  condition  of  an  S^-wave  and  was  used  here  to  aide  in  the 
identification  of  the  direct  S^-wave.  Presence  of  the  polarity  reversal 
eliminates  the  first  arrival  as  being  a  compressional  wave,  a  vertically 
polarized  shear  wave,  or  a  Stone! ey  wave;  however,  a  Love  wave  will  also 
produce  a  polarity  reversal  with  the  employed  geophone  orientation,  since 
the  Love  wave  is  a  low  frequency  surface  wave  of  the  type,  formed  from 
the  superposition  of  internal  reflections  within  a  surface  layer. 

\ 
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Horizontal  Scale  ••  5iHsec 

cm 


Vertical  Scale  «  10  — 

cm 

Figure  91.  Signal  wave  forms,  trace  1  and  3:  signals  from  horizontal  geo¬ 
phones,  2  in  and  4  in  from  source,  trace  2  and  4:  signals  from  vertical 
geophones ,  2  in  and  4  in  from  source  (SH:  direct  Su-wave ;  p  direct  p-wave 
Sy:  direct  Sy-wave) . 

A  layer  over  a  half  space  or  a  velocity  gradient  produced  from  ver¬ 
tical  inhomogeneities  is  a  necessary  condition  for  the  existence  of  the 
Love  wave.  For  a  source  and  receiver  located  in  the  distance  between  the 
source  and  boundary  of  the  half  space,  the  direct  SH  wave  will  be  the 
first  arrival.  This  is  attributed  to  insufficient  time  (or  distance)  for 
the  superposition  of  reflections  to  form  the  Love  wave.  Beyond  the  criti¬ 
cal  uistance  (the  distance  at  which  the  direct  wave  and  head  wave  arrive 
at  the  same  time),  the  Love  wave  will  be  the  first  arrival  due  to  the 
superposition  of  the  head  wave  with  other  waves.  For  the  lake  sediments, 
no  visible  change  in  the  sediment  type  was  observed  within  the  top  1.1  m, 
with  the  probes  implanted  0.45  m,  the  next  interface  below  the  water- 
sediment  interface  is  situated  at  a  minimum  depth  of  0.65  m  beneath  the 
probes.  The  minimum  critical  distance  is  thus  1.3  m;  therefore,  it  is 
reasonable  to  assume  that  the  direct  SH  wave  is  the  first  arrival  on  the 
traces  of  horizontal  hydrophones  shown  in  figure  92. 

i/i/i 


Horizontal  Scale*  5 “ crrT 

m  v 

Vertical  Scale  *  ,0^rrT 

Figure  92.  Signals  from  four  horizontal  geophones 
(2'j  3’,  4't  S’  in  from  source). 


In  addition  to  velocity,  the  attenuation  of  the  direct  S^-wave  was 
measured  from  records  shown  in  figure  92.  Four  horizontal  geophones  at  2, 
3,  4,  and  5  in  from  the  source  were  implanted  at  the  same  depth  as  the 
exit  port  on  the  source  probe,  lh  in.  Note  the  polarity  reversal  of  trace 
1  due  to  reversing  the  sensitive  axis  180°  with  respect  to  the  other  three 
geophones.  From  the  measured  attenuation,  the  specific  attenuation  factor 
(Q"1)  was  calculated  to  be  0.40.  The  error  in  using  the  low  loss  approxi¬ 
mations  (i.e.,  Hookean  equations  of  elasticity)  for  this  sediment  is  ap¬ 
proximately  5  percent  for  the  shear  velocity. 

The  arrival  of  the  direct  p-wave  was  difficult  to  measure.  The 
problem  stemmed  primarily  from  the  low  response  of  the  vertical  geophone 
to  the  direct  p-wave.  For  a  vertical  geophone  implanted  further  than  2 
in  from  the  source,  the  actual  beginning  of  the  congressional  signal  was 
absent.  Receivers  within  2  in  of  the  source  detected  the  beginning  of  the 
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compressional  signal  until  an  appropriate  oscilloscope  time  scale  was  set 
to  measure  the  arrival  time  differences.  This  problem  was  rectified  some¬ 
what  by  implanting  the  receivers  and  source  at  different  depths. 


7.2.2  Results  from  Monterey  Bay  Sediments 

The  in  situ  measurements  and  the  calculated  properties  are  listed 
in  table  10.  The  mass-physical  properties  of  the  stations  listed  are  in¬ 
dicated  in  Appendices  A-5  and  A-6.  The  attenuation  of  the  shear  wave  was 
not  measured  because  of  the  previously  mentioned  difficulties  with  the  geo¬ 
phone  probes.  The  source  and  receivers  were  at  different  depths.  As  the 
radiation  pattern  varies  with  the  angle  from  the  source,  attenuation  could 
not  be  measured,  therefore,  the  low-loss  approximation  for  velocity  had  to 
be  used  in  order  to  determine  the  calculated  values  listed  in  table  10.  As 
discussed  earlier,  low-loss  approximation  for  the  shear  velocity  could  re¬ 
sult  in  errors  as  high  as  43  percent  for  a  sediment  with  Q"1  of  2. 

Examples  of  the  better  quality  records  obtained  are  shown  in  figure 
93.  Figure  94  illustrates  a  record  of  poorer  signal  quality.  The  depth 
of  the  source  and  the  receivers  varied.  Generally,  the  source  probe  and 
the  receivers  were  inserted  lh  in  and  1  in,  respectively.  The  four  geo¬ 
phone  units  were  spaced  1  in  apart  with  the  nearest  geophone  unit  1  in  from 
the  source.  Distances  between  the  exit  port  on  the  source  probe  and  the 
receivers  are  noted  on  the  records.  Usually,  only  the  two  nearest  geophone 
units  to  the  source  provided  information  on  the  direct  waves.  The  first 
arrival  at  the  other  two  geophone  units  was  generally  a  head  wave. 


Table  10.  In  Situ  Measured  and  Calculated  Sediment 


Core 

Site 

Vs 

(m/sec) 

VP 

(m/sec) 

* 

(g/cm3) 

Mooring 

567 

1622 

2.06 

XI 

504 

2030 

2.24 

X2 

521 

1512 

2.08 

X3 

— 

1530 

2.23 

94 

438 

1485 

1.99 

103 

308 

1400 

1.85 

118 

524 

— 

2.33 

121 

— 

1440 

2.03 

122 

475 

2.03 

123 

371 

14-"0 

2.09 

dynes/cm2  x  IQ10 


u 

K'  ' 

X 

.66 

4.53 

4.09 

.57 

8.46 

8.09 

.57 

3.97 

3.59 

— 

— 

_ 

.38 

3.37 

3.61 

.18 

.64 

3.39 

3.27 

AC, 

.... 

.... 

.  HO 

.31 

4.14 

3.95 

Properties 


1.91  .431 
1.66  .467 
1.62  .433 

1.11  .452 
.52  .475 


.90  .465 


«* 


Signals  from  four  vertical  geophones. 


Figure  92.  Two  records  obtained  at  the  mooring  site.  Distance  between 
exit  port  on  source  probe  and  receiver  is  provided  on  the  trace. 


a.  Signals  from  four  horizontal  geophones . 
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Figure  94.  Signals  from  four  horizontal  geophones  at  station  122  (dis¬ 
tance  between  exit  port  on  source  probe  and  receiver  is  provided  on 
trace) . 


Difficulty  in  obtaining  the  arrival  of  the  direct  p-wave  also  occur¬ 
red  at  this  test  area.  The  problem  stemmed  from  the  absence  of  a  clearly 
defined  first  arrival  and  from  information  loss  when  setting  the  oscillo¬ 
scope  time  scale  to  an  appropriate  level.  Aim  sec/cm  or  0.5  m  sec/cm 
time  scale  were  used  in  recording  the  direct  p-wave;  therefore,  signifi¬ 
cant  errors  in  arrival  time  differences  could  have  been  introduced. 

With  no  knowledge  of  the  attenuation  of  the  first  arrival  at  the 
horizontal  geophone,  it  is  difficult  to  state  whether  this  first  arrival 
is  the  direct  S^-wave  or  the  Love  wave.  The  feasibility  of  being  the  Love 
wave  was  investigated  for  the  sediment  at  the  mooring  station.  At  this 
station  a  vibracore  was  obtained  which  provided  information  on  the  sediment 
properties  in  the  top  1.7  m.  Extending  to  0.8  m  was  a  very  fine  sand  den¬ 
sity  2.14  g/cm3  overlaying  a  dense  clayey  silt  of  density  2.4  g/cm3.  As 
the  source  probe  was  implanted  0.45  m,  the  minimum  critical  distance  is 
0.70  m.  It  is  therefore  unlikely  the  first  arrival  at  the  nearest  two  hy¬ 
drophones  (0.3  m  and  0.61  m  from  the  source)  is  the  Love  wave.  At  the 
other  stations  a  core  length  of  less  than  0.6  m  was  taken. 


As  arrival  time  measurements  were  generally  made  only  on  the  near¬ 
est  receivers,  1  in  and  2  in  from  the  source,  the  time  for  the  direct 
wave  to  propagate  to  these  receivers  is  very  short.  Therefore,  if  the 
layer  thickness  is  large,  the  direct  wave  can  arrive  before  the  Love 
wave.  If  the  layer  thickness  is  small  or  if  there  exists  a  velocity  grad¬ 
ient  in  the  surface  sediment,  the  Love  wave  can  arrive  first;  but  the  time 
required  to  forri  the  Love  wave  must  be  short.  Consequently,  the  energy 
forming  the  Love  v'ave  cannot  penetrate  further  than  2  in  from  the  source. 
For  this  case,  the  Love  wave  velocity  would  be  comparable  to  the  shear 
wave  velocity.  Therefore,  it  will  be  assumed  that  the  first  arrival  is  the 
direct  S^-wave  (bearing  in  mind  that  further  mathematical  and  experimental 
investigations  are  needed  to  clearly  establish  the  extent  to  which  this 
can  be  assumed). 

7.2.3  Correlation  Between  Acoustical  and  Mass-Physical  Properties 

Previous  investigations  by  Hamilton  (1970b),  Horn  et  al .  (1968), 
Schreiber  (1968),  and  others  have  found  predictable  relationships  for  com- 
pressional  velocity  vs  porosity  and  compressional  velocity.  Porosity  is 
the  ratio  of  the  volume  of  voids  occupied  by  the  water  to  the  total  volume. 
As  the  rigidity  is  of  low  magnitude  in  saturated  sediments,  the  compressi¬ 
bility  of  the  water  largely  determines  the  compressional  wave  velocity 
rather  than  the  compressibility  of  the  grains.  Therefore,  a  strong  cor¬ 
relation  is  expected  to  exist  between  the  compressional  wave  velocity  and 
porosity.  Figure  95  illustrates  the  relationship  obtained  between  the 
compressional  velocity  and  porosity. 

Also,  the  compressional  velocity  is  strongly  related  to  the  sedi¬ 
ment's  grain  size.  This  affect  is  due  to  the  fact  that  density,  porosity, 
and  other  sediment  porperties  are  a  function  of  grain  size.  In  marine 
sediments,  a  decrease  in  grain  size  produces  an  increase  in  porosity,  and 
consequently  a  decrease  in  compressional  velocity.  Velocity-mean  grain 
size  relationship  is  shown  in  figure  96. 

In  discussing  the  relationship  of  the  shear  wave  data  with  the  mass- 
physical  properties,  it  was  more  convenient  to  use  rigidity  as  the  shear 
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Figure  96.  Dynamic  rigidity  vs  mean  grain  size. 


parameter.  Hamilton  (1971)  and  Hardin  and  Richart  (1963)  give  excellent 
discussions  on  the  correlation  of  rigidity  with  other  sediment  properties 
such  as  porosity  and  grain  size.  From  these  reports,  the  following  con¬ 
clusions  were  derived  for  sands: 

(1)  dynamic  rigidity  should  be  a  minimum  in  the  coarse  sands 
as  they  are  rounder  and  possess  fewer  intergrain  contacts; 

(2)  an  increase  in  porosity  will  cause  a  decrease  in  rigidity 
for  sands  of  the  same  grain  size  and  angularity; 

(3)  although  porosity  increases  as  the  grain  sizes  decrease, 
rigidity  may  increase  as  the  grains  become  finer  because 
finer  grains  have  more  intergrain  contacts  and  higher  an¬ 
gularity; 

(4)  at  some  fine  sand-silt  sizes,  the  intergrain  friction  and 
interlocking  will  be  maximum,  and,  therefore,  rigidity  will 
be  maximum  for  natural  uncemented  sediments;  and 

(5)  with  an  increase  in  amounts  of  silt  and  clay,  the  sand 
grains  are  no  longer  In  contact;  the  intergrain  friction 
and  interlocking  are  no  longer  effective,  and  rigidity  be¬ 
comes  dependent  on  the  cohesion  between  finer  particles. 


As  mentioned  previously,  the  determined  empirical  relationships  be¬ 
tween  rigidity  and  porosity  and  between  rigidity  and  mean  grain  size  are 
Illustrated  in  figures  95  and  96.  These  relationships  generally  follow  the 
conclusions  of  Hamilton  (1971)  and  Hardin  and  Richart  (1963)  discussed 
above.  Maximum  rigidity  occurs  for  a  fine  sand  of  grain  size  2.5  phi  unit's 
(or  .24  mm).  For  the  larger  grains  sizes,  the  rigidity  is  less  due  to 
fewer  intergrain  contacts  and  rounder  grains.  Moreover,  the  rigidity  de¬ 
creases  rapidly  when  silt  is  introduced  into  the  sediment.  This  phenomenum 
is  better  illustrated  in  figure  97.  Included  in  figure  98  is  the  variations 
of  rigidity  with  percent  of  sand  and  gravel.  As  shown,  the  increase  in  the 
amount  of  silt  reduces  the  rigidity.  This  is  caused  by  an  increase  in  por¬ 
osity  and  a  decrease  of  density. 


Figure  97.  Dynamic  rigidity  vs  porosity. 
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7.3  Direct  Current  Resistivity  Experiment 


7.3.1  Data  Interpretation 

For  maximum  flexibility  in  the  interpretation  of  field  data,  it 
would  be  desirable  to  have  available  sets  of  interpretation  curves  for  a 
seafloor  array,  similar  to  those  published  for  land  use.  Such  curves  are 
computed  by  calculating  the  potential  field  generated  by  a  single  current 
electrode  placed  over  a  layered  medium,  with  each  layer  having  arbitrary 
resistivity  and  thickness.  A  geometrical  factor,  depending  on  the  array 
configuration  and  spacing,  then  is  used  to  obtain  the  voltage  difference, 
aV,  appearing  between  the  potential  electrodes  for  a  given  current  I13. 

Terekhin  (1962a  and  b)  and  Van'yan  (1956)  have  examined  the  theo¬ 
retical  potential  generated  by  a  current  electrode  on  a  layered  sea  floor, 
but  their  results  are  not  in  a  form  suitable  for  immediate  computation. 
Other  workers  in  the  Soviet  Union  apparently  have  published  interpretation 
curves  for  seafloor  arrays,  but  we  have  not  yet  been  able  to  obtain  copies 
of  their  work.  Equations  In  a  form  suitable  for  computer  solution  are 
given  in  Barnes  et  al .  (1972).  A  simpler,  although  slightly  more  time- 
consuming  method  of  interpreting  sea-floor  resistivity  data  has  been  used 
for  all  of  our  work  to  date. 

For  analytical  purposes,  an  array  on  the  sea  floor  differs  from  one 
on  land  only  in  that  a  layer  of  conductive  sea  water  is  present  above  the 
array.  If  the  depth  and  resistivity  of  the  water  layer  (which  are  easily 
measured)  are  known,  the  apparent  resistivity  pu  of  the  water  layer  may  be 
obtained  from  the  standard  interpretation  curves.  The  water  layer  then  may 
be  considered  as  a  known  resistance  in  parallel  with  the  sea  floor,  and  the 
true  apparent  resistivity  of  the  sea  floor,  pa,  obtained  from  the  total 
measured  resistivity  p^.  by  using  the  formula 


1  _  1  1 


where 


(40) 
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Once  the  value  of  pa  has  been  calculated,  published  standard  interpreta¬ 
tion  curves  may  be  used  to  obtain  the  true  bottom  resistivity  profile. 

The  curves  of  apparent  bottom  resistivity  pa  vs  electrode  separa¬ 
tion  a  obtained  in  Puerto  Escondido  are  shown  in  figure  100.  Locations 
B-5a  and  B-5b  were  about  30  m  apart,  in  water  12  m  deep,  and  location  B-3 
was  in  an  area  about  one  km  away  in  water  6  m  deep.  Bottom  water  salinity 
was  32  0/00  in  both  locations,  and  temperature  ranged  between  21  and  23°C. 

Figure  101  shows  the  geologic  interpretation  obtained  by  matching 
the  curves  of  figure  100  against  the  Wenner  interpretation  curves  published 
by  Mooney  and  Wetzel  ( 1956a ) . '  The  resistivity  profiles  of  the  two  nearby 
locations,  8-5a  and  B-5b,  show  no  evidence  of  bedrock.  The  bottom  sedi¬ 
ment  layers  do  not,  of  course,  extend  to  infinity,  but  the  bedrock  in  this 
area  apparently  is  too  deep  to  be  "seen"  with  a  maximum  value  of  electrode 
spacing  of  300  cm.  The  presence  of  bedrock  is  evident  at  a  depth  of  4.8  m 
at  location  B-3,  which  was  in  shallower  water,  close  to  shore.  It  is  sig¬ 
nificant  that  even  this  low-powered  array,  with  relatively  small  electrode 
spacings,  was  able  to  provide  information  to  a  depth  of  almost  5  rn  into 
the  sub-bottom. 


10  30  100  300 


a  (cm.) 

Figure  100.  Data  from  Puerto  Escondido,  Mexico. 
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B-5b  B-3 


Figure  101.  Data  interpretation  —  seafloor  array 
Puerto  Escondido  test  site. 

The  tests  in  Monterey  Bay  using  the  in  situ  expandable  survey  were 
run  in  water  6  m  deep,  at  a  temperature  of  15°C  and  salinity  of  34  0/00. 
The  formation  factor  vs  electrode  separation  curve  of  figure  39  shows  an 
initial  decrease  in  apparent  resistivity  with  depth.  A  core  (vibracorer) 
of  1.7  m  taken  at  this  location  showed  compact  silty  sand  overlying,  a 
loosely  packed  layer  of  sand  and  gravel  beginning  30  cm  below  the  surface. 
The  sand  and  gravel  graded  into  well  compacted  sandy  silt  over  the  next 
meter.  Unfortunately,  the  sediment  was  disturbed  considerably  during  the 
the  process  of  coring  and  handling,  so  no  resistivity  measurements  on  the 
core  itself  were  possible.  It  is  reasonable  to  assume,  however,  that  the 
resistivity  of  the  loosely  packed  sand  and  gravel  was  less  than  that  of 
the  fairly  compact  overlying  silty  sand,  and  that  of  the  well -compacted 


sandy  silt  beneath.  This  interpretation  agrees  with  that  obtained  by 
curve  matching,  shown  on  figure  9  (section  3.3)  and  figure  102.  The 
high-resistivity  layer  at  a  depth  of  1.7  in  may  have  been  the  cause  of 
the  core  tube  stepping  at  that  depth. 

Surface  Towed  Array.  As  mentioned  in  section  3.3,  a  surface-towed 
horizontal  resistivity  array  could  provide  a  rapid  method  of  measuring  bot¬ 
tom  resistivity  values  over  large  ore as  of  the  seafloor.  A  surface-towed 
array  cannot  offer  the  resolution  of  bottom  structure  available  from  a 
bottom-deployed  array,  but  if  electrode  spacings  are  large  enough  to  rorce 
a  sufficient  amount  of  current  into  the  bottom,  and  if  the  impressed  cur¬ 
rent  is  large  enough  to  produce  signal  levels  well  above  the  background 
noise,  considerable  information  about  bottom  properties  may  be  obtained. 

This  is  apparent  from  figures  102  and  103.  The  ratio  aV/I  (the 
voltage  appearing  between  the  potential  electrodes  divided  by  the  impressed 
current  is  plotted  against  the  formation  factor  of  the  bottom,  with  the  sea 
floor  considered  as  a  single  homogeneous  layer  and  the  ratio  of  electrode 
separation  a  to  water  depth  d  is  equal  to  4/3.  For  an  impressed  current 
I  of  3  amps,  a  change  in  sea-floor  formation  factor  from  2  to  3  will  re¬ 
sult  in  a  change  of  about  0.7  millivolts  in  aV,  and  a  formation  factor 
change  from  9  to  10  gives  a  change  of  about  0.2  millivolts  in  aV.  A  0.2 
millivolt  change  is  detectable  in  practice,  so  unit  changes  in  sea-floor 
formation  factor  should  be  distinguishable  using  this  electrode  arrange¬ 
ment. 

The  water  was  calm  near  shore,  buc  waves  increased  to  about  5  ft 
(1.5  cm)  in  deeper  water,  resulting  in  considerable  ship  motion  (which  did 
not  impede  the  profiling  operation).  Water  surface  temperature  was  15°C, 
and  salinity  34  0/00,  over  the  entire  area,  giving  resistivity  of  0.24 
ohm-meters.  Time,  core  site,  current  I,  and  water  depth  were  noted  on 
the  chart  recorder  output  of  aV  (fig.  23). 

Previous  measurements  of  resistivity  over  water-covered  areas  have 
been  by  Schlumberger  and  Leonardon  (1934);  Volker  and  Dijkstra  (1955);  and 
Bonlos  (1972),  but  these  were  done  with  stationary  rather  than  towed  ar¬ 
rays.  Unz  (1959)  discusses  some  of  the  theory  of  water-surface  resistiv¬ 
ity  measurements. 
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Figure  102.  Data  interpretation  —  seafloor  array 
Monterey  Bay  test  site. 
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Figure  103.  A V/I  vs  formation  factor  of  seafloor  sediments 


: 


In  contrast  to  the  bottom-deployed  array,  data  obtained  from  a 
surface -towed  array  is  handled  the  same  as  data  obtained  from  an  identi¬ 
cal  array  on  land.  Differences  between  the  land  and  marine  cases  arise 


mainly  during  actual  application: 

(1)  Compared  to  most  land  situations,  the  resistivity  of  the 
first  (salt  water)  layer  in  the  marine  case  is  very  low, 


so  a  large  fraction  of  the  impressed  current  will  flow 
through  the  water  layer.  Therefore,  higher  current  lev¬ 
els  must  be  used  in  offshore  work  to  assure  that  suffi¬ 
cient  current  will  penetrate  the  seafloor  to  give  accu¬ 
rate  readings. 

The  contact  resistance  of  sea  water  is  much  lower  than 
that  of  most  earth  materials,  so  high  voltages  are  not 
needed  to  force  a  large  current  flow  into  sea  water. 
dhile  several  hundred  volts  may  be  required  to  force  one 
amp  into  the  earth,  twelve  volts  are  sufficient  to  force 
10  amps  into  sea  water,  using  electrodes  of  the  same  size. 
Noise  sources  for  the  land  and  marine  cases  are  somewhat 
different  in  magnitude  and  origin.  While  the  use  of  com¬ 
mutated  direct  current  is  helpful  in  removing  the  effect 
of  dc  drift  of  the  potential  electrodes  on  land,  exces¬ 
sive  changes  ot  dc  level  often  are  difficult  to  compen¬ 
sate.  At  sea,  the  drift  of  the  silver-silver  chloride 
potential  electrodes  is  very  small,  but  the  effects  of 
waves,  currents,  and  other  marine  phenomena  may  be 
troublesome.  Typically,  the  noise  level  of  a  towed 
electrode  pair  is  less  than  a  few  tenths  of  a  milli¬ 
volt. 

In  contrast  to  the  onshore  situation,  the  offshore 
resistivity  profile  is  continuous  with  every  point 
along  the  ship's  course  sampled.  Whereas  one  or  more 
electrodes  (usually  two  or  all  four)  must  be  dug  up 
and  re-implanted  for  each  resistivity  measurement  on 
shore,  all  electrodes  are  in  continuous  contact  with 
the  water  in  the  offshore  case.  This  allows  much 
faster  (and  thus  more  economical)  coverage  offshore, 
as  a  towing  speed  of  at  1  feast  3  knots  may  be  maintained 
in  most  sea  states,  with  speeds  of  up  to  10  knots  feasi¬ 
ble  in  some  situations . 


7.3.2  Results 


The  towed  resistivity  data  was  analyzed  by  the  following  procedure: 
(1)  The  apparent  resistivity  pa  is  calculated  using  the 


equation 


_  o  a  AV 
pa  -  2ira  T 
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where  a  is  the  electrode  separation  (20  m),  aV  is  the 
measured  potential  difference  at  th^  core  site,  and  I 
is  the  measured  current. 

(2)  Using  the  calculated  value  of  pa  and  the  known  water 
depth,  the  sediment  resistivity  ps  may  be  found  by 
reference  to  a  standard  set  of  two-layer  resistivity 
curves,  such  as  those  given  by  Mooney  and  Wetzel  (1956a) 
or  Van  Nostrand  and  Cook  (1966,  p.  91). 

(3)  The  formation  factor  (F)  is  calculated  by  dividing  the 

sediment  resistivity  ps  by  the  water  resistivity  Pll 
(0.24  ohm-m) :  w 


The  data  is  summarized  in  Appendix  A-5.  The  data  was  recorded  con¬ 
tinuously  (except  for  short  periods  when  the  current  was  turned  off  to 
check  the  zero  level  of  the  potential  electrodes).  For  analysis,  however, 
the  data  has  been  tabulated  only  at  individual  core  sites.  The  values  of 
formation  factor  F  are  considered  accurate  to  ±  10  percent,  except  at  the 
deepest  core  sites  (71,  71,  83,  and  84),  where  the  amount  of  current  pene¬ 
trating  the  bottom  may  have  been  insufficient  for  accurate  readings. 

The  surface  towed  resistivity  data  shown  in  table  11  covers  the  sur¬ 
vey  area  shown  in  figure  104.  The  data  is  shown  graphically  in  figures  105 
through  107.  Figure  105  shows  the  formation  factor  (F)  and  the  porosity  (<j>) 
measured  on  core  samples  for  the  three  east-west  traverse  lines  shown  in 
figure  106.  The  expected  inverse  correlation  between  F  and  <{>  is  readily 
apparent  in  lines  124-117  and  107-99,  but  does  not  appear  on  line  116-108. 
The  overall  pattern  of  sediment  porosity  distribution  along  the  three  east- 
west  survey  lines  is  evident  in  figure  106:  soft  (high-porosity)  sediments 
toward  the  center  of  the  area,  grading  into  harder  material  at  the  eastern 
and  western  boundaries,  with  a  probable  bedrock  outcrop  in  the  northwestern 
corner.  This  pattern  also  may  be  seen  in  the  formation  factor  contour  map 

(fig.  107),  which  is  quite  similar  in  general  outline  to  the  porosity  con¬ 
tour  map  (fig.  108). 


160 


LOCATION  MAP 


LEGEND 


-  8-23-72 

.  8-24-72 

°  •  Care  Site  (  Diver  Operated  Carer ) 


Figure  104.  Index  map  of  towed  resistivity  survey 
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Figure  105.  Direct  current  resistivity  test  —  east-west 
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Figure  107.  Formation  factor  contour  map. 
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Figure  108.  North-south  formation  factor  and  porosity 
profiles  8ur face-towed  resistivity. 


Formation  factor  and  porosity  along  the  two  north-south  traverse 
lines  (118-71  and  121-72)  are  shown  in  figure  108.  The  inverse  correla¬ 
tion  between  F  and  $  again  is  apparent  on  these  profiles,  although,  as 
mentioned  above,  the  deep-water  values  may  be  somewhat  suspect.  The  low 
porosity  at  core  site  X-2  is  reflected  as  a  peak  in  the  formation  factor 
curve  at  that  point.  As  the  depth  at  X-2  was  64  ft  (19.5  m)  roughlv  equal 
to  the  electrode  separation  a,  it  appears  that  valid  data  may  be  obtained 
at  a  water  depth  at  least  equal  to  the  electrode  separation. 

A  statistical  analysis  of  formation  factor  and  porosity  over  the 
entire  survey  area  showed  no  mathematical  correlation  between  the  two 
parameters,  although  reasonable  correlation  along  individual  survey  lines 
is  evident  in  figures  105  through  108.  There  are  two  possible  explanations 
for  the  lack  of  mathematical  correlation: 

(1)  It  is  quite  possible  that  the  porosity  measured  on  some 
of  the  cores  was  not  representative  of  true  in  situ  por¬ 
osity.  The  sediments  are  disturbed  to  some  extent  when 
the  core  is  taken,  and  the  degree  of  disturbance  probably 
varied  from  core  to  core.  A  few  "wild  points",  due  to  er¬ 
roneous  porosities,  could  have  a  strong  negative  effect 

on  the  mathematical  correlation  coefficient. 

(2)  The  diver-operated  coring  device  has  a  maximum  penetration 
of  70  cm,  and  most  cores  in  the  survey  area  did  not  achieve 
maximum  penetration.  The  electric  current  penetrates  the 
bottom  to  a  depth  considerably  greater  than  one  meter,  and 
so  measures  sediment  properties  at  depths  greater  than  those 
sampled  by  the  cores.  If  sediment  porosity  changes  signifi¬ 
cantly  with  depth  at  a  given  location,  or  if  bedrock  is 
within  a  few  meters  of  the  sediment  surface,  this  informa¬ 
tion  will  be  reflected  in  the  towed  resistivity  values,  but 
may  not  be  evident  in  the  cores. 


8.  CONCLUSIONS  AND  RECOMMENDATIONS 
8.1  Reflection  Coefficient  Mapping  Experiment 

The  correlation  matrices  and  other  results  of  this  experiment  indi¬ 
cate  that  there  is  some  potential  in  using  remote  sensing  acoustical  ap¬ 
paratus  by  a  ship  underway  for  mapping,  at  least  qualitatively,  the  mass 
physical  properties  of  seafloor  sediments.  Correlations  obtained  in  this 
study  were  inconclusive,  but  in-depth  evaluation  of  the  data  would  provide 
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information  on  alternate  approaches  to  the  problems  involved  in  data  ac¬ 
quisition  and  analysis.  It  was  obvious  from  the  preliminary  analysis  of 
the  results  that  there  are  many  problems  in  applying  surface  towed  systems 
to  map  sediment  properties,  even  with  the  utmost  navigation  control  and 
advanced  acquisition  instrumentation.  Improved  automatic  data  processing, 
using  real-time,  on-line  displays  is  absolutely  essential.  It  was  learned 
that  at  least  four  problems  exist  which  may  affect  the  reliability  of  the 
data: 

(1)  Determining  what  is  a  representative  statistical  sample. 

(2 )  Obtaining  accurate  near-surface  sediment  velocity  data. 

(3)  Control  of  the  distance  between  the  source-receiver  sys¬ 
tem  and  the  seafloor  to  maintain  the  same  differential 
area  that  the  energy  impinges  upon. 

(4)  Maintaining  the  stability  of  the  sensor  package  in  a 
horizontal  orientation,  for  attaining  normal  incidence. 

It  is  recommended  that  further  research  be  conducted  in  data  pro¬ 
cessing  and  analysis  procedures,  using  the  data  derived  from  this  study. 
Investigation  of  all  the  permutations  of  data  analysis  in  this  study  was 
beyond  the  capability  and  funds  allotted.  Development  should  be  pursued 
in  the  engineering  design  of  self-contained  submerged  tow  systems  capable 
of  powering  and  accurately  navigating  at  any  depth  and  at  a  fixed  eleva¬ 
tion  in  a  horizontal  orientation.  This  development  is  particularly  desir¬ 
able  in  view  of  the  requirements  for  more  information  on  deep  sea  sediment 
properties  and  the  need  to  obtain  these  measurements  in  as  rapid  and  eco¬ 
nomical  manner  as  possible. 

8.2  Shear  Wave  Experiment 

In  the  past,  the  detection  of  the  shear  wave  in  saturated  marine 
sediments  was  accomplished  indirectly  by  measuring  the  Stoneley  wave  velo¬ 
city  or  employing  a  viscoelastometer.  The  high  energy  losses  of  the  shear 
wave  and  the  low  sediment  rigidity  prevented  direct  measurements  of  the 
shear  wave.  It  Is  felt  that,  on  the  basis  of  the  data  and  discussions  pre¬ 
sented,  the  described  seismic  method  successfully  detected  the  direct  SH 
wave  in  gravelly  coarse  sand  at  a  lake  test  site.  For  most  of  the  stations 
in  the  Monterey  Bay  test  area,  it  could  not  be  clearly  established  that  the 
detected  wave  was  the  shear  wave  or  the  Love  wave,  which  has  been  attributed 
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to  insufficient  core  lengths  to  assertain  sediment  stratigraphy.  As  dras¬ 
tic  changes  in  the  sediment  properties  and,  thus,  shear  velocities  are  not 
generally  anticipated,  and  because  the  separation  distances  between  source 
and  receivers  is  small  the  Love  wave,  if  detected,  would  have  comparable 
velocity  to  the  shear  wave.  Therefore  it  was  assumed  that  the  velocity 
obtained  is  the  shear  velocity. 

A  refined  seismic  method  is  planned  for  future  in  situ  tests.  A 
multi-channel  magnetic  tape  recorder  will  be  used  instead  of  the  oscillo¬ 
scope.  This  will  increase  the  accuracy  in  time  measurements  and  reduce 
the  number  of  source  shots.  Hydrophones  will  be  employed  to  detect  the 
compressional  wave.  Moreover,  both  receivers  probes  and  source  probe  will 
be  inserted  mechanically  to  obtain  accurate  distances  between  probes.  A 
sampler  capable  of  taking  longer  cores  is  also  a  necessary  addition. 

The  source  proved  adequate  in  generating  shear  waves;  however,  the 
air  gun  will  introduce  air  into  the  sediment.  Some  air  might  enter  the 
sediment  region  between  the  source  and  receivers  although  the  source  en¬ 
ergy  is  directed  away  from  this  sediment.  This  will  reduce  the  value  of 
compressional  wave  velocity  largely  through  the  reduction  in  the  bulk 
modulus.  Therefore,  consideration  has  been  given  to  develop  a  new  source 
such  as  a  bi-directional  electro-magnetic  hammer.  Being  bi-directional, 
the  new  source  can  produce  stacked  records  to  magnify  the  waves  or  rec¬ 
ords  showing  the  polarity  reversal  of  the  waves  through  reversing  the 
direction  of  the  source. 

As  discussed,  the  dynamic  values  of  rigidity,  bulk  modulus,  Young's 
modulus,  and  Poisson's  ratio  can  be  computed  from  known  values  of  the  com¬ 
pressional  velocity,  shear  wave  velocity,  and  density.  As  the  time  dura¬ 
tions  of  acoustical  stresses  last  only  for  a  thousandth  of  a  second,  the 
dynamic  constants  will  vary  considerably  with  those  determined  by  conven¬ 
tional  static  tests  in  the  laboratory.  Clearly,  in  situ  acoustical  mea¬ 
surements  can  provide  useful  information  about  these  constants,  particu¬ 
larly  if  acoustical  properties  can  be  related  with  those  resulting  from 
static  laboratory  tests. 

Empirical  relationships  do  exist  between  acoustical  and  mass  physi¬ 
cal  properties.  These  relationships  can  be  used  in  two  ways.  If  porosity, 
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textural  properties,  shear  strength,  bearing  capacity,  or  settlement  can 
be  related  to  the  acoustical  properties  for  a  particular  marine  environ¬ 
ment  (bay,  shelf,  abyssal  pi ai ns/hi  1 1 s ) ,  one  need  only  secure  the  acousti¬ 
cal  properties  and  a  disturbed  core  sample  to  ascertain  the  sediment  type 
to  correlate  these  relationships  and  obtain  the  approximate  mass  physical 
and  engineering  properties  of  the  site.  Reversing  the  procedures,  one 
could  obtain  the  acoustical  properties  by  consulting  these  relationships 
if  the  mass  physical  properties  are  measured  from  an  undisturbed  core  sam¬ 
ple.  Future  work  should  entail  securing  these  relationships  and  determin¬ 
ing  the  accuracy  that  can  be  expected  in  using  these  relationships. 

In  conclusion,  knowledge  of  the  acoustical  properties  of  seafloor 
sediments,  as  derived  from  shear  wave  analysis,  can  provide  useful  infor¬ 
mation  in  obtaining  the  geotechnical  variations  within  the  marine  environ¬ 
ment.  This  information  would  be  extremely  helpful  in  solving  engineering 
problems  related  to  excavation,  tunneling,  or  for  construction  of  platforms, 
dams,  and  bridges  in  the  marine  environment. 


8.3  Direct  Current  Resistivity  Experiment 

From  the  above  discussion  of  results,  it  is  evident  that  a  towed 
offshore  resistivity  array  can  give  information  about  seafloor  properties, 
even  with  the  use  of  unsophisticated  equipment.  Several  improvements  in 
equipment  and  data  handling  could  result  in  even  more  rapid,  more  accurate, 
and  less  expensive  resistivity  data  acquisition: 

(1)  Other  electrode  arrays,  such  as  the  Schlumberger  and  Ani- 
pole,  should  be  investigated  theoretically  in  the  field. 

It  should  be  possible  to  optimize  array  configuration  and 
electrode  spacing  for  expected  water  depths  and  bottom 
properties. 

(2)  The  use  of  a  multiple  electrode  array  would  enable  multi¬ 
ple-layer  interpretations  of  bottom  geology  to  be  made. 
Electrodes  could  be  switched  sequentially  to  vary  their 
spacing,  or  arrays  could  be  devised  in  which  a  single  pair 
current  electrodes  is  continuously  energized  and  the  sig¬ 
nal  between  multiple  potential  electrode  pairs  recorded. 

(3)  Improved  accuracy  and  speed  would  result  from  electronic 
processing  of  the  data,  with  the  current  signal  divided  elec¬ 
tronically  into  the  filtered  potential  signal.  This  infor¬ 
mation  could  be  stored  on  tape,  along  with  water  depth,  tem¬ 
perature,  and  salinity,  for  automatic  processing  of  the  data. 
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In  summary,  the  results  of  this  study  indicate  that  a  towed  off¬ 
shore  resistivity  system  can,  especially  when  used  in  conjunction  with 
other  geophysical  techniques,  provide  valuable  data  about  seafloor  prop¬ 
erties.  System  improvements,  such  as  those  suggested  above,  could  make 
the  technique  even  more  accurate  and  economical. 
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APPENDIX  A.l.  CALCULATIONS  FOR  UNDERWAY  ACOUSTICS  EXPERIMENT 
A. 1.1  Calculation  of  Reflection  Coefficient  (R)  for  Sparker  (0-5  kHz  Source) 


I- 


t  =  time 

a  =  Sound  velocity 

u  =  Particle  motion  amplitude  or  pressure 


Streamer  hydro  signal 

e-m  couplinq- 
pulse 


4- 


direct  arrival 

- time 

bottom  Rf 1 . 


Assuming  a  point  source: 


•u  -  u  •  const ■ 

monitor  usource  ’  d  -d„ 


m  s 

Assume  d$  and  df  «  D  so  source  to  receiver  reflector  path  length  is  approximate 


(A-2) 


V(D-ds)2  +  (•C/2)2  +  ^(D-dr)2  +  U/2)2  =  reflection  path. 
So  almost  exactly 


reflector  path  =  2  V(  D-ds)2  =  U/2)2 

=  V4(°-ds)2  +  t* 

Note:  In  deeper  water  £«D  so  reflector  path  2D  -  (d  +  d  ). 
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rcvr  source 


/  4(D-ds)2  +  l2 


Eliminating  usource  between  (1)  and  (2) 


u  ..  .  a/4  (D-d  )2  +  l2 

momtor  1  v  '  s' 


Solving  for  R, 


/4 ( D-d  )2  +  l2  u 
f__  s'  #  rcvr 

dm  "ds  umonitor 


emonitor  umonitor  ’  ^impulse  hydro 


ercvr  urcvr  *  ^streamer  hydrophone 


•  urcvr  _  ercvr  .  Gimp.  hyd. 
umonitor  emonitor  Gsta.  hyd. 


Eliminating  u  /umon1tor  between  (A-5)  and  (A-6)  we  have  for  R  for  the  sparker  system: 


/4 (D-ds ) +  l  Gimp.  hydro. 


Streamer  hydro  monitor 


geometrical 

factor 


hydrophone  ratio  of  VoHages 

gain  ratio  recorded  on  Mag. 

tape,  i.e.,  Sr/Sx 


Substituting  the  numerical  data  from  below  we  have 


R  =  .532  \/402  +  832  •  (.0233) 


Numerical  Data: 


Sparker : 

ds  =  10"  =  .832' 
dm  =  2'  8V  =  2.71' 
dr  =  15". 

Pinger: 

d  =  22" 

Constants:  V  =  4.66  ft/m-sec  t  =  6.2  m-sec 
,\Z  =  vt  =  28.8' 


- =  .0233 

“streamer  hyd. 


Geometrical  Factor  = 


_V4(D-.832)2  +  (28.8) 


=  .532  V4°2  +  832  • 


Streamer  Calibration  Flow  Chart: 


.5  rnv 


Rcvr 


MC-115  (Standard) 


A.  1.2  Calculation  of  Reflection  Coefficient  for  Pinger  (12  kHz  Source) 


u  =  particle  amplitude  or  pressure 

In  the  pinger,  unlike  the  sparker,  the  source  and  receiver  are  the  same  transducer;  how¬ 
ever,  to  aid  analysis  assume  they  are  two  separate  transducers. 


Also  assume  that  the  calibrated  hydrophone,  MC-115,  is  located  at  the  reference  point  3  feet 
below  the  pinger.  From  previous  tank  calibrations,  when  ex  =  183  v.,  the  MC-115  develops  lv. 
The  pinger  receiver  in  the  same  location  would  develop  4.65. 

Assuming  a  point  source. 


ux  =  c  3~ 


( A- 13 ) 


ur  =  c 


21^37 


(A-14) 


Eliminating  c  u 


ur  =  ux  mzr  * R  • 


(A-15) 


Let  ur  =  er/4.65  and  ux  =  ex/183. 


Then 


=  C®3]  [2 (D-d)]  4,65  R- 


(A-16) 


183 


\ 


APPENDIX  A. 2  DERIVATION  OF  BASIC  EQUATIONS  OF  LINEAR  VISCOELASTICITY 

(after  Borcherdt,  R.  D.,  1972) 


Preceding  page  blank 
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APPENDIX  A. 2  DERIVATION  OF  BASIC  EQUATIONS  OF  LINEAR  VISCOELASTICITY 
(after  Borcherdt,  R.  D.,  197?) 

A. 2.1  Body  Waves  In  a  Linear  Viscoelastic  Medium 

This  section  Is  IncludpJ  for  completeness,  although  a  detailed  study  will  not  be  given. 
The  development  of  the  basic  equations  will  follow  the  approach  giver  by  Borcherdt  (1972), 
and  the  reader  Is  referred  to  that  thesis  for  a  more  comprehensive  study. 

Viscoelasticity  Is  a  property  of  solids  and  liquids  which  exhibit  both  viscous  ano  elas¬ 
tic  behavior  under  deformation.  The  dissipation  and  storage  of  mechanical  energy  Is  observed 
from  this  behavior.  In  addition,  the  constant  that  links  stress  and  strain  In  the  theory  of 
elasticity  becomes  a  time  dependent  function  In  the  constitutive  equation  of  viscoelastic 
theory.  For  sufficiently  small  strains,  the  behavior  of  a  viscoelastic  material  can  be  des¬ 
cribed  by  the  linear  theory  of  elasticity.  Employing  the  superposition  principal,  the  stress 
strain  relation  can  be  written  as  follows: 


V'1 


f 


Wt_T)  dek£(l) 


(A-19) 


where  j ( t ) 

eu(t) 

r1jk£ 


=  components  of  the  2nd  order  time  dependent  stress  tensor, 

:  components  of  the  2nd  crder  time  dependent  strain  tensor,  and 

=  components  of  the  4th  order  time  dependent  relaxation  function. 


This  equation  can  be  written  In  a  more  condensed  form  as  a  convolution: 


P1j  =  r1  jk£  dek£  *  (A_20) 

If  the  material  is  Isotropic,  a  further  simplification  can  be  made: 

Pjj-r,  \ti  (A-il) 

Pkk  '  rk  dekk  (#-22’ 


where 


The  primes 


rs  =  relaxation  function  describing  shear  behavior  and 
rk  =  relaxation  function  describing  bulk  behavior, 
denote  the  deviator  components  of  stress  and  strain: 


P1j  =  P1j  ’  I61j  Pkk 
e1j  *  e1j  "  7  5ij  ekk’ 


(A-23) 

(A-24) 
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It  would  be  convenient  to  simplify  these  stress  strain  relationships  by  considering 
the  steady  state  problem.  Moreover,  It  will  be  assumed  sufficient  time  has  elapsed  so 
initial  conditions  may  be  neglected.  Let  P^  =  P^e1^  and  e^.  =  E^e1  *  where  P^  and 
are  complex  constants  independent  of  time.  Equations  (A-21)  and  (A-22)  can  now  be  written 


p-ij  =  1(0  RS(M  E1J 


(A-25) 


P-kk  =  iu)  Rk(1w)  Ekk 


(A-26) 


where  Rs  and  represent  the  Fourier  transforms  of  the  shear  and  bulk  relaxation  functions, 
respectively.  These  quantities,  R$  and  R^,  can  be  formulated  into  a  more  familiar  form: 


(A-27) 


(A-28) 


MVRS, 


(A-29) 


It  follows 


Pij  -  2  us  Eij  1  M 


(A-30a) 


P-kk  =  3  K  Ekk  • 


(A-30b) 


Other  useful  relationships  are 


Pij  =  fiij  »  Ekk  +  2  Eij 


(A-31) 


P-ij  "  6ij  (K  -  I  Ekk  +  2  11 


(A-32) 


Attention  will  now  be  focused  on  the  equation  of  motion  and  itrain  tensor, 
conservation  of  linear  momentum,  the  equation  of  motion  is: 


:rom  the 


PU.J  +  f1 


(A-33) 


where  f.  =  body  force  per  unit  volume, 

<p  =  mass  density,  and 
u^  e  component  of  the  displacement  vector. 


The  strain  tensor,  assuming  infinitesimal  displacement  is  given  by 


e1j  ‘  ?  (ui.J  *  uj.i)  -  (A-3' 

Steady  state  simplification  will  be  employed  again.  Let  u.  =  U1e1“)t  where  Ui  is  a  complex 
constant  independent  of  time.  Therefore,  the  equation  of  motion  and  strain  tensor  simplify 
to 


and 


P-1j,j  +  ♦  “  Ui  =  0 


(A-35) 


Eij  =  7 


(A-36) 


(where  body  forces  have  been  neglected).  Using  (A-36)  and  (A-32)  the  equation  of  motion, 
(A-35),  can  be  written  in  terms  of  u^  (in  vector  form): 

(<  +  J  u)v  6  -  pV  x(Vxu)  +  4>  w2  u  =  0  (A-37) 

where  / 

e  =  v  u  . 

Note  that  this  equation  is  equivalent  to  the  equation  of  motion  for  an  elastic  medium  except 
for  the  Lame  constants  which  may  now  be  complex  (i.e.,  y  =  pR  +  ipj  and  k  =  kr  +  i<j).  Also, 
the  expression  for  the  time  rate  of  energy  and  2nd  law  at  thermodynamics  reveal  that  pj  and 
Kj  are  non-negative  (Borcherdt,  1972,  p.  43). 

By  taking  the  divergence  and  the  curl  of  the  equation  of  motion,  (A-37),  one  can  ob¬ 
tain,  respectively,  the  equations  for  dilatational  and  rotational  waves: 


and 


v2e  +  —  e  =  o 

~  2 


(A-38) 


(A-39) 


where  the  complex  compressional  velocity  (a)  and  the  complex  shear  velocity  (p)  are  given  by 


F  - 


(A-40) 


B  =  P.V .  (H.)^  .  (A-41 ) 

(P.V.  represents  the  principal  value).  The  respective  velocities  of  the  shear  and  compres¬ 
sional  wave  in  homogeneous  isotropic  linear  viscoelastic  medium  are 


Vp  =  (A-42) 

vs  "  Hf5]"1-  (A-43) 
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The  absorption  coefficient  of  these  two  waves  are 


(A-44) 

(A-45) 


Little  will  be  said  concerning  the  quantity  Q'1  except  to 
specific  attenuation  factor  in  bulk  will  be  denoted  as  Q'1  and 


define  Q^1,  Qj1 ,  and  Q^1 . 
is  expressed  as 


The 


(A-46) 


Similarly,  the  attenuation  factor  for  general  dilatation  and  shear  is  given  respectively  as 


and 


+  T  UT 
n-i  =  J - £ — L 

WP  74 

^R  7  UR 

Ur 

Q->  =  -1  . 
s  UR 


(A-47) 


(A-48) 


It  is  convenient  to  express  the  compressional  and  shear  velocity  in  terms  of  Qp1  or  Qs*. 


From  (A-43 ) 


or 


^  ^  Re  ^  ^W\.  . 

Further  reauction  using  (A-48)  yields 


5  ^ 

'Mr — — — - 


vr+  q-*  + 1 


similarly, 


Vi  +  Qp2  + 


4 

T  UR 


(A-49) 


(■  *50) 


(A-51) 


(A-52) 


This  concludes  the  discussion  on  the  basic  equations  of  body  waves  in  linear  viscoelastic 
medium.  For  further  discussions,  refer  to  Porcherdt  (1972). 
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We  can  check  the  correctness  of  the  formula  from  the  following  consideration.  The  term 
2 f D-d )  is  the  distance  from  source  to  receiver.  If  we  move  the  receiver  to  the  reference 
point  as  shown  this  term  equals  3'  and  there  is  no  reflection  (R=l).  Then  (A-16)  becomes: 


for  ex  =  183v.  This  yields  er  =  4.65  volts  which  is  the  value  from  tank  measurements.  Solv 
ing  for  R  with  d  =  22"  =  1.83' 


R  =  26.2  (D1  -  1.83') 


ercvr 

cxmtr 


( A- 18 ) 


Tank  Calibrations  of  Pinger  —  Block  Diagram: 


X  mtr 


183  V 


T 


3' 

IV 


Pinger  Xmtr 


MC  -  115 


Rcvr. 


MC-II5 


Source 


IV 


-I  \  \ 


J-9 


Pr.nger 

Rcvr. 


4.65  V 
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3  SEDIMENT  PROPERTY  EQUATIONS 


191 


Property 

Shear-wave 

dissipation 

function 


APPENDIX  A. 3  SEDIMENT  PROPERTY  FILIATIONS 
Symbol 


Method  of  Calculation 


Q"1 

vs 


2  V  a 

q-i  =  — Li. 

VS  u 


Rigidity 

( \\ 2 

representing 
elastic  response 

UR 

WR  = 

M) 

Riqidity 

:  yR  V 

representing  wave 

UI 

UI  = 

energy  damping 

Bulk  modulus 

A 

representing 

kr 

kr 

:  ♦  Vp  ~  JVR 

elastic  response 

Lame  constant 

2 

representing 
elastic  response 

XR 

XR  = 

!  KR  '  7  UR 

Poisson's  ratio 

<V=A's)2  -  2 

0 

0  = 

2(Vpqs/Vs)2-  2 

Young's  modulus  E  E  =  uR(l+a) 
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A-3.2  Notation  Summary 


0 


ft 


Symbol 


$ 


vl 

yR 

a 

w 


_ Description _ _ 

Shear-wave  attenuation  coefficient 
Compressional -wave  attenuation  coefficient 
Young's  modulus 

Bulk  modulus  representing  wave  energy  damping 
Bulk  modulus  representing  elastic  response 
y2(l+Q-yj/v/l+QpT  +  1 

J2(i*q;2)/i/i*q;t  +  i 

Shear-wave  dissipation  function 
Compressional -wave  dissipation  function 
Snear-wave  velocity 
Compressional -wave  velocity 

Bulk  density 

Lame  constant  representing  wave  energy  damping 
Lame  constant  representing  elastic  response 
Rigidity  representing  wave  energy  damping 
Rigidity  representing  elastic  response 
Poisson's  ratio 
Angular  frequency 
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OF  SAMPLING  PROGRAM 


95 


Site 

No 

Corrected 
Water  Depth 
(ft) 

Core 

Length 

(cm) 

Predominant 

Sediment  Type 

Corer  Type 

Driver  Sampler: 
Piston  type-P 
Open  type-w 

23 

102 

60.4 

silty-sand 

P 

24 

106 

43.8 

silty-sand 

P 

25 

106 

63.2 

sandy-sll t 

P 

26 

113 

65.1 

silty-sand 

P 

38 

124 

61.7 

s  1 1  ty-sand 

P 

39 

117 

62.2 

silty-sand 

P 

40 

111 

65.6 

silty-sand 

P 

41 

105 

60.9 

s 11 ty-sand 

P 

42 

99 

65.0 

silty-sand 

P 

43 

100 

62.9 

s  11  ty-sand 

P 

44 

97 

63.6 

silty-sand 

P 

45 

92 

64.0 

silty-sand 

P 

46 

92 

65.4 

silty-sand 

P 

47 

95 

63.3 

silty-sand 

P 

48 

96 

61.7 

silty-sand 

P 

49 

1C2 

62.0 

silty-sand 

P 

50 

105 

27.9 

silty-sand 

P 

51 

114 

54.5 

silty-sand 

P 

52 

121 

47.9 

sandy-silt 

P 

57 

123 

42.2 

silty-sand 

P 

58 

118 

55.6 

silty-sand 

P 

59 

111 

62.3 

silty-sand 

P 

60 

105 

64.4 

sllty-sand/sandy-sllt 

P 

61 

96 

65.2 

silty-sand 

P 

62 

93 

62.9 

silty-sand 

P 

63 

91 

62.3 

sllty-sr.nd 

P 

64 

89 

40.2 

silty-sand 

P 

65 

87 

60.9 

silty-sand 

P 

66 

85 

51.1 

silty-sand 

P 

68 

83 

13.8 

silty-sand 

P 

69 

83 

45.2 

silty-sand 

P 

69 

83 

90.2 

silty-sand 

0 

70 

83 

62.6 

silty-sand 

P 

71 

84 

63.8 

silty-sand 

P 

72 

85 

63.2 

silty-sand 

P 

73 

90 

54.5 

silty-sand 

P 

73 

90 

104.0 

silty-gravelly  sand 

0 

81 

73 

6.3 

silty-sand 

P 

82 

70 

32.6 

silty-sand 

P 

83 

68 

4.5 

sandy-silt 

P 

84 

55 

54.2 

silty-sand 

P 

85 

19 

63.7 

silty-sand 

P 

86 

73 

52.2 

silty-sand 

P 

87 

75 

10.4 

silty-sand 

P 

88 

73 

11.0 

silty-sand 

P 

92 

45 

9.0 

silty-sand 

P 

93 

48 

63.9 

silty-sand 

P 

94 

53 

59.7 

sandy-silt 

P 

95A 

56 

38.7 

sandy-silt 

P 

958 

56 

5.1 

sandy-silt 

P 

196 
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Corrected 

Core 

Corer  fype 

Predominant 

Driver  Sewpltr: 

”  ,  Water  Deoth 

No  (ft) 

Length 

(cm) 

Sediment  Type 

Piston  typ«-P 
Open  type-0 

100 

43 

34.1 

gravelly-sand 

P 

102 

47 

39.9 

silty-sand/sandy-sllt 

P 

103 

48 

51.0 

sandy-sllt 

P 

103 

48 

38.1 

sandy-sllt 

0 

104 

49 

63.0 

silty-gravelly  sand 

P 

105 

49 

41.2 

silty-gravelly  sand 

P 

109 

38 

54.3 

silty-sand 

P 

109 

38 

36.3 

silty-sand 

0 

110 

40 

51.2 

silty-sand 

p 

111 

42 

23.7 

sandy-sllt 

p 

112 

42 

41.7 

silty-sand 

p 

113 

43 

59.7 

silty-gravelly  sand 

p 

114 

45 

51.7 

silty-sandy  gravel 

p 

115 

46 

6.0 

silty-sand 

p 

118(1) 

32 

47.2 

silty-sand 

p 

118(2) 

32 

56.8 

silty-sand 

p 

119 

33 

14.9 

silty-sand 

p 

120 

34 

22.4 

silty-sand 

p 

121 

34 

54.9 

silty-sand 

p 

121 

34 

31.4 

silty-sand 

0 

122 

36 

55.8 

silty-sand 

p 

123 

39 

52.7 

silty-sand 

p 

123 

39 

34.3 

silty-sand 

0 

124 

38 

24.1 

silty-sand 

0 

IX 

56 

15.2 

silty-sand 

p 

2X 

60 

31.2 

silty-sand 

p 

3X 

63 

6.0 

silty- sand 

p 

l.H.l 

17 

64.6 

silty-sand 

p 

T.H.2 

17 

62.8 

silty-sand 

p 

S.C.  Buoy 

17 

62.0 

gravelly-silty  sand  and  sandy-sllt 

p 

S.C.  8uoy 
7/10 

17 

259.0 

gravelly-silty  sand 

0 

S.C.  8uoy 
8/28 

17 

170.2 

silty-gravelly  sand 

0 
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APPENDIX  A. 5  SUMMARY  OF  ACOUSTIC,  ELECTRICAL  RESISTIVITY 
AND  SEDIMENT  ANALYSES  DATA 
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APPENDIX  A. 6  SEDIMENT  ANALYSES 
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APPENDIX  A. 


7  DATA  SUMMARY  SHEETS  -  REFLECTION  COEFFICIENT 
MAPPING  EXPERIMENT 
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Best 

Available 

Copy 


^ I 


VJ1  OlCJD 


DATE  06/?W/ 


18.45.5-* 

18.46.57 


18.57.51 

18.58.51 


19.03.5"* 


DATE  06/30/72 


.  c 


Reproduced  from 
best  available  copy. 


TIMF 

pi  f' :i  rH*  .|VF.L 
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.  1  ■> 

6 . 6 

10.12.53 

.  /I  /•  - 

.  1  3 

6.5 

19.13.53 

.  4  .<  -• 

.3  17 

6.4 

19.14.57 

,4 

•  n : 

6.1 

19.15.57 

. r , 

. 

5.7 

19.JLA.tiI  ... 

-  -i.i . . 

.  1  >* 

0 . 6 

19.17.53 

•  **  *■" 

j.9 

19.05.57 

•  *’  i 

.  7  V 

5.9 

1°. 06.53 

.)r 

6.3 

19.07.53 

/  . 

.  1  lr 

6 . 4 

I9.ne.53 

c 

<n  0 

•  / 

6.1 

19. r9.p3 

.t;  7 

1  . .  "• 

5.3 

19.10.53 

•  '< 

5 . 6 

SP  A.RK  £Q  fHMA'cl  4  T  F  R  DcPTM 


AVE 

S I  GH  A 

a  -  L  r  S  c 

A  VF 

8  I0,MA 

.590 

.158 

4.9 

137.5 

.  7  __ 

."577 

.097' 

4.9 

1J6 .9 

.  0  ' 

.850 

.164 

4.3 

108.8 

1.9 

.879 

.147 

0 . 2 

1  J5.6 

.5 

.  "'61 

.  1  67 

5 . 6 

1  J  4 . 7 

.6 

.563 

.  074 

4.8 

1  j  7 . 4 

1.6 

•  62  A 

.  1  58 

4.7 

1.97.1 

.8 

.C1 7 

.169 

6 . 4 

1  0  7 . 4 

.8 

.  1  5e 

4.9 

1  ,j1  .  n 

.5 

.C57 

.138 

5 . 6 

101.3 

3.9 

.  8|>  7 

.  050 

4.7 

1  0  0 . 8 

.8 

.  *7,0 

.  r  v  f 

4 . 1 

10°.? 

.7 

.  r'  L  1 

.217 

7.4 

71.0 

.  7 

.588 

.090 

5.2  ' 

VH  .  7 

"  '  '  .3 

•  *  63 

.1  J7 

3.7 

y  8 , 8 

.6 

.c54 

.19  7 

0 . 2 

98.0 

.5 

.651 

.3«1 

9.2 

47.4 

.  5 

. c  96 

.127 

4  7 

93.3 

.5 

.867 

,14/> 

->  .  4 

917 .5 

5 . 2 

.690 

.130 

3 . 4 

9*  .  7 

.5 

.64* 

.  1  45 

4.0 

98.6 

1.6 

.705 

.126 

3.2 

45.6 

.8 

.60  3 

.947 

4.3 

94.? 

4.3 

.*1* 

.12° 

4.4 

9C  .  0 

.7 

.625 

.701 

3.4 

97.6 

2  9 

.629 

.161 

4 . 6 

9-1.5 

.  6 

.  *59 

.93^ 

0 . 6 

97.5 

2.7 

.471 

.  1  8f 

7.2 

93.9 

.5 

.  *  3° 

.  1  i>7 

4  .  { 

97.6 

.6 

.*31 

.717 

4.7 

V7 .5 

.6 

.'o7 

.  1 25 

3 . 6 

y7 .  n 

.6 

,72* 

.149 

1 / 

97.5 

1.0 

.  143 

.047 

17.3 

97.5 

.9 

.133 

.  036 

17.8 

45.1 

1.9 

.  1  31 

.  ('6  8 

M.4 

91.5 

1.1 

.14° 

.061 

1.  6 .  ? 

97.7 

1  •  4 

.157 

.983 

1  0 .  7 

97.4 

1.1 

.tbi 

.05a 

lo.6 

47 .  n 

1.7 

.159 

.  "58 

16.4 

97 .  n 

.5 

.  1  45 

.  08c 

’5.7 

9^.6 

4.0 

.121 

.050 

16.9 

91  .  7 

.5 

.123 

.028 

18.4 

•■1.3 

.5 

.135 

jJ'S? 

17.6 

yi  .  7 

.6 

.15? 

.108 

15.7 

■M  .  9 

.4 

.  1  9  7 

.161 

1 1 . 2 

.6 

1  -  C 

.1  j3 

.10  6 

1->.9 

97.5 

1.  J 

.144 

.03s 

17.0 

47.1 

l.B 

.  1  21. 

.0  71 

1.3.3 

97.6 

1.0 

•  1 

.15° 

_  iUJJ. _ 

97.8 

.1-4 

.160 

.  0c  4 

1  0 . 5 

97.1 

1.1 

253 


>viwiojisit-*^h‘'£ooo-oO'Ui'sio>’fls/Cy(ai(>  oiw  Jiuivrocwsisiaioji  ai'ai  o  a  aio  o-isi 


i 


DATE  07/?7/7?  ROLL* 0  ?9 


P1\R-rR  rHAvlVEL 
M'r  S  I G  M  A  8-LOSS 


SPARKER  CHANNEL _ 

AVE  "  SIGMA  S-LOSS 


WATFR  DEPTH 
A VF  SIGMA 


19.05.57 
19 . 0  6.57 
19.07.5’ 
19.08.57 
19^09.57 
19.10.5’ 

19.11.57 
19.12.5’ 
19.13.5’ 

19.14.57 
19.15.5’ 

19.16.57 

19.17.57 

19.18.57 

19.19.57 

19.20.57 

19.21.57 

19.22.57 

19.23.57 
19. 24.5’' 
19.25.5’ 

19 . 26 . 57 
19.27.5’ 

19.28.57 

19.29.57 

19.70.57 

19._31.57 

19.32.5’ 

19.33.5’ 


92.2 

92.2 

91.8 
92.2 

92.4 

92.1 
91  .7 

91 .9 

90.5 

91.1 
91.0 
91.0 
90.5 

90 . 7 

90.4 

“8977“ 

80.8 

89.8 

89.2 
89.0 

88.8 
~BW77 

38.5 

88.4 
88.1 
87.9 

87.4 

8773” 

87.0 


■  -•  * 


xj  ro  xj:  O  CX>lvJ1  Xj'  OJ  ro  A  Nil ho  VJl;  03 


1 


DATS  OR/'1!//''  r  I  LVI0  3? 


TIMP 

n  CMN 

'EL 

SPARKEP  CHANNEL 

WATPR  DEPTH 

AV  - 

0  I  iiM  A 

3-LOSS 

AVE 

SIGMA  B 

-L"SS 

AVE 

SIGMA 

17.17.53 

.701 

.06® 

3.1 

.606 

.094 

4.4 

124.7 

.4 

17.18.53 

.  6  2') 

.23° 

3.9 

.470 

.181 

7.5 

122.6 

5.2 

17.19.5'* 

.6®* 

.070 

3.3 

.694 

.157 

5.3 

123.2 

•  6 

1 7 . 2  0 . 5  3 

.*77 

.  in17 

3.5 

.640 

.  1 1? 

3.9 

122.5 

1.6 

17 . 21 . 57 

.*  V 

.  1  )* 

3.7 

.466 

.184 

7.6 

127.0 

.9 

17. 22.53 

.  7.17 

.121 

3.2 

.627 

.173 

4.9 

1  21 .4 

1.1 

17.23.5-* 

.14  7 

3.8 

.662 

.071 

3.6 

120. 0 

2.7 

17.24.57 

•  iu9 

3.8 

.581 

.537 

4.9 

119.8 

.6 

17.25.53 

.100 

3.4 

.  657 

.089 

3.7 

119.1 

.9 

17.26.53 

.67* 

.'’59 

3.5 

.60« 

.171 

5.2 

118.4 

.8 

17.27.53 

.67* 

.'57 

3.4 

.  616 

.077 

4.3 

117.7 

.6 

17.28.5'* 

,  f*  t\ 

1 1*5 

4.0 

.5  8* 

.168 

5.4" 

116.8 

.  4 

17. 29. 5-* 

.  a  7  •  ■ 

.  M  J 

3.4 

.579 

.138 

5 . 0 

116.1 

.7 

1 7 . 30 . 5  7 

.67  r. 

.  036 

3.4 

.  63n 

.164 

4.8 

115.5 

1.4 

17.31.57 

.6  3" 

.10* 

4.0 

.  597 

.106 

4 . 7 

114.7 

.8 

17.32.57 

.636 

.0  74 

4.0 

.486 

.218 

6.2 

113.9 

1.0 

17.33.57 

.  6  4" 

.  067 

3.9 

.528 

.^11 

6.9 

1 17 . 2 

.5 

1 7 . 3  4 , 5  7 

.  6  "  7 

.ot?7 

3.7 

.*27 

.248 

7.5 

117.5 

.8 

17.35.57 

*  .1 

.1  ?7 

4.7 

.46? 

.  1  69 

7-1 

111.4 

__LlZ_ 

17.36.57 

.  6  /  >: 

."77 

4.1 

.626 

.166 

4.9 

111.1 

.7 

17.37.57 

.64® 

."7* 

3.9 

.591 

.0  73 

4.6 

1.10.4 

.8 

17.38.53 

.6  37 

.  0  71 

4.0 

.555 

.139 

5.5 

ll]9.9 

.8 

17.39.57 

.616 

.  2  21 

4.5 

.549 

.153 

6.0 

109.0 

2.2 

17.40.57 

.65  4 

.06" 

4.2 

.600 

.153 

5.2 

108.6 

1.1 

17.41.57 

.6.7 

.14® 

4.8 

.541 

.201 

6.8 

106.8 

4.5 

17.42.53 

.63  ' 

.071 

4  .  1 

.577 

.170 

5.6 

107.7 

1.4 

17.43.57 

.67" 

.  0  4  » 

3.5 

.601 

.061 

4.5 

107.0 

.7 

17.44.57 

.  6  P  1 

.  oe° 

4.1 

.59® 

.046 

4.6 

105.9 

•  6 

17.45.57 

.65  = 

.  057 

3.7 

.519 

.156 

6.5 

10*. 8 

.9 

17.46.57 

.  6  j.  " 

.*34 

4.6 

.50? 

.135 

6.7 

104.1 

4.7 

17.47.57 

.64  * 

.10  6 

4.0 

.609 

.142 

6.6 

103.2 

4.7 

17.48.57 

.6  74 

.056 

3.4 

.526 

.197 

7.1 

104.2 

1.2 

17.49. 57 

.  6.11 

. n  7r> 

4.1 

.617 

.0  73 

4.2 

10  5.6 

.5 

17.50.57 

#  ~  1 

.  r  in 

3.8 

.491 

.187 

7.4 

103.0 

.6 

17.51.57 

•  *  *r r 

.077 

4.2 

.604 

.167 

6.9 

102.4 

.7 

17.52.57 

•  A  "  0 

.  0  4  4 

3.8 

.481 

.228 

8.5 

0  102.2 

1.1 

17.55.57 

.6 

,  •  g 

3.9 

.  *116 

.5  5* 

9.9 

10C .  6 

.7 

17.56.57 

.077 

3.8 

.285 

.214 

13.0 

100.1 

.6 

1 7 . 5  8 . 5  7 

.6  7 

.15° 

4 . 7 

.19* 

.  167 

16.1 

98.5 

2-2 

17.  *53. 57 

.66'- 

.*37 

4.1 

.14* 

,09« 

17.3 

100.0 

5.9 

17 . 53 . 5  7 

.67' 

.6  21 

3.4 

.547 

.054 

5.3 

101.6 

.3 

17.54.57 

.  6  7  ’1 

.0  4" 

3.5 

.626 

.117 

5.8 

101.2 

.5 

17.55.57 

*  JC. 

.031 

3.4 

.531 

.140 

6.3 

100.6 

.7 

17.56.57 

.6-0 

.  "3? 

3.6 

.4  84 

.189 

7.6 

100.1 

.5 

17.57.57 

.6:  " 

n  T 
•  <. 

3.6 

.521 

.90? 

7.2 

99.8 

1.0 

17.58.57 

.641 

.1  5 7 

4.3 

.  50° 

.204 

7.4 

98.5 

2.2 

17.59.57 

.6- 

.  n27 

3.6 

.479 

.131 

7.1 

9®. 7 

.4 

18.00.57 

•  A  *' 

•  °  J‘~ 

3.9 

.561 

.  121 

5.3 

9*. 4 

.4 

jjs.ox.*7 

£  « 

4 . 6 

.531 

.154 

6.1 

97.5 

c  •  6 

18. 02.57 

•  ^  ^  J 

;■»  %  *7 

4  .  J 

.  40H 

.  1  99 

9 . 6 

97.6 

.  7 

260 


HATE  0*/nl/7~ 


TImp  f 1 

AV-  S I  a ‘A  "Vi -LOS? 


19.36.5 


19.37.5'’  .a 

19 . 78 . 5T 

'  1  .5 

.? 

br 

A® 

19.39.57 

19.40.57 

.c*-  .?1< 
.5.7  m*y* 

19 . 4 1 , 5  T 

.10-' 

19.44.5 

7  . r  M  .If' 

19.45.5 

19.46.5 

.71'  .7?7 

7  ,  4  >  7  .904 

19.47.57  .'-7- 

1 9 . 48 . 57  .  ®  7  ■  ,9.y« 

19.49.5 

19.50.5 

7  .4;---  •  T  * .  a 

7  .■'>54 

19.51.5 

19.52.5 

7  . .  7  . 1 l;  1. 
•r<t  -71 

19.53.5 

7  .r7'  .'97 

19.55.5 
19. *6. 5 

19.58.5 

19.59.5 
20.00.5 


°P*R«E°  CHA\ 

-|  =  |.  *ATCR 

DFPTh 

VE  ' 

S I C  M  A  3 

-LOS? 

WE 

SJrm'a 

476 

.148 

7.1 

10  0.6 

5.0 

77P 

.924 

11.2 

96 . 0 

9.7 

?  4  7 

.90° 

15.6 

97 . 3 

12.1 

731 

.  ?3® 

13.6  “ 

"99.2 

10-9 

.796 

.938 

10.3 

1 J 1  .  4 

10.1 

40? 

.  ?21 

flT.  3 

1  0  ?  .  7 

6.2 

411 

.938 

10.4 

101.6 

8.1 

46(i 

.  1  7  fl  " 

6.2' 

104.9 

~7J 

35? 

.906 

11.5 

104.3 

4.8 

77  p 

.  ?4  9 

12  '.  0 

95.8 

12.0 

49? 

.14° 

7.0 

100.1 

12.9 

46® 

.175 

6.2 

10?.  9 

9 '.'8 

459 

.17® 

6.3 

1  J4 . 7 

7.6 

100.3  11.7 
10*.*  _1 . 2 
10?.  9  "  1.6 
105.5  9.1 
107.1  9.0 
107.8  8.7 


10.9 


IIP  5 . a 07 


17.21.5T 


17.22.5T 

17.23.57 


17.26.5T 

7.27.5T 


17.28.5T 

17.29.5T 


17.30.57 

17.31.57 


17.34.5'' 
1 7 . 3  5 . 5  T 


17 . 40 . -  7 
17.41.5T 


17.42.5T 


17.48.57 


17.50.5T 
U_.  51. 57 
17  .  s>2 . 5  7 
17.^3.5T 


17.54.57 

17.55. 57 


17.  *8. 5 
17.59.5 
18.00.5 
18.01.5 


18.03.5 


Pir>:Ocr  rHt  HNEL 


SP4RKER  CHANNEL 


WATER  DEPTH 


.OOP  .190 


.  0 4  .144 

.7,^  ,rbo 


.6 

.  * 


.7(,9  .  n  91 

.  7  f  7  .'fir 


'  .  '19r‘ 

2.6 

.  695 

.  I 

.141 

3.2 

.718 

.'S' 

2.8 

.657 

~.l 

L— 

5.3 

.  684 

.  I 

.1^/1 

4.3 

.  72? 

,  1 

r 

3.1 

.  734 

.  I 

100.1 


101.1 


90 

3 

74 

4 

35 

3' 

15 

3 

9 

10  0.3 

.5 

0 

100.0 

.5 

7 

99.6 

.4 

9 

99.6 

.8 

9 

99.2 

.9 

9 

98.4 

2.7 

7 

98.7 

.5 

1 

98.4 

.6 

7 

97.8 

1.7 

3 

97.9 

.5 

3 

96.0 

5.2 

5 

96.6 

1.2 

4 

96.9 

.  8 

2 

96.9 

.8 

1 

94.0 

‘5.4 

4 

95.9 

.5 

6 

95.8 

.  7 

3 

95.4 

.6 

5 

95.3 

.6 

1 

9T .  7 

4.7 

65 

3.2 

9 

73 

3.1 

9 

u 


DATE  0fl/*M/7? 


32110.527 


19 . 48 . 51 
19.49.57 

1.  ’  S'  3.7 

.utn  .  *%■»  ,j 

.9  1  .Tb-'  .0 

19.50. 5"1 

•C'V  .J&9  1.4 
.»tl  .?«■»  1.8 

19.51._V> 

19.52. 5-* 
19.53.51 

rrrs  ror*  yrf - 

•°7«  .1?*  .1 

19.54.51 

19.55.51 

•h-,H  •  r  7  5  1.6 

.731  .-'ll®  1.3 

19.56.51 

19.57.51 

.9in  r5 

.Vj5  1.8 

19.58.51  .  .ljo  1,4" 

-  19.59.5-*  .*7"  .?J3  1.5 

20.00.51 

.AM  .n;t  i  t 

19.26.51 

19.27.5-* 

19.28.5.1 
19. 29. 5"* 

19.30.51 
19. 31. 5"* 

19.32.5"* 

19.13.5i 

pir.r.~°  r,HA\MbL 


SPARKER  CHANNEL 


WATER  DFPTH 


21.21.53 

21.22.5? 

21.23.5? 

21.24.5~* 

21.25.5? 

21.26.5~* 


21.27.5~* 

21.28.5'* 


21.29.5* 

21.30.5? 


21.31.5* 

21.32.5* 


21.33.5* 

21.34.5* 


21.35.5* 

21.36.5* 


21.37.5* 

21.38.5? 


VJ1  o  04  r\jf  O  O  vJiM  O  <51  CD 


I 


DATE  o/"*/7" 


"  "l  L  >‘0 


PI  : °  "h-UHtL 


SPARKER  CHAflNEL 


*  ATFR  DEPTH 


A  \i  - 

SI,'.  H  4 

8-L0S5 

ave 

SIGMA 

b-LDSR 

A  VF 

SIGMA 

20.00.0° 

•  7  v" 

.  <11 

2.8 

.477 

.098 

6.6 

44 . 7 

.7 

20.01.0° 

"7  "1 

•  t 

.  loS 

3.1 

.4  73 

.09? 

6.7 

44 . 7 

.6 

CM 

c 

o 

CM 

• 

J 

.  >3  i 

2.6 

.50  0 

.075 

6.1 

44.6 

.6 

20.03.0" 

.6>;7 

.  123 

3.7 

.4  21 

.027 

7.5 

43.8 

.3 

20.04.00 

.1  *  1 

3.9 

.378 

.073 

a .  6 

43.7 

.6 

20. 05.0" 

.5/4 

.  4 1 

4  .  0 

.  332 

.056 

9.7 

43.7 

.9 

20.06.0° 

.0." 

.345 

1.7 

.40’ 

.063 

8 . 0 

43.3 

.  7 

20.07. j° 

.5*1  < 

M  -  c: 

4.5 

.355 

.164 

6.9 

43.8 

.5 

20.08.0" 

.7/- 

.145 

2.9 

.  329 

.0  78 

9.9 

4’.9 

.  4 

20.10.0° 

.  t  1 

<  A 

•  J 

4.4 

.  435 

.091 

7.4 

47.7 

.9 

20.11.00 

.60" 

.  725 

3.8 

.431 

.068 

7.4 

42.1 

.6 

20.12.Ji 

.1 

.  M" 

3.3 

.42« 

.083 

7.5 

42.4 

1.0 

20.13.0" 

.7.1 

."4  9 

3.6 

.477 

."57 

6 . 5 

41.5 

.5 

20.14.0° 

.6  1 

.’25 

4.1 

.415 

.177 

8.4 

4". 2 

.8 

20. 15.o" 

.4  " 

.  1  09 

8.1 

.  385 

.066 

a .  4 

41.1 

.5 

20.16.0" 

.•»  '4 

.  U7 

8.7 

.501 

.  059 

6.1 

40.9 

.3 

20.17.3" 

.’17 

.*•■  ;7 

8.6 

.484 

.07? 

6.4 

42.3 

5.6 

20.18.0" 

.7-r 

.  1  >9 

8.8 

:  479 

.099 

6 . 6 

47.5 

6*3 

20.19.u-l 
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